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VECTOR AND TENSOR ANALYSIS 



Chapter I 


ELEMENTARY OPERATIONS 

1. Definitions. Quantities which have magnitude only are called 
scalars. The following are examples: mass, distance, area, volume. A 
scalar can be represented by a number with an associated sign, which 
indicates its magnitude to some convenient scale. 

There are quantities which have not only magnitude but also direc¬ 
tion. The following are examples: force, displacement of a point, velocity 
of a point, acceleration of a point. Such quantities are called vectors if 
they obey a certain law of addition set forth in § 2 below. A vector can 
be represented by an arrow. The direction of the arrow indicates the 
direction of the vector, and the length of the arrow indicates the 
magnitude of the vector to some convenient scale. 

Let us consider a vector represented by an arrow running from a 
point P to a point Q,, as shown in Figure 1. The straight line through 
P and dis called the line of action of the vector, the point P is called the 
origin of the vector, and the point Qis called the terminus of the vector. 

To denote a vector we write the letter indicating its origin followed 
by the letter indicating its terminus, and place a bar over the two let¬ 
ters. The vector represented in Figure 1 is then represented by the 
symbols FQ. In this book the superimposed bar will not be used in 
any capacity other than the above, and hence its presence can always 

Q. 
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be interpreted as denoting vector character. This notation for vectors 
is somewhat cumbersome. Hence when convenient we shall use a 
simpler notation which consists in denoting a vector by a single symbol 
in bold-faced type. Thus, the vector in Figure 1 might be denoted by 
the symbol a. In this book no mathematical symbols will be printed 
in bold-faced type except those denoting vectors.* 

The magnitude of a vector is a scalar which is never negative. The 
magnitude of a vector PQ, w iU be denoted by either PQor |PQ,|. Simi¬ 
larly, the magnitude of a vector a will be denoted by either a or |a|. 

Two vectors are said to be equal if they have the same magnitudes 
and the same directions. To denote the equality of two vectors the 
usual sign is employed. Hence, if a and b are equal vectors, we write 

a = b. 

A vector a is said to be equal to zero if its magnitude a is equal to 
zero. Thus a = 0 if a = 0. Such a vector is called a zero vector. 


2. Addition of vectors. In § 1 it was stated that vectors are quanti¬ 
ties with magnitude and direction, and which obey a certain law of 
addition. This law, which is called the law of vector addition , is as follows. 



R 



Let a and b be two vectors, as shown in Figure 2. The origin and 
terminus of a are P and £). A vector equal to b is constructed with 

* It is difficult to write bold-faced symbols on the blackboard or in the exercise 
book. When it is desired to write a single symbol denoting a vector, the reader 
will find it convenient to write the symbol in the ordinary manner, and to place 
a bar over it to indicate vector character. 
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its origin at Q. Its terminus falls at a point R. The sum a+b is the 
vector PR, and we write 

a+b = PR. 


Theorem 1. Vectors satisfy the commutative law of addition; that is, 

a+b = b+a. 


Proof. Let a and b be the two vectors shown in Figure 2. Then 
(2.1) a+b = PR. 


We now construct a vector equal to b, with its origin at P. Its ter¬ 
minus falls at a point S. A vector equal to a is then constructed with 
its origin at S. The terminus of this vector will fall at R, and Figure 
3 results. Hence 



(2.2) b+a = PR. 

From (2.1) and (2.2) it follows that a+b = b+a. 

Theorem 2. Vectors satisfy the associative law of addition; that is, 

(a+b) +c = a+ (b+c). 

Proof. Let us construct the polygon in Figure 4 having the vectors 
a, b, c as consecutive sides. The corners of this polygon are labelled P, 
Q, R and S. It then appears that 

(a+b)+c ='PR+c 

= PS, 

a+(b+c) = a +QS 
= PS. 

Hence the theorem is true. 

According to Theorem 2 the sum of three vectors a, b, and c is 
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independent of the order in which they are added. Hence we can 
write a+b+c without ambiguity. 



Figure 4 shows the construction of the vector a-f-b+c. The sum 
of a larger number of vectors can be constructed similarly. Thus, to 
find the vector a+b+c-fd it is only necessary to construct the 
polygon having a, b, c and d as consecutive sides. The required vector 
is then the vector with its origin at the origin of a, and its terminus 
at the terminus of d. 

3. Multiplication of a vector by a scalar. By definition, if m is a posi¬ 
tive scalar and a is a vector, the expression ma is a vector with mag¬ 
nitude ma and pointing in the same direction as a; and if m is nega¬ 
tive, ma is a vector with magnitude \m \ a, and pointing in the direction 
opposite to a. 

We note in particular that -a is a vector with the same magnitude 
as a but pointing in the direction opposite to a. Figure 5 shows this 
vector, and as further examples of the multiplication of a vector by a 
scalar, the vectors 2a and -2a. 



Theorem. The multiplication of a vector by a scalar satisfies the 
distributive laws; that is, 
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(3.1) (m+n) a = ma+Ka, 

(3.2) m(a+b) = ma+mb. 

Proof of (3.1). If m-\-n is positive, both sides of (3.1) represent 
a vector with magnitude ( m+ri)a and pointing in the same direction 
as a. If m+n is negative, both sides of (3.1) represent a vector with 
magnitude \m + n\a and pointing in the direction opposite to a. 

Proof of (3.2). Let m be positive, and let a, b, ma and mh be as 
shown in Figures 6 and 7. Then 



(3.3) m(a-j-b) = mPR, ma+mb = SU. 

The two triangles PQR and STU are similar. Corresponding sides 
are then proportional, the constant of proportionality being m. Thus 

(3.4) mPR=SU. 

Since PR and SU have the same directions, and since m is positive, 
then mPR = SU. Substitution in both sides of this equation from 
(3.3) yields (3.2). 

Now, let m be negative. Then Figure 7 is replaced by Figure 8. 
Equations (3.3) apply in this case also. The triangles PQR and STU 
are again similar, but the constant of proportionality is \m\, so \m\PR 
= SU. Since PR and SU have opposite directions and m is negative. 
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Figure 8 

then mPR = SU. Substitution in both sides of this equation from 
(3.3) again yields (3.2). 

4. Subtraction of vectors. If a and b are two vectors, their difference 
a - b is defined by the relation 

a-b = a+(-b), 

where the vector -b is as defined in the previous section. Figure 9 
shows two vectors a and b, and also their difference a-b. 



Figure 9 


5. Linear functions. If a and b are any two vectors, and m and n 
are any two scalars, the expression ma-f »b is called a linear function 
of a and b. Similarly, ma + nb +pc is a linear function of a, b, 
and c. The extension of this to the cases involving more than three 
vectors follows the obvious lines. 

Theorem 1. If a and b are any two nonparallel vectors in a plane, 
and if c is any third vector in the plane of a and b, then c can be 
expressed as a linear function of a and b. 

Proof. Since a and b are not parallel, there exists a parallelogram 
with c as its diagonal and with edges parallel to a and b. Figure 10 
shows this parallelogram. We note from this figure that 
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(5.1) c = PQ.+QP. 

But PQ is parallel to a, and QR is parallel to b. Thus there exist 
scalars m and n such that 

PCI = ma, QR = nh. 

Substitution from these relations in (5.1) yields 

c = ma+nb. 

Theorem 2. If a, b and c are any three vectors not all parallel to 
a single plane, and if d is any other vector, then d can be expressed 
as a linear function of a, b and c. 

Proof. This theorem is the extension of Theorem 1 to space. Since 
a, b and c are not parallel to a single plane, there exists a parallele¬ 
piped with d as its diagonal and with edges parallel to a, b and c. 
Hence there exist scalars m, n and p such that 

d = m’a.-pnh-ppc. 

6. Rectangular cartesian coordinates. In much of the theory and appli¬ 
cation of vectors it is convenient to introduce a set of rectangular 
cartesian coordinates. We shall not denote these by the usual symbols 
x,y and z, however, but shall use instead the symbols x t , x 2 and x 3 . 
These coordinates are said to have “right-handed orientation” or to be 
“right-handed” if when the thumb of the right hand is made to point 
in the direction of the positive x 3 axis, the fingers point in the direction 
of the 90° rotation which carries the positive x 1 axis into coincidence 
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with the positive x 2 axis. Otherwise the coordinates are “left-handed”. 
In Vector Analysis it is highly desirable to use the same orientation 
always, for certain basic formulas are changed by a change in orien¬ 
tation. In this book we shall follow the usual practise of using right- 
handed coordinates throughout. Figure 11 contains the axes of such 
a set of coordinates. 



It is also convenient to introduce three vectors of unit magnitude, 
one pointing in the direction of each of the three positive coordinate 
axes. These vectors are denoted by i x , ig and i 3 , and are shown in 
Figure 11. 

Let us consider a vector a. It has orthogonal projections in the 
directions of the positive coordinate axes. These are denoted by a x , a 2 
and a 3 , as shown in Figure 11. They are called the components of a. 
It should be noted that they can be positive or negative. Thus, for 
example, a x is positive when the angle between a and the direction of 
the positive x x axis (the angle QPS in the figure) is acute, and is nega¬ 
tive when this angle is obtuse. 

From Figure 11 it also appears that a is the diagonal of a rectangular 
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parallelepiped whose edges have lengths |< 2 X |, \a 2 \ and ja 3 |. Hence the 
magnitude a of the vector a is given by the relation 

(6.1) a = Va^+a^+a^. 

From the figure it also appears that 

(6.2) a = pq+qr+rs. 

Now the vector PQ is parallel to i v Because of the definitions of 
and of the product of a scalar by a vector, we then have the 1 elation 
PQ = <Zjii. Similarly QR = a 2 i 2 and RS = a^. Substitution in (6.2) 
from these relations yields 

(6.3) a = . 

This relation expresses the vector a as a linear function of the unit 
vectors , i 2 and i 3 . We note that the coefficients are the components 
of a. 

Theorem. The components of the sum of a number of vectors are 
equal to the sums of the components of the vectors. 

Proof. We consider two vectors a and b with components , a z , 
b 1 , b 2 and b 3 . Then 

a = ++I 2 4- a 3 i 3 , 

b = +A 2 i 2 A 3 i 3 . 

Addition of both sides of these equations leads to the relation 

a+b = +a 2 i 2 ++ ^ 1 ++ b.^ + b 3 i 3 . 

Now the sum of a number of vectors is independent of the order in 
which the vectors are added, by Theorem 1 of § 2. Hence we may 
write the above equation in the form 

a+b = tfjix+ ipi+ <2 2 i 2 + b 2 i 2 +a 3 i 3 + +i 3 . 

By the theorem in § 3 we may then write this in the form 
a+b = (a. 1 -\-b^)i 1 -\-(a 2 -\-b^)i 2 -\-(a 3 -\-'b^)i 3 . 

Hence the components of a+b are a 3 +b 3 , a 2 -\-b 2 and a 3 -\-b 3 . This 
proves the theorem when two vectors are added. The proof is similar 
when more than two vectors are added. 
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7. The scalar product. Let us consider two vectors a and b with 
magnitudes a and b, respectively. Let oc be the smallest nonnegative 
angle between a and b. as shown in Figure 12. Then 0° < a < 180°. 



The scalar ab cos a arises quite frequently, and hence it is convenient 
to give it a name. It is called the scalar product of a and b. It is also 
denoted by the symbols a • b, and hence we have 

( 7 . 1 ) ab = ab cos a. 

The scalar product Is sometimes referred to as the dot product. 

If the components of a and b are denoted by , a 2 , a 3 , b ls b 2 and b 3 
in the usual manner, the direction cosines of the directions of a and b 
are respectively 

a-± a% a 3 b ^ b 3 b 3 

T 7 ’ 7 ’ J’ T T 


By a formula of analytic geometry, we then have 


cos a 


«1 ^2_|_^3 ^3 

a b a b a b 


Substitution in (7.1) of this expression for cos a yields 
(7.2) a-b = a x b 1 -\-a 2 b 2 -\-a z b 3 . 

This relation expresses the scalar product of two vectors in terms of 
the components of the vectors. 

Theorem 1. The scalar product is commutative; that is, 


a-b = b ■ a. 


Proof. Because of (7.2), we have 

a-b = a 1 b 1 -\-a 2 b 2 -\-a 3 b 3 , 
b-a = b l a 1 -(- b 2 a 2 T b 3 a 3 . 
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Since a 1 b 1 — , etc., the truth of the theorem follows immediately. 

Theorem 2. The scalar product is distributive; that is, 

a- (b+c) = a-b+a-c. 

Proof. According to the theorem in § 6, the components of b+c 
are ii+c ls b 2 J r c i and b 3 -\-c 3 . Hence, by (7.2) we have 

a-(b + c) = a^b^c^) -\-a 2 (b 2 -\-c^) -\-a 3 [b 3 -\-c^) 

= aib 1 -\-a 2 b 2 J r a 3 b 3 J ra 1 c 1 -{-a 2 c 2 -\-a 3 c 3 

= ab+ac. 

This completes the proof. 

If a and b are perpendicular, then 

a-b = 0. 

However, if it is given that a • b = 0, it does not necessarily follow 
that a is perpendicular to b. It can be said only that at least one of the 
following must be true :a=0;6=0;ais perpendicular to b. Similar¬ 
ly, if it is given that 

a-b = a-c, 

it does not necessarily follow that b = c. For this relation can be 
written in the form a ■ (b - c) =0, and hence it can be said only that 
at least one of the following is true :a=0;b=c;ais perpendicular 
to the vector b-c. 

We note the following expressions, in which a is any vector and 
i x , L and i 3 are the unit vectors introduced in § 6: 

a-a = a 2 , 

ii-ii = 1, i 1 -i 2 = 0, i 1 -i 3 =0, 

(7.3) i 2 -i a =0, 12-12 = 1, V i 3 = 0, 

V»i = °> V h = 0, i 3 i 3 = 1- 

8 . The vector product. Let us again consider two vectors a and b, 
the smallest nonnegative angle between then being denoted by <x, as 
shown in Figure 12. Then 0° < a < 180°. The vector product of a 
and b is a third vector c defined in terms of a and b by the following 
three conditions: 
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(i) c is perpendicular to both a and b; 

(ii) the direction of c is that indicated by the thumb of the right 
hand when the fingers point in the sense of the rotation a from 
the direction of a to the direction of b; 

(iii) c = ab sin a. 

These conditions define c uniquely. Figure 13 shows c. The vector 



Figure 13 


product of a and b is also denoted by a X b. Hence 

(8.1) c = aXb. 

The vector product is also called the cross product. 

Theorem 1. The area A " of the parallelogram with the vectors a 
and b forming adjacent edges is given by the relation 

(8.2) A = |aXb|. 


Proof. Figure 14 shows the parallelogram. If p is the perpendicular 
distance from the terminus of b to the line of action of a, then A = ap. 
But p = b sin a. Hence 

A = ab sin a 
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Figure 14 


We shall now determine the components of the vector product c 
in (8.1) in terms of the components of a and b. Because of condition (i) 
above, we have ac = 0, b e =0. Because of (7.2), these equations 
can take the form 

a \ c \~\~ a ip2~\~ iVs = 0, 

T b 2 C 2 “t - ^3^3 = 

If these equations are solved for c x and c 2 in terms of c 3 , it is found 
that 

c i = ^2 = _ g 3 

^2^3 ~ fl 3^2 ~ a A &\b 2 ^2^1 

In order to preserve symmetry, we denote the common value o 
these three fractions by K , whence we have 

Cj = K[a 2 b 3 —a 3 b 2 ) 1 
(8.3) c 2 — K[a 3 b±—, 

c 3 = K(a 1 b 2 - a 2 bj). 

Now c 2 = c 1 2 +r 2 2 + f 3 2 - Hence 

c 2 = K 2 [{a 2 b 3 -a 3 b 2 ) 2 +{a 3 b x -a 1 b 3 ) 2 +( K a 1 b 2 -a 2 b^ 2 '] 

= K 2 [a 1 2 (b 2 2 +b 3 2 )+a 2 2 (b 3 2 + b 1 2 ) + a 3 \b 1 2 +b 2 2 ) 

—2 (a 2 b 2 a 3 b 3 -f- a 3 b 3 a 1 b 1 + a 1 ^ 1 a 2 ^ 2 )] • 

The first term inside the square brackets can be written in the form 
a 2 {b 2 -b 2 ). If the second and third terms are treated similarly it is 
found that 

c 2 = K 2 [{a 2 +a 2 +a 3 2 ) b 2 ~(a 1 b 1 +a 2 b 2 +a 3 b 3 ) 2 ] 

= K 2 [a 2 b 2 -(ab cos a) 2 ] 

= K 2 a 2 b\ 1-cos 2 a) 

= K 2 a 2 b 2 sin 2 a. 

But by condition (iii) above, c 2 = a 2 b 2 sin 2 «. Thus K = ± 1. If 
these two values of K are inserted in (8.3) two vectors c result with 
the same magnitude but pointing in opposite directions. Only one of 
these vectors satisfies condition (ii) above. Now both values of K are 
numerical, and are hence independent of a and b. Thus the same 
value of K will satisfy condition (ii) for all vectors a and b. Hence it 
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is only necessary to find K for any one special case in which c can be 
found directly and with ease from conditions (i) - (iii) above. If we 
take a = i x and b = i 2 , it is found from these conditions that c = ig. 
Thus (ii = b 2 = c 3 = 1 , a 2 = a 3 = b 1 = b 3 = = c 2 = 0, and substi¬ 

tution in (8.3) yields K — 1. From (8.3) we then have in general, 

( 8 -f) — a 2 b 3 -a 3 b 2 , c 2 — a 3 by~dyb 3 , c 3 — o 3 b 2 o 2 by* 

Thus the components of the vector product axb are a 2 b 3 -a 3 b 2 , o 3 b 1 —a 1 b 3> 
a 1 b 2 -a 2 b 1 . Hence 

(8.5) axb = {o- 2 b 3 —d 3 b^)\i-\-(^ii 3 b- l — df>^)\ 2 -\-^o,- i b 2 —a. 2 b-P)i 3 , 
or, in determinant form 

»i h h 

(8.6) axb = flj d 2 a 3 . 

by b 2 b 3 

Theorem 2. The vector product is not commutative, because 

(8.7) axb = -bxa. 

Proof. By ( 8 . 6 ) it follows that 

*i *2 i 3 

bxa = by b 2 b 3 . 

dy Og d 3 

Since this determinant differs from the determinant in ( 8 . 6 ) only 
in that two rows are interchanged, the two determinants differ only 
in sign. Hence (8.7) is true. The truth of this theorem can also be 
seen easily by examining the three conditions which define the vector 
product. According to these conditions the effect of interchanging the 
order of a and b is only to reverse the direction of the vector product. 
Theorem 3. The vector product is distributive; that is, 

( 8 . 8 ) ax(b+c) =axb+axc. 

Proof. Let us write 

d=ax(b+c), e = axb, f = axc. 
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Then 


d \ — ^2(^3 I ^3) _ ^3 (^2 d - ^2) 

= (^2^3 “ ^3^2) "b (^2^3 ~ ^3^2) 

= %+/i- 

Similarly d 2 = e 2 -\-f 2 and d 3 = « 3 +Ja • Hence d = e-ff, and so ( 8 . 8 ) 
is true. 

We note the following expressions, in which a is any vector and i t> 
i 2 and 13 are the unit vectors introduced in § 6 : 



a X a = 0, 


*iX>i — 0, 

*i XI 2 = 13 3 

liXia = -i 2 , 

.r- 

X 

►i'¬ 

ll 

cT- 

o' 

II 

X 

CQ 

• H 

*2 ^ *3 = *1 3 

hX'l =i 2, 

i 3 Xi 2 = -i 1; 

*3 X^3 — 0. 


If a and b are parallel, then 

axb = 0 . 

Also, if it is given that a X b =0, then at least one of the following 
must be true :a=0;6=0;ais parallel to b. Similarly, if 

a X b = a X C, 

then aX (b-c) =0 and at least one of the following must be true: 
a = 0 ; b=c; a is parallel to b-c. 


9. Multiple products of vectors. Let a, b and c be any three vectors. 
The expression 

a • (b X c) 

is a scalar, and is called a scalar triple product of a, b and c. 

If the components of a, b and c are denoted in the usual way, then 
the components of b X c are b 2 c 3 - b 3 c 2 , b 3 c x - b x c 3 , b^-b ^, and we have 

by (7.2) 

a • (b X c) = ^3^2) — ^1^3) H~ ^3 (. ^1^2 


or 

(9.1) 


a- (bxc) 


^2 ^3 

&i b% 

C 1 C 2 C 3 


15 



Theorem 1. The permutation theorem for scalar triple products. If the 
vectors in a scalar triple product are subjected to an odd number of 
permutations, the value of this product is changed only in sign; and 
if the number of permutations is even the value of the product is not 
changed. 

Proof. A permutation of the vectors in a scalar triple product is 
defined as the interchange of any two vectors which appear in the 
product. From (9.1) it appears that a single permutation produces an 
interchange of two rows in the determinant. Since such an interchange 
of rows results in a change of sign only, the truth of the theorem is 
established. 

Because of this theorem we have 

a-(bxc) =b-(cxa) =c-(axb) 

= -c ■ (b X a) = -a- (cxb) =-b-(aXc). 

Theorem 2. The volume V of the parallelepiped with the vectors a, 
b and c forming adjacent edges is given by the relation 

(9.2) V = |a- (b x c) |, 

where the vertical lines here denote the absolute value. 

Proof. Figure 15 shows the parallelepiped. Let d =bxc. Then 



d = |bxc|, and by Theorem 1 of § 8 the area of the parallelogram 
forming the base of the parallelepiped is then d. Hence V = hd, 
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where h is the altitude of the parallelepiped. But d is perpendicular 
to the base, and if /? is the angle between a and d, then h = a | cos /? |. 
(The absolute value signs are necessary here, since /? lies in the range 
0° < /? < 180° and hence cos ft may be negative.) Thus 

V = I ad cos /? | 

= I ad | 

= | a- (bxc)|. 

The expression 

ax(bxc) 

is a vector, and is called a vector triple product of a, b and c. Let us write 

d = b X c, e=axd. 


Then e is equal to the vector triple product ax(bxc). By (8.5) 
we have 

Cl 3 d 3 Q 3 d 2 

= a 2 (b i c 2 — b 2 c x ) — o. 3 {b 3 c 3 — b x c 3 ) 

= bp d 2 c 2 -p ci 3 c 3 j — ^i(^2^2 “f ^ 3 ^ 3 ) • 


Because of (7.2), this can be written in the form 

= Ma-c-qcJ-qfa'b-aA) 

= ^i(a-c) -c x (a-b). 

Similarly 

«2 = i 2 (a-c)-c 2 (a-b) 3 
= b 3 (a.-c)-c 3 (n-h). 

Hence e — b(a-c)-c(a-b), and since e = ax(bxc), we have 
(9.3) ax(bxc) = b(a-c)—c(a-b). 

This is a rather important identity. It will be used frequently. 

We note that the right side of (9.3) is a vector in the plane of b 
and c. This is to be expected, since the vector a X (b X c) is perpen¬ 
dicular to the vector b X c which is itself perpendicular to the plane 
of b and c. 

Let us now consider the expression 

(axb) X (cxd). 
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It is a vector. If we regard it as a vector triple product of a X b, 
c and d, then by (9.3), 

(9.4) (axb)x(cxd) = c[(axb)-d]-d[(axb)-c]. 

Since an interchange of the order of the vectors in a vector product 
only changes the sign, 

(axb)x(cxd) =-(cxd)x(axb). 

If we regard the right side of this equation as the vector triple pro¬ 
duct of c X d, a and b, then by (9.3), 

(9.5) (axb)x(cxd) =-a [(cxd)-b]+b [(cxd)-a], 

We next consider the expression 

(axb) • (c xd). 

It is a scalar. If we consider it as the scalar triple product of a X b, 
c and d, and subject these three vectors to two permutations, then 
according to Theorem 1 of § 9, we have 

(axb)-(cxd) = c-[d X (axb)]. 

If the vector triple product on the right-hand side of this equation is 
expanded by the identity in (9.3), we obtain 

(axb) - (c xd) = (c• a)(d■ b)—(c• b)(d• a), 

or in a form more easily recalled, 

(9.6) (axb)-(cxd) = (a-c)(b-d)-(b-c)(a-d). 

There are many other multiple products of vectors. In general these 
can be simplified by means of the theorems and formulas above. In 
dealing with multiple products of vectors care must be exercised to 
avoid writing down expressions which are ambiguous or have not 
been defined. Thus, for example, the expression axbxc is ambi¬ 
guous since ax(bxc) ^ (axb)xc, and the following expressions 
have not been defined: 

ab, a • (b • c), a X (b • c). 

10. Moment of a vector about a point. Let a be a vector with origin at 
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Figure 16 

a point X, and let x be a vector with origin at a point 0 and terminus 
at X, as shown in Figure 16. The moment of a about the point 0 is by 
definition the vector P given by the relation 

P = xXa. 

Theorem 1. If P is the moment of a about a point 0, then 

P —pa, 

where p is the perpendicular distance from 0 to the line of action of a. 

Proof. Now P — xa sin a, where a is the angle between x and a. 
But p = x sin a. Hence P = pa. 

Theorem 2. The moments about a point of any two equal vectors 
with the same line of action are equal. 

Proof. Let a and a’ be two equal vectors with the same line of ac- 
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Figure 17 
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tion, and with origins X and X', as shown in Figure 17. Let 0 be any 
point, and let OX = x, OX' = x'. Let P and P' be the moments of 
a and a' about 0. Then 

P = xXa, P' =x'xa'. 

But x' = x-f-XX". Hence 

P' = (x+ZZ')xa' 

= xXa'+S'Xa'. 

Since XX' is parallel to a', XX' X a' = 0. Hence, since a' = a we 
have finally 

P' = x X a = P . 

11. Moment of a vector about a directed line. Each line defines two 
directions which are opposite. A line is said to be directed when one of 
these directions is labelled the positive direction and the other the nega¬ 
tive direction. 

Let us consider a directed line L, and let b denote a unit vector 
pointing in the positive direction of the line, as shown in Figure 18. We 



also introduce a point 0 on L and a vector a with origin at a point X. 
If P denotes the moment of a about 0, the moment of a about L is by 
definition the orthogonal projection of P on L. It is a scalar, and if it is 
denoted by Q,we have 

Q, = P cos ®, 
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where fi> is the angle between P and the unit vector b. Hence Q, = b • P. 
But P = x X a, where x = OR, Thus 

(H.l) & = b.(xxa). 

Theorem 1. The above definition of the moment of a vector about 
a directed line L is independent of the position of the point 0 on L. 
Proof. Let O' be a second point on L , as shown in Figure 19. Also, 



let P' be the moment of a about O', and let Q be the corresponding 
moment of a about L. Then 

Q,= b-(xxa) , (y = b-(x'Xa) , 
where x' is as shown. But x' = 0"0+x. Thus 

£L' = b [(0 77 0+x) xa] 

= b-(0' r 0Xa)+Q,. 

Since b and 0"0 have the same line of action L, then b • (0"0 Xx) =0 
by Theorem 2 of § 9. Thus Ql = Q, which proves the theorem. 

Theorem 2. If P denotes the moment of a vector a about the origin 
of the coordinates, then the three components of P are equal respec¬ 
tively to the moments of a about the three coordinate axes. 

Proof. The truth of this theorem follows immediately from the 
above definitions of the moments of a vector about a point and about 
a line. 
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12. Differentiation with respect to a scalar variable. Let u be a scalar 
variable. If there is a value of a vector a corresponding to each value 
of the scalar u, a is said to be a function of u. When it is desired to in¬ 
dicate such a correspondence, we write a (u). 

Let us consider a general value of the scalar u and the corresponding 
vector a(u). Let the vector OP in Figure 20 denote this vector. We 



Figure 20 

now increase the scalar u by an amount A u. The vector corresponding 
to the scalar h+Ah is a(H+Au). Let the vector 0Q,in Figure 20 denote 
this vector. The change in a(u) corresponding to the change A u in u 
is then a(n + An) — a(a). In the usual notation of calculus we denote 
it by Aa, so that 

Aa = a(a +Ah) - a(u). 

From the figure it is seen that Aa = PQ_. Since Ah is a scalar, the 

^ 

vector —- has the same direction as PO . The vector 
A u 

limit 
a u->o Am 

is the rate of change of a with respect to u. It is also called the deriva- 

fj *| 

tive of a with respect to u, and is denoted by the symbol —, so that 



In precisely the same way, we define the derivative with respect to a 

d.2L 

of the vector - 7 - . This vector is denoted by 
du 


d [da. d 2 a 

du du ^ ° r du 2 


Higher derivatives of a with respect to u are defined similarly. 

Let a(a) and b(a) be any two vectors which are functions of a scalar 
u, and let m be a scalar function of u. We shall now derive the following 
formulas: 


( 12 . 1 ) 

( 12 . 2 ) 


d 1 1 \ da. , d b 

& (a+b) = da + ST 

d . . da. , dm 

*("•) =»*■+**• 


(12.3) 



dh da 

du du 


b, 


(12.4) 


d_ 

du 


X b) = a X 


dh , da , 

da + da X 


Proof of (12.1). When a increases by an amount Aa, the change in 
the sum a+b is 

(12.5) A(a+b) = (a+Aa+b + Ab) - (a+b). 

According to Theorem 2 of § 2 the sum of a number of vectors is inde¬ 
pendent of the order of summation. Thus the right side of (12.5) can 
be written in the form a - a+b - b+Aa+Ab, which reduces to 
Aa+Ab. Thus 

A(a+b) = Aa+Ab. 

If both sides of this equation are divided by Aa, and if Aa is then made 
to approach zero, ( 12 . 1 ) is obtained. 

Proof of (12.2). When a increases by an amount Aa, the change in 
ma is 

(12.6) A(ma) = (m+Am)(a+Aa) — ma. 

According to the theorem in § 3 the multiplication of a vector by a 
scalar satisfies the distributive laws, as exemplified by Equations (3.1) 
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and (3.2). Because of the law exemplified by (3.1) we can then write 

(12.6) in the form 

(12.7) A (ma) = m(a-f-Aa)+Am(a + Aa) -ma, 

and because of the law exemplified by (3.2), we can then write (12.7) 
in the form 

A(ma) = ma+mAa+ Am a + Am Aa - ma 
= m Aa+Am a+Am Aa. 

If both sides of this equation are divided by A u, and if A u is then made 
to approach zero, (12.2) results. 

Proof of (12.3). When u increases by an amount A u, the change in 

a b is 

A(a-b) = (a+Aa) • (b + Ab) - a-b. 

Since the scalar product is distributive, this equation may be written 
in the form 

A(a-b) = a-b + a- Ab + Aa-b+Aa-Ab - a-b 
= a-Ab+Aa-b+Aa-Ab. 


If both sides of this equation are divided by the scalar A u, we have 

=a .^+^. b+ ^.Ab. 

A u A u ' Au 'au 

If we now let A« approach zero, (12.3) results. 

Proof of (12.4). This proof follows exactly the same pattern as the 
proof of (12.3), and hence will not be given here. 

It is important to note that in (12.4) the order in which the vectors 
a and b appear must be the same in all terms, since a x b is not equal 
to b X a. 

If a(a) is a vector with components a u a 2 and a 3 , then 


a — apj+agig+ajig. 
By (12.1) and (12.2) we then have 

da. d . • , , d . • . - d . .. 
^ = Tu^ )+ Tu M 


_ n dii 

du 




di 2 

du 


d*z j d&i . dd 2 . ( d&2 • 


+ U.M.1 • M.U-2 • , U 

w ■‘ + * ■*+&■■ 
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Now a 1 , a 2 and a 3 are scalar functions of u. Also i l5 i 2 and i 3 are unit 
vectors pointing in the directions of the positive coordinate axis. If they 
are the same for all values of u, then 

_ dh _ dh n 

du du du ’ 

and so 

fk _ da \ • j_da 2 . da 3 . . 
du du 1 du 2 du * 3 

From this equation we see that the components of the derivative of a 
vector are equal to the derivatives of the components, provided the 
directions of the coordinate axes are independent of the variable of 
differentiation. 


13. Integration with respect to a scalar variable. ' Let a be a given func¬ 
tion of a scalar u. We introduce orthogonal unit vectors i x , i 2 and i 3 
with directions independent of u. Then 

a(u) = 

We make the definition 

(13.1) / a (“) d u = hf a ii u ) du+i 2 Ja 2 (u) du+i 3 Ja 3 (u) du. 

From this definition it follows that 

d r 

— I a (a) du = %(«) +i 2 a 2 (a) +i 3 a 3 (a) = a, 

as expected, since integration is the inverse of differentiation. It is to be 
noted that each integral on the right side of (13.1) gives rise to a con¬ 
stant of integration. 

Theorem. If a (a) is a linear function of constant vectors, with coef¬ 
ficients which are functions of u, then J a(a) du can be obtained by for¬ 
mal integration in which constant vectors are treated as are constants 
in ordinary integration, and arbitrary constant vectors are inserted 
where arbitrary constants would appear in ordinary integration. 
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Proof. We have 


a (“) = p/(«)+qg(«H— > 

where p, q, are constant vectors and /(a), g(u), ••• are given 
functions of a. By (13.1) it then follows that 

/ a (“) du =% [Piff(u) du+q 1 Jg{u ) du+ ■ • •] 

//(«) -]+» 3 [As//(«) du -1-] ' 

= P//(a) <&+q/£(«) -. 

If k, l, • ■ ■ denote the integration constants of the integrals in the last 
line, then the total contribution of these constants to J a (a) du is the 
single arbitrary constant vector c such that 

c = p£+q/+ • • ■ . 

The following examples illustrate the above theorem: 

J (pu-\-q)du = |pa 2 +qa-|-c, 
f p cos a rfa = p sin a+c. 


14. Linear vector differential equations. The equation 


(14.1) 


(A> 


i"' 1 
~^ 1 ia ”" 1 


+A-i ^+A) x = a, 


in which a and po, pi, • • -,p„ are given functions of the scalar a and x 
is an unknown vector, is a linear vector differential equation of order n. 
Let F denote the differential operator in (14.1). Then (14.1) can be 
expressed in the form 

(14.2) .F[x]=a. 


Theorem. The general solution of the linear vector differential 
equation F[x] = a is x = Y+A, where A is a particular solution of 
this differential equation, and 

Y = CjA+c^+CsjgH-l-c n y n , 
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C 13 c 2 , c 3 , being arbitrary constant vectors andy 1} y 2 , y 3) ■ ■ - ,y„ being 

n linearly independent solutions of the homogeneous scalar differential 
equation F[jt] = 0 . 

Proof. Let us introduce the unit vectors q, i 2 and i 3 with directions 
independent of a. Then 

a = cpj+fljij+fljig, 

X ^llx I H - ^3*3 , 

F[x] =FMi 1 +%]i 2 +%]i 3 . 

Hence, from (14.2) we have 

(14.3) F\ x^=ay, F[x 2 \ = a 2 , F[> 3 ] = a 3 . 

Let , A 2 and A 3 denote particular solutions of these three scalar dif¬ 
ferential equations, and lcty 1 ,y 2 , ■ ■ ■ ,y„ denote u linearly independent 
particular solutions of the scalar differential equation F[j] =0. Then 
the general solutions of Equations (14.3) are 

*i= %JVi+<a 2 Ja+WaH- 

*2= C 2lJ'l + C 22j'2 + C 23J>'3d-- + C 2nJn+^2> 

*3= %J>'l + (r 32J>’2+^33J , ’3H-- “H <r 3nJ'n+^3 > 

where the c’s are arbitrary constants. Let us multiply these three equa¬ 
tions by i l5 i 2 and i 3 , respectively, and then add. The result can be 
written in the form 

(14.4) x=Y+A, 
where 

Y = c 1 j 1 +c 2 j 2 +c 3 y 3 H-he „y n , 

A = T 1 i 1 +T 2 i 2 +T 3 i 3J 

the vectors c l3 c 2 , c 3 , being arbitrary constant vectors. Equa¬ 

tion (14.4) gives the general solution of Equation (14.1). We note that 
Y is the general solution of the homogeneous equation F[x] = 0, and 
that A is a particular solution of Equation (14.1). The particular solu¬ 
tion A can be found by procedures very similar to those used to find 
particular solutions of linear scalar differential equations. This is de¬ 
monstrated below. 
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As an example, let us find the general solution of the differential 
equation 


(14.5) 


d 2 x 


dx 

du 


-2x = lOp sin a+q(2a + l), 


where p and q are constant vectors. We must first find two linearly 
independent solutions of the equation 


(14.6) 


d 2 y dy 
du 2 du 


- 2y = 0. 


The auxiliary equation of this differential equation is 

m? - m-2 =0. 


It has roots -1, 2, whence the required solutions of (14.6) are e'“ and 
e 2 “. Thus 

Y = <v-“+<v 2 “. 

We now use the method of undetermined coefficients to find a particu¬ 
lar solution A of Equation (14.5). The function on the right side of 
(14.5), and the derivatives of this function, are linear functions of sin u, 
cos a, u, 1, none of which are particular solutions of (14.6). Hence, we 
look for a particular solution A in the form 

A = b sin a+c cos a+da+e, 

where b, c, d and e are constant vectors. By substitution in (14.5) we 
readily find by equating coefficients that 

3b-c.= -10p, b + 3c = 0, 
d = q, d+2e = -q. 

Solving these four equations for b, c, d, e, we find that 
b = -3p, c = p, d = -q, e = 0. 

The general solution of Equation (14.5) is then 

x = c 1 e'“-)-c 2 e 2 “-|-p(-3 sin a-f-cos a) -qa. 


Problems 

1. The vectors a, b, c and d all lie in a horizontal plane. Their 
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magnitudes are 1, 2, 3 and 2, and their directions are east, northeast, 
north and northwest, respectively. Construct these vectors. 

2. If a, b and c are defined as in Problem 1, construct the vectors 
(a+b)+c, (b+a)+c, c + (a+b), and by measuring their magnitudes 
and directions verify that they are equal. 

3. Ifa and b are defined as in Problem 1, contruct the vectors a+2b, 
2a+b, 3a-b, -2a-2b. 

4. If a, b and c are defined as in Problem 1, express each of these 
vectors as a linear function of the other two, determining the coeffi¬ 
cients graphically to two decimal places in each case. 

5. Given that 

a 2b = m, 2a b = n, 

where m and n are known vectors, solve for a and b. 

6 . If a and b are vectors with a common origin 0 and terminuses 
A and B, in terms of a and b find the vector OC, where C is the middle 
point of AB. 

7. The vectors a and b form consecutive sides of a regular hexagon, 
the terminus of a coinciding with the origin of b. In terms of a and b 
find the vectors forming the other four sides. 

8 . If a, b, c and d have a common origin and terminuses A, B, C 
and D, and if b-a = c-d, show that ABCD is a parallelogram. 

9. The vectors a, b and c have a common origin and form adjacent 
edges of a parallelepiped. Show that a+b+c forms a diagonal. 

10. Vectors are drawn from the center of a regular pentagon to its 
vertices. Show that their sum is zero. 

11. Consider the vectors a, b, c and d defined in Problem 1. Intro¬ 
duce rectangular cartesian coordinate axes such that the four vectors 
lie in the x x x % plane with the % axis pointing east and the x 2 axis point¬ 
ing north. Find the components of the vectors a, b, c, d, a+2b, and 
3a-2b; also, express these vectors in terms of their components and 
the unit vectors i x , i 2 and i 3 . 

12. Do Problem 4, making an exact determination of the coefficients 
analytically by the use of components. 
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13. Given that 


a = i 1 +2i 2 +i 3 , 
b = 2i 1 +i 2) 
c = 3i x - 4i 2 - 5i 3 , 

verify that a-(b+c) = ab+ac, ax(b+c) =axb+aXc. Also, 
find |axb|, |axb+aXc|. 

14. Show that (a+b) • (a - b) =a 2 -b 2 . 

15. Show that (a-b) X (a+b) =2axb. 

16. The two vectors a = i 1 +3i 2 -2i 3 , b = 3i 1 +2i 2 -2i 3 have a 
common origin. Show that the line joining their terminuses is parallel 
to the x x x 2 plane, and find its length. 

17. Show that the vectors a = i 1 +4i 2 +3i 3 , b = 4i x +2i 2 -4i 3 are 
perpendicular. 

18. If a, b and c are as defined in Problem 13, find a-(bxc), 
(bxa)-.c, ax(bxc), (axb)-(axc), (axb)x(aXc). 

19. If the vectors drawn from the origin to three points A, B and C 
are respectively equal to the three vectors a, b and c defined in Prob¬ 
lem 13, find a unit vector n perpendicular to the plane ABC. Hence 
find the distance from the origin to this plane. 

20. Show that 

ax (bxc)+bx (cXa)+cX (axb) =0. 

21. Show that 

ax[bx (cxd)] = [a-(cXd)]b-(a-b) (cxd) 

= (b-d) (axc)-(b-c) (axd). 

22. Show that 

[a X b] ■ [(b X c) X (c X a)] = [a • (b X c)] 2 . 

23. Show that 

a[b- (c xd)] - b[c- (dxa)]+c[d- (axb)] - d[a- (b X c)] = 0. 

24. Show that 

a b c 

[a-(bxc)](fxg) = f-a f-b f-c . 

g a g b g c 
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25. Show that 

n-(bxc) n-(cXa) n-(axb) 

n =-a-|-b-|-c . 

a-(bxc) b-(cxa) c-(axb) 

This formula can be used to express any vector n as a linear function of 
any three vectors a, b and c not lying in the same plane. 

26. Express the vector n = i 1 H-2i 2 -)-3i3 as a linear function of the 
vectors a, b and c defined in Problem 13. (See Problem 25.) 

27. Express the vector n = 2i x - 2i 2 - 3i 3 as a linear function of the 
vectors a, b and c defined in Problem 13. (See Problem 25.) 

28. Show that 

a X [(f X b) X (g X c)] +b X[(f X c) X (g X a)] 

+cx[(fxa) X(gXb)] =0. 

29. If 0 is the origin of the coordinates and A, B and C are three 
points such that 

OA = 2i 1 +2i 2 -i 3 , 

AB = i 1 -i 2 +2i 3) 

BC =-2ii+2i 2 -3i 3 , 

find (i) the moment of BC about A, (ii) the moment of CB about 0, 
(iii) the moment of BC about the directed line OA, (iv) the moments 
of BC about the coordinate axes. 

30. If a and b are two vectors, prove that a times the moment of b 
about the line of action of a is equal to b times the moment of a about 
the line of action of b. 

31. If a (a) has a constant magnitude, show that 



32. If a = p cos a+q sin a, where p and q are constant vectors and 
a is a variable, show that 




= 0 v 


33. If a is a function of a variable a, show that 
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34. Given that the unit vectors i l3 i 2 and i 3 are independent of a 
variable u, evaluate f a du when 


(i) 

(ii) 

(iii) 


a 

a 

a 


= i 1 +2Mi 2 4-8M 3 i 3) 

= cos u +i 2 sin u - i 3 
2 ^ 2i 2 

4-u 2 ^ 4+a 2 


35. Evaluate the following integrals, in which p and q are constant 


vectors: 

/ (p+q u 2 )du, 

(i) 

(ii) 

/(pcosa +q sec 2 u)du, 

(iii) 

I tT *• 

J 4-i i 2 


36. Find the vector x(u) in each of the following cases, given that 
p, q and r are constant vectors: 


(i) 


dx 

du 


= pu 2 +q* 2 “, 


(ii) 

(iii) 


d 2 x 


= p cos u+ q sin u, 


d 2 x 

du 2 


(p sintt-q cos u) Xr. 


37. Find the general solutions of the differential equations 


(i) 


d 2 x 

did' 


^-6x = 0, 

du 


(ii) 



(iii) 



(iv) 

d i x _ d 3 x .. d 2 x 

__ b _ + n _ 

- dx 
6 Hu 


38. Find the general solutions of the following differential equations, 
given that p and q are constant vectors: 
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(i) 

3x = p(3u 2 +l), 

<f 2 x 

(ii) 

■^2 -4x = 16p cos 2 u, 

iii) 

d 2 x . dx „ „ 

M +2 Tu = - 6pe +5< * sm “ 


39. Find the general solutions of the following differential equa¬ 
tions, given that p and q are constant vectors : 


(i) 


d 2 x dx 

— 2 ~ + x = 2 p«“, 


du 2 


du 


(ii) 

//2v dx 

^-3| = 2K + 18q a 2 , 

(iii) 

„ ^ 2 x dx 0 c 
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Chapter II 


APPLICATIONS TO GEOMETRY 

15. Introduction. This chapter contains a treatment by vector 
methods of various elementary topics in geometry. This treatment is 
included in the present volume for two reasons. First, it indicates the 
ease and power which vector methods lend to studies in geometry. 
Secondly, it affords the student an opportunity to gain additional skill 
in the use of the vector operations introduced in the previous chapter. 

Proofs of some well-known theorems of plane geometry will first be 
given. Then a fairly broad treatment of solid analytic geometry will 
be presented, in which some of the more familiar formulas will be de¬ 
duced in vector form and by vector methods. Finally, the differential 
geometry of curves in space will be considered briefly. For more com¬ 
plete treatments of analytic geometry and differential geometry by 
vector methods, the reader is referred to the excellent books by F. D. 
Murnaghan, 1 W. C. Graustein 2 and C. E. Weatherburn. 3 

16. Some theorems of plane geometry. In this section we shall consider 
proofs by vector methods of two well-known theorems of plane ge¬ 
ometry. 

Theorem 1. The diagonals of a parallelogram bisect each other. 

Proof. Let us consider the parallelogram OABC in Figure 21. The 
diagonals cut each other at the point D. We must prove that D bisects 
both of the line segments OB and AC. 

For convenience we denote the vectors drawn from 0 to the points 

1 F. D. Murnaghan, Analytic geometry, Prentice-Hall, New York, 1946. 

2 W. C. Graustein, Differential geometry, The Macmillan Company, New York, 
1935. 

3 C. E. Weatherburn, Differential geometry of three dimensions, Cambridge 
University Press, Cambridge, England. Vol. 1, 1927. Vol. 2, 1930. 
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Figure 21 

A, B, C and D by a, b, c and d, respectively. Since D lies on the line 
OB, there exists a scalar a such that 

(16.1) d=ab. 

Also, from the figure we see that 

(16.2) d - OA+AD. 

Since D lies on the line AC there exists a scalar v such that 


AD - vAC. 

Hence we can write (16.2) in the form 

(16.3) d = a A-vAC- 

We now equate the above two expressions given for d in (16.1) and 

(16.3) , obtaining 

(16.4) a+p,4C=ab. 

The next step is to express all vectors in this equation as linear functions 
of any two vectors in the plane, say a and c. From the figure we see 
that 

AC = -a-)-c, b - a+c, 
whence (16.4) becomes 

a-f-r (—a+c) =a(a+c), 


35 



Since a and c do not have the same line of action it then follows that 
1-k-d=0, u - v — 0. 

We solve these equations for u and v, obtaining u = v = Thus Equa¬ 
tion (16.1) becomes d = \ b, from which it follows that D is the middle 
point of OB. 

We have now proved that the point of intersection D of the diago¬ 
nals is the middle point of one of these diagonals. From symmetry, D 
must also be the middle point of the other diagonal. 

Theorem 2. The medians of a triangle meet in a single point which 
trisects each of them. 

Proof. Let us consider the triangle OAB in Figure 22. The points 


B 



C, D and E are the middle points of the sides, and F is the point of 
intersection of the medians AD and BE. We must prove that F is a 
point of trisection of each of the three medians AD, BE and OC. 

For convenience we denote the vectors drawn from 0 to the points 
A, B, C, D,-E and F by a, b, c, d, e and f, respectively. Now F lies on 
the median AD. Thus there exists a scalar u such that 

(16.5) DF = uDA 

and we then have 

f = OD+uDA. 
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Similarly, since F lies on the median BE, we have 

f = OE+vEB, 

where v is some scalar. We now equate these two expressions for f, 
obtaining 

(16.6) OD+uDA = OE+vEB. 

The next step is to express all vectors in this equation as linear func¬ 
tions of any two vectors in the plane, say a and b. From Figure 22 we 
see that 

/jg y\ OB | b, OE = J a, 

DA = b-j-a, EB = a-f-b . 

Thus (16.6) becomes, after substitution from (16.7) and collection of 
like terms, 

{\-u-\v)a. = (|-| a _„)b. 

Since a and b do not have the same line of action, it follows that 
\-u-\v= 0, —»=0. 

We solve these equations for u and v, obtaining u = v = Thus 
Equation (16.5) becomes DF = J DA, from which it follows that F is 
a point of trisection of the median DA. 

We have now proved that the point of intersection F of two medians 
is a point of trisection of one of these medians. From symmetry, F must 
also be a point of trisection of the other of these medians. From sym¬ 
metry F must also be a point of trisection of the third median. 


Solid Analytic Geometry 

17. Notation. We shall now consider the analytic geometry of 
points in space. A right-handed set of rectangular cartesian coordi¬ 
nates is introduced, with origin at a point 0. Just as in Chapter I, we 
denote these coordinates by the symbols x 1} x 2 and x 3 , as shown in 
Figure 23. Unit vectors i 1; i 2 and i 3 pointing in the directions of the 
positive coordinate axes are also introduced, as shown. 
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Figure 23 


To denote a general point in space we shall use the letter X, and we 
shall call the vector OX the position-vector of X. For convenience we 
shall denote this vector also by the symbol x, and its components by 
the symbols (x 1} x 2 , x 3 ). We then have 

x = Xjii+^+^sis- 

The quantities x x , x 2 and x 3 are also the coordinates of the point X. 

We shall use the letters A,B,C, ■ ■ ■ to denote specific points in space, 
and shall denote the position-vectors of these points by a, b, c, • • •. The 
component of these vectors will be denoted in the usual way by the 
symbols (a 1 , a 2 , a 3 ), (b x , b 2 , b 3 ), (c x , c 2 , c 3 ), ■ ■ ■. We note that these 
quantities are also the coordinates of the points A, B,C, 

18. Division of a line segment in a given ratio. Let us suppose that we 
are given two points A and B, and that it is desired to find a third 
point C which divides the line segment AB in the given ratio m to n. 
Figure 24 illustrates the problem. If C lies between A and B, then 
0 < m/n<A~ oo; if C lies beyond B then -oo<m/n<-l; if C lies 
beyond A, then -1 < m/n < 0. In any event we have 
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A 


Figure 24 


(18.1) 


AC = CB 
m n 


If we now denote the position-vectors of A, B and C by a, b and c, 
respectively, then 

AC = c - a, CB = b - c, 

and (18.1) can then be written in the form 


n(c - a) = m{ b - c). 

Solving this equation for c, we obtain 

mb-t-na 

c =-'-- 

m-\-n 


(18.2) 


This formula expresses the position-vector c of the desired point C in 
terms of the known quantities a, b, m and n. 

In books on Analytic Geometry, formulas are usually given which 
express the coordinates of C in terms of m, n and the coordinates of 
A and B. It should be noted that (18.2) is entirely equivalent to these 
formulas, for these formulas can be deduced from (18.2) simply by 
equating the components of the left side of (18.2) to the components of 
the right side of (18.2). 


19. The distance between two points. Let us suppose that A and B are 
two given points, and that it is desired to find the distance d between 
A and B in terms of the position-vectors a and b of A and B. Figure 
25 illustrates the problem. Now 
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d = \AB\ . 

But AB = b - a . Thus 

d 2 = \AB\ 2 = AB-AB = (b - a) • (b - a) 

whence 


d = -\/(b - a) • (b - a). 



20. The area of a triangle. Let us suppose that A, B and C are three 
given points, and that it is desired to find the area A abc of the triangle 
ABC in terms of the position-vectors a, b and c of A, B and C. Figure 
26 illustraties the problem. 

c £ 



We first construct the parallelogram of which AB and AC form two 
adjacent edges, as shown. Then A abc is equal to one half the area of 
this parallelogram. But, by Theorem 1 of § 8, the area of this parallelo¬ 
gram is | AB X AC |. Hence we can write 

(20.1) A aic = * 9 , 
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where 

(20.2) c p=ABXAC. 

Now AB — b - a, AC = c - a. Thus 

<p = (b - a) x (c - a). 

This simplifies to 

(20.3) <p = bxc+cxa+axb. 

The required area of the triangle is thus given by (20.1), 9 being the 
magnitude of the vector given by (20.3). 

A property of the vector <p will now be recorded, for future use. 
Since this vector is equal to AB X AC we conclude that the vector cp is 
perpendicular to the plane of the triangle ABC , and its direction is that indicated 
by the thumb of the right hand when the fingers are placed to indicate the direction 
of the passage around the triangle from A to B to C. 

21. The equation of a plane. There are several ways in which a plane 
can be specified. For example, three points which are on the plane and 
do not lie on a single straight line can be given, or a line in the plane 
and a point on the plane but not on the line can be given. In each of 
several such cases we shall now deduce the equation which must be 
satisfied by the position-vector x of every point X on the plane. This 
equation will be referred to simply as the equation of the plane. In books 
on analytic geometry the equation of a plane usually appears as an 
equation which involves scalars only, and is satisfied only by the co¬ 
ordinates of points on the plane. We shall refer to this latter equation 
as the cartesian form of the equation of the plane. 

(i) To find the equation of the plane through a given point and perpendicular 
to a given vector. Let A be the given point and b be the given vector. 
Figure 27 illustrates the problem, the plane P being the plane in 
question. 

Let X be a general point on P, and let a and x denote the position- 
vectors of A and X, respectively. Now AX is perpendicular to b. Thus 

AX- b =0. 
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But AX = x - a, whence it follows that 
(21.1) (x - a) -b = 0. 

This is the desired equation of the plane P. 



The cartesian form of the equation of the plane P can be obtained 
readily from (21-1). It is only necessary to express (21.1) in terms of 
the components of the vectors involved. In this way we obtain the 
equation 

(*1 — fl l) (*2 — ^2"f" (*3 ~ %) — 0 . 

(ii) To find the equation of the plane through three given points. Let A, B 
and C be three given points. It is desired to find the equation of the 
plane P containing these three points. Figure 28 illustrates the problem. 

Let X be a general point on the plane P, and let a, b, c and x denote 
the position-vectors of the points A, B, C and X, respectively. In § 20 
we saw that the vector tp given by the relation 
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<p = bxc+cxa+axb 

is perpendicular to the plane P. Hence, by Problem (i) above the 
equation of P is 

(21.2) (x-a)-q>=0 . 

Since a-«p = a-(bxc+cxa+axb) = a-(bxc). 
Equation (21.2) can be written in the equivalent form 
(21-3) x-«p = a-(b x c). 



22. The vector-perpendicular from a point to a plane. The vector-perpen¬ 
dicular from a point D to a plane P is the vector with origin at D and 
terminus at the point on P nearest D. 

(i) To find the vector-perpendicular from a point D to a plane P through a 
given point and perpendicular to a given vector. Let A be the given point 
and let b be the given vector. Figure 29 illustrates the problem. We 
denote the position-vectors of the points A and D by a and d, respec¬ 
tively. If the point E is the foot of the perpendicular from the point D 
to the plane P, then DE is the vector-perpendicular from the point D to 
the plane P. We shall denote it by the symbol p. It is required to find 
p in terms of a and b. 
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Now p and b are parallel. Thus there exists a scalar constant K 
such that 

(22.1) p=ifb. 

From the figure we see that 

OD+DE+EA+AO = 0, 
or 

(22.2) d+rb + £d-a=0. 

Now b is perpendicular to EA. Hence b • EA = 0, and so scalar multi¬ 
plication of (22.2) by b yields 

(22.3) (d — a) • h-{-Kb 2 = 0. 

Thus 

r (a-d)-b 
A =-To-’ 


44 



and substitution for K in Equation (22.1) then yields 

( 22 . 4 ) p= (a - d) b b 

In particular, if the point D is at the origin, then d = 0 and 

(22.5) P = b 

(ii) To find the vector-perpendicular from a point D to a plane P through three 
given points. Let A, B and C be the three given points, with position- 
vectors a, b and c, respectively. Figure 30 illustrates the problem. It 
is desired to find the vector-perpendicular p in terms of a, b, c and d. 



In § 20 we saw that the vector cp given by the relation 

<p = bxc+cxa+axb 

is perpendicular to the plane P. Hence we may regard P as the plane 
through the given point A and perpendicular to the given vector <p. 
Thus, from Equation (22.4) it follows that the required vector-perpen¬ 
dicular is given by the relation 
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(a-d)-tp 

P= -^2 


Since a-«p = a - (b X c) it follows that 

( 22 . 6 ) 


a-(bXc) -dtp 

P =-li-<P- 


cp* 


23. The equation of a line. There are several ways in which a line in 
space can be specified. For example, two points on the line can be 
given, or two planes through the line can be given. In each of several 
such cases we shall now deduce the equation which must be satisfied 
by the position-vector x of every point X on the line. This equation will 
be referred to simply as the equation of the line. 

(i) To find the equation of the line through a given point and parallel to a 
given vector. Let A be the given point, with position-vector a, and let b 
be the given vector. Figure 31 illustrates the problem, L being the line 
in question. 




Figure 31 


Let X be a general point on L, and let x denote the position-vector 
of X. Now AX is parallel to b. Thus 


AX X b = 0. 
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But AX = x - a, whence it follows that 

(23.1) (x-a) xb = 0. 

This is the desired equation of the line L. 

(ii) To find the equation of the line through two given points. Let A and B be 
the given points, with position-vectors a and b, respectively. Figure 32 



illustrates the problem, L being the line in question. Now L is parallel 
to the vector AB, and AB = b - a. Thus, by Problem (i) above, the 
desired equation of L is 

(23.2) (x - a) x (b - a) = 0 . 

(iii) To find the equation of the line through a given point and perpendicular 
to two given vectors. Let A be the given point with position-vector a, and 
let b and c be the given vectors. Figure 33 illustrates the problem, L 
being the line in question. Now L is parallel to the vector b x c. Hence, 
by Problem (i) above, the desired equation of L is 

(23.3) (x- a) X (bxc) = 0. 

(iv) To find the equation of the line through a given point and perpendicular 
to the plane through three given points. Let A be the given point on the line, 
and let B, C and D be the given points on the plane. Figure 34 illustrates 
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0 


Figure 33 

the problem, L and P being the line and plane in question. We denote 
the position-vectors of A, B, C and D in the usual manner. Let us con¬ 
sider the vector tp given by the relation 


<p = cxd-fdxb+bxc. 



According to § 20, this vector is perpendicular to the plane P. Thus L 
is the line through A and parallel to <p, and hence by Problem (i) above 
the desired equation of Z, is 
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(23.4) (x-a)xcp = 0. 

24. The equation of a sphere. Let S be a sphere of radius a with 
center at a point C, as shown in Figure 35. If X is general point on 
the sphere S, then 

CX-CX = \CX\ 2 = a\ 

But CX = x - c. Thus 

(24.1) (x-c)-(x-c) = a 2 . 

This is an equation satisfied by the position-vector of a general point X 
on the sphere. It is thus the equation of the sphere. 



We shall now prove the following well-known property of a sphere: 
the angle at the surface of a sphere subtended by a diameter of the 
sphere is a right angle. For convenience, the origin of the coordinate 
system is chosen at the center of the sphere, as shown in Figure 36. Let 
D and E be points at the ends of a diameter, and let X be a general 
point on the sphere. We denote the position-vectors of these points in 
the usual manner. From the figure 

DX — x - d, EX = x - e. 
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Thus 


DX-EX = (x-d)-(x-e) = x-x - x-(d+e)+d-e. 

But d+e=0, and if a denotes the radius of the sphere then x- x — a 2 , 
d- e = -a 2 . Hence it follows that 

DX-EX = 0, 

and so DX is perpendicular to EX. 

x 



Figure 36 

25. The tangent plane to a sphere. Let S be a sphere of radius a with 
center at a point C, as shown in Figure 37. We shall now find 



Figure 37 


50 



the equation of the plane P which touches S at a given point X. 

Let Y be a general point on P. We denote the position-vectors of the 
various points in the usual manner. From the figure it follows that 

CY = CX+XY, 
or 

y-c = CX-\-XY. 

We now multiply this equation scalarly by CX NowCXCX = a 2 , and 
CX-XY = 0 because CX is perpendicular to XY. Thus we have 

(25.1) CX- (y-c) = a 2 . 

But CX = x- c. Thus Equation (25.1) becomes 

(x - c) • (y - c) = a 2 . 

This is the desired equation of the plane P. 

Differential Geometry 

26. Introduction. We shall consider here only a small portion of the 
differential geometry of curves in space. Rectangular cartesian co-or¬ 
dinates x lt a; 2 and x 3 are introduced, with origin at a point 0. The quan¬ 
tities x 1 , x 2 and x 3 denote the coordinates of a general point X with 
position-vector x. If q, i 2 and i 3 are unit vectors in the directions of 
the positive coordinate axes, then as before, 

(26.1) x = Xjii+Xjjia+Xsig. 

A curve consists of the set of points the position-vectors of which satisfy 
the relation 

x — x(«), 

wtiere x(u) is a function of a scalar parameter u. We shall consider only 
those parts of the curve which are free of singularities of all kinds. 

If the set of points comprising a curve all lie in a single plane, the 
curve is said to be a plane curve. If this set of points does not lie in a 
single plane, the curve is said to be a skew curve. 

It is convenient to choose as the scalar parameter the length s of the 
arc of the curve measured from some fixed point A. The quantity s is 
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positive for points on one side of A, and negative for points on the 
other side of A. The equation of the curve may then take the form 

X = x(s). 

The derivatives with respect to s of the function x(j) will be denoted by 
x', x", x'", etc. 

27. The principal triad. Let us consider a general point X on a curve 
C. The position-vector of X is x. We shall now define a set of three or¬ 
thogonal unit vectors j : , j 2 and j 3 at X. They are functions of the para¬ 
meter Sj and their derivatives respect to j- will be denoted in the usual 



Figure 38 


way by the symbols j/, j 2 ' and j 3 '. They are shown in Figure 38, and 
are defined by the conditions: 

(i) j x is tangent to the curve C, and points in the direction of 5 
increasing; 

(ii) j 2 lies in the plane of the vectors j x and j/, and makes an acute 
angle with j/; 

(iii) j 3 is such that the vectors jj,, j 2 and j 3 form a right-handed triad 1 . 

1 At points on the curve where j,' is equal to zero, these conditions are not suf¬ 
ficient for a unique determination of j 2 and j 3 . We exclude such points from con¬ 
sideration here. 
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The straight line through the point X and parallel to j 2 is called the 
principal normal to the curve . The straight line through X and parallel 
to j 3 is called the binormal to the curve. The vectors j x , j 2 and j 3 are 
called the unit tangent vector, unit normal vector, and unit binormal vector, 
respectively. The triad formed by these vectors is called the principal 
triad. The plane through X and perpendicular to j r is called the normal 
plane. The plane through X and perpendicular to j 3 is called the 
osculating plane. 

28. The Serret-Frenet formulas. Let Y be a point on the curve C and 
near the point X, as shown in Figure 38. We denote the length of the 
arc XY by the symbol As, and the vector XY by the symbol Ax. Let uf 
now consider the vector 

x' = limit — • 

As -*■ 0 As 

Now 

limit ^ X 1 = i. 

As -*• 0 A S 

Thus x' is a unit vector. Further, the vector Ax/A s lies along XY, and 
its direction then tends to that of j x as As tends to zero. Since j x is a unit 
vector, we can then write 

(28.1) = x'. 

The vectors j^, j 2 ' and j 3 ' can each be expressed as a linear func¬ 
tion of any three non-coplaner vectors. In particular, they can be ex¬ 
pressed as linear functions of the vectors jj, j 2 and j 3 , and we then 
have relations of the form 

jl ==fl lljld' a l2j2T a l3j3 J 

(28.2) j2'= a 2ljl + fl 22j2 + a 23 j 3! 

J3 =a 3lJl + a 32.kT a 33j3 j 

where the scalar coefficients are functions of the parameter s. Since 
the vectors jj, j 2 and j 3 are orthogonal unit vectors, they satisfy the 
relations 
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(28.3) 


ji-ji =1 > jrja=0, jrj 3 =0, 

j 2 \jl=0, j 2 ’j 2 =l, j 2 j3=0, 

j3jl=0» j 3 -j 2 = 0 > j3'j 3 = 1 • 

We differentiate with respect to s the first equation in the first line 
of (28.3). This yields the relation 

ji-ji'+ji'\h=0. 

Since in a scalar product the order in which the vectors appear is im¬ 
material, we can interchange the vectors in the second scalar product. 
It then follows that 

jrji' = 0. 

If we substitute here for j x ' from the first equation in (28.2), and then 
make use of Equations (28.3),we find that« n =0. Similarly a 22 =a 33 =0, 
and we may write 

(28.4) a n ==a 2 2 ==a 33 == 0. 

We now differentiate with respect to s the second relation in the first 
line of (28.3). This yields 

jrj 2 '+ji'-j 2 =0. 

If we substitute here for j/ and j 2 ' from the first two equations in 
(28.2), and then make use ofEquations (28.3), we find that a 12 +a 21 = 0. 
Similarly we can find two similar relations, and we have altogether 

(28.5) • a 12 -)-a 21 = 0, a 23 -\-a 32 = 0, a 3 i _ t~ a i 3 = 0 - 

So far, only conditions (i) and (iii) above have been used. By con¬ 
dition (ii) the vector j/ is to be in the plane of j* and j 2 . This can be 
true only if 

(28.6) fll3 =0. 

By condition (ii), the vector j x ' is to make an acute angle with j 2 . If 
this angle is denoted by a, then cos a must be positive. But 

Iji'l cos a=j 2 -j 1 '. 

If we substitute here for j x ' from Equations (28.2) and then use Equa¬ 
tions (28.3) we find that 
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Iji'l cos a =a 18 . 

Thus 

(28.7) a 12 > 0. 

We now define two quantities x and t by the relations 

(28.8) x=<z 12 , v=a 23 . 

Then, by (28.7) it follows that 

(28.9) x > 0, 

and because of Equations (28.4), (28.5), (28.6) and (28.8), we can now 
express Equations (28.2) in the form 

ji' = *ja > 

(28.10) j 2 '= Tj 3 -xj 1; 

j3=-^2 

These are the Serret-Frenet formulas. They were given originally by 
Serret (1851) and Frenet (1852) in an equivalent form which did not 
involve vectors. The quantities x and t, which are functions of the arc 
length s of the curve C, will be considered in some detail in the next 
section. 

29. Curvature and torsion. The quantity x appearing in Equations 

(28.10) is called the curvatures of the curve. It can be shown that 

.. A6 

x = limit — j 
An- o A s 

where A0 is the angle between the tangents to the curve C at the points 
X and Y in Figure 38. Thus x is the rate at which the tangent at the 
point X rotates as X moves along the curve. The reciprocal of x is 
called the radius of curvature, and will be denoted by the symbol p. 

The quantity t appearing in Equations (28.10) is called the torsion 
of the curve. It can be shown that 

t = limit —, 

Ar -*■ 0 As 

where AO is the angle between the binormals to the curve C at the 
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points X and Y in Figure 38. Thus t is the rate at which the unit bi¬ 
normal at the point X rotates as X moves along the curve. The re¬ 
ciprocal of t is called the radius of torsion , and will be denoted by a. 

To find x, we note from Equation (28.1) that 

(29.1) jx = x', jx' = x". 

Substitution from the second of these relations for j x ' in the first of the 
Serret-Frenet formulas then yields the equation 

(29.2) x"=xj 2 . 

We now multiply each side of this equation scalarly by itself, obtaining 


Because of (28.9), x is positive, and so 

(29.3) x = Vx"-x". 

To find the torsion r we differentiate Equation (29.2) with respect 
to s. This yields 

nr • r i /• 

X =XJ 2 +XJ 2 . 

We now substitute for j 2 ' from the second Serret-Frenet formula, ob¬ 
taining 

(29.4) x'" = x (Tj 3 -xj 1 )+x'j 2 . 

From Equations (29.1), (29.2) and (29.4) it now follows that 

x'-(x"xx'") = jx- [x j 2 X (xt ]3 - X^J+X 'j 2 )] 

= ji-[^'rj l +x 3 j 3 ] 

= x 2 t. 

Thus 


(29.5) 


T = -4 x'-(x"Xx'’ 


The curvature and torsion can be computed from Equations (29.3) 
and (29.5). We can express these equations in different forms. Now 
x = Aqij+Xjij+Xgia, 
x' =x 1 'i 1 +x 2 'i 2 +x 3 'i 3 , 
x" = *x"ix+* 2 "i 2 +* 3 "i 3 , 




56 



Substitution from these relations in Equations (29.3) and (29.4) then 
yields 


(29.6) 

X = Vx 1 " i + X 2 ' 

' S + * 3 " 2 

(29.7) 

-l"x 

II 

X- 

*1 X 2 

X " Y " 

X I X 2 

X Y 

*1 X 2 

V 

x a ” 

X 

X 3 


Since x can now be found, we can obtain the unit tangent vector j x 
and the unit normal vector j 2 by use of Equations (29.1) and (29.2). 
The unit binormal vector j 3 can then be found easily, since it is equal 
to jiXj 2 . We have the collected results 


(29.8) 


Jx = x 



j 3 = 


— x' X x" . 

X 


Let us now find the equation of the tangent to the curve at the point 
X. If Y is a general point on this tangent, the desired equation is easily 
seen to be 

(y - x ) x ji = o. 

Because of Equation (29.8), this can be written in the form 

(29.9) (y-x)xx' = 0. 

In the same way, the equations of the principal normal and binormal 
can be found in the forms 


(29.10) , (y - x) x x" = 0, 

(29.11) (y-x)x(x'xx") =0. 

The equation of the normal plane at the point X is easily seen to be 


(y- x )-ji = o. 


Because of Equation (29.8), this can be written in the form 

(29.12) (y — x) • x' = 0. 

In the same way we can find the equation of the osculating plane in the 
form 


(29.13) 


(y - x) • (x'Xi") = 0. 
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Problems 


1. Prove that the line joining the middle points of any two sides of a 
triangle is parallel to the third side, and is equal in length to one half 
the length of the third side. 

2. Prove that the lines joining the middle points of the sides of a 
quadrilateral form a parallelogram. 

3. If 0 is a point in space, ABC is a triangle, and D, E and F are 
the middle points of the sides, prove that 

OA + OB+OC = OD+OE+OF. 

4. If 0 is a point in space, ABCD is a parallelogram, and E is the 
point of intersection of the diagonals, prove that 

OA+OB+OC+OD - 4 OE. 

5. Prove that the line joining one vertex of a parallelogram to the 
middle point of an opposite side trisects a diagonal of the parallelogram. 

6. Prove that it is possible to construct a triangle with sides equal 
and parallel to the medians of a given triangle. 

7. Prove that the lines joining the middle points of opposite sides of 
a skew quadrilateral bisect each other. Prove also that the point 
where these lines cross is the middle point of the line joining the middle 
points of the diagonals of the quadrilateral. 

8. Prove that the three perpendiculars from the vertices of a triangle 
to the opposite sides meet in a point. 

9. Prove that the bisectors of the angles of a triangle meet in a 
point. Hint: the sum of unit vectors along two sides lies along the bi¬ 
sector of the contained angle. 

10. Prove that the perpendicular bisectors of the sides of a triangle 
meet in a point. 

11. If 0 is a point in space and ABC is a triangle with sides of lengths 
l, m and n, then 

l OA+m OB-Pn OC = (l+m+n) 65 , 

where D is the center of the inscribed circle. 

12. If ABC is a given triangle, the middle points of the sides BC, CA 
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and AB are denoted by D, E and F respectively, G is the point of inter¬ 
section of the perpendiculars from the vertices to the opposite sides, 
and H is the center of the circumscribed circle, prove that 

AG =2 HD, BG = 2 HE, CG = 2HF. 

Hence prove that 

GA+GB+GC= 2 GH. 

In problems 13-22 the points A, B, C and D have the following co¬ 
ordinates: A{-1, 2, 3), B(2, 5, -3), C{ 4, 1, -1), .0(1, 3, -3). 

13. Find the position-vectors of the points of trisection of the line 
segment AB. 

14. Find the distance between the points A and B. 

15. Find the area of the triangle ABC. 

16. Find the cartesian form of the equation of the plane through A 
and perpendicular to OB. 

17. Find the cartesian form of the equation of the plane through 
(i) the origin and the points A and B, (ii) the points A, B and C. 

18. Find the distance from the point D to the plane through A and 
perpendicular to OB. 

19. Find the distance from D to the plane ABC. 

20. Find the cartesian form of the equation of the line through A 
and parallel to BC. 

21. Find the cartesian form of the equation of the line through 
A and B. 

22. Find the cartesian form of the equation of the line through 
D and (i) perpendicular to the plane through A, B and C, (ii) perpen¬ 
dicular to BC and OC. 

23. A plane passes through a given point A with position-vector a, 
and is parallel to each of two given vectors b and c. Derive the equa¬ 
tion of this plane in the form 

(x - a) • (b X c) = 0 . 

24. A straight line L passes through a point A with position-vector 
a, and is parallel to a vector b. A vector p has its origin at a point C 
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with position-vector c, its line of action along the perpendicular from 
C to L, and its terminus on L. Show that 


P = 


- c - 


(a - c) -b 
b 2 


b. 


25. A straight line L passes through a point A with position-vector 
a, and is parallel to a vector b. A second straight line L' passes 
through a point A' with position-vector a', and is parallel to a vector 
b'. The vector p runs from L to L‘ along the common perpendicular. 
Show that 


where c = bxb'. 

26. Prove that if the torsion of a curve is equal to zero, the curve is 
a plane curve. 

27. Prove that 


=-3xx'j 1 +(x"-x 3 -xT 2 )j 2 + (2 x'T+XT')j s . 


28. If the position-vector x of a general point on a curve is given as 
a function of a parameter t, and if primes denote differentiations with 
respect to t prove that 


x =- 7 iA/x"-x"-r" 3 , t = 3 -iX'-(x" X x"') , 




; *' • 1 / „ s" v . x'xx" 

Ji= 7 > J 2 = ^(* - 7 x'), J3=^pr-- 


Also, derive the equations of the tangent, principal normal, binormal, 
normal plane and osculating plane in the following forms: 


tangent, 

principal normal, 

binormal, 
normal plane, 
osculating plane, 


(y-x) Xx' = 0^ 

(y " x) X (x" - — x') = 0; 

(y- x) X (x'xx") =0; 
(y-x)-x' =0; 

(y- x) - (x'xx") = 0. 


60 



29. The position-vector x of a general point on a circular helix is 
given by the relation 

x = a cos £ ij+a sin ti 2 -\-at cot ai 3 , 

where a and a are constants, and £ is a parameter. Find p, a and the 
principal triad. Answer: p = a cosec 2 a, cr = 2 a cosec 2a, j* = sin a 
(-ij sin £+i 2 cos £+i 3 cot a), j 2 = -ij cos t— i 2 sin £, j 3 = cos a (i 3 sin £ 
- i 2 cos £+i 3 tan a). 

30. The position-vector x of a general point on a curve is given by 
the relation 

x = a(3t - £ 3 )i 1 + 3ai 2 i 2 -(-fl(3£-|-£ 3 )i3, 

where a is a constant and £ is a parameter. Find p, a and the principal 
triad. Answer: p = a = 3a V2 jj = Y _1 (a ii+?i 2 +Yi 3 ) ; j 2 = Y'H-pii 
+ai 2 ), V2 j 3 = Y _1 (-“i] -Ma+Yis)* where a = 1 - £ 2 , p = 2£, 
Y = 1+^ 2 - 

31. The position-vector x of a general point on a curve is given by 
the relation 

x = a[(t - sin £)ii + (l - cos £)i 2 +£ i 3 ], 

where a is a constant and £ is a parameter. Find p and cr. Answer: 
p = a a 3/2 , fT l/2 , a = -a (J, where a = 3 - 2 cos £, fl = 2 - 2 cos £ 
+cos 2 £. 


61 



Chapter III 


APPLICATION OF VECTORS TO MECHANICS 
Motion of a particle 

30. Kinematics of a particle. The phrase ,,kinematics of a particle” 
refers to that portion of the study of the motion of a particle which is 
not concerned with the forces producing the motion, but is concerned 
rather with the mathematical concepts useful in describing the motion. 

Let us consider a moving particle. It is necessary to introduce a 
„frame of reference” relative to which the motion of the particle can 
be measured. For a frame of reference we take a rigid body. Such a 
body is one having the property that the distances between all pairs 
of particles in it do not vary with the time. We then introduce a set of 
rectangular cartesian coordinate axes fixed in the frame of reference. 
Figure 39 shows these axes and the associated unit vectors i x , i 2 and i 3 . 



*1 
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Figure 39 



Let the curve C in this figure be the path of the particle, and let the 
point X denote the position of the particle at time t. The vector OX 
is the position-vector of the particle. We denote this vector also by x. 
It is a function of the time t. 

The velocity v of the particle relative to the frame of reference, and 
the acceleration a of the particle relative to the frame of reference, 
are defined by the relations 


(30.1) 


dx _ dv _ d 2 x 

dt ’ a dt dt 2 ‘ 


The magnitude v of the velocity v is called the speed of the particle. 
Thus, velocity is a vector and speed is a scalar. We shall now compute 
various sets of components of the vectors v and a. 

(i) The components of the velocity and acceleration in the directions of rec¬ 
tangular cartesian coordinate axes. Let x lt x 2 , x 3 denote the rectangular 
cartesian coordinates of the point X in Figure 39. Then 

x = *ii i+* 2 i 2 -|-* 3 i 3 - 


If we now adopt the convention that a single superimposed 
notes a first time derivative, and a pair of superimposed dots 
a second time derivative, then 

v = ^ = *]il + *2»2+*3 i 3 , 


dot de¬ 
denotes 


a 


dv 

dt 


= i 3 i 3 . 


Thus the desired components of v and a are 


(30.2) *1,*2,* 3 ; *1, *2)*3- 

(ii) The components of the velocity and acceleration in the directions of the 
principal triad of the curve traced out by the particle. The curve C in Figure 
39 is the path of the particle. Let j x , j 2 and j 3 denote the principal 
triad at the general point X on C. The principal triad was discussed 
in § 27. If s denotes the arc length of C, then from Equations (28.1) 
and (28.10) we have 


(30.3) 
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where x is the curvature of C. Now 


dx dx . 
v = !t = Ts s ’ 

and because of the first equation in (30.3) we then have 

(30.4) v = s ji. 

Thus the velocity of the particle is directed along the tangent to its path, and the 
speed is v = s. 

Because of (30.4) we have 

dv .. . . djj. 

a= It =5Jl+f It' 

But because of the second equation in (30.3) we have 

(Si — s = xii 
dt ~ds KSj2 ’ 

and hence 

(30.5) a = 

Thus the acceleration a lies in the osculating plane of C. Also, the components 
of a in the directions of the tangent, normal and binormal are 


(30.6) 


s = v 


dv 

Is 5 


xi 2 = xv 2 = 


0 , 


where p is the radius of curvature of C. 

(iii) The components of the velocity and acceleration in the directions of the 
parametric lines of cylindrical coordinates. Let r, 0, x 3 be cylindrical coor¬ 
dinates of a general point X on the path C of the particle. Figure 40 
shows these coordinates. We introduce a triad of unit vectors k 1; k 2 , 
k 3 at O as shown; k x points toward the point X’ which is the projec¬ 
tion of X on the xyc^ plane, k 3 is equal to i 3 , and k 2 is such that the 
triad is right-handed. It will be noted that k ls k 2 , k 3 point in the 
directions of the parametric lines 1 of the cylindrical coordinates r, 6, 
x 3 at X. 


1 It will be recalled that the directions of the parametric lines of a coordinate 
system at a point X are those directions in which one of the coordinates increases 
while the other coordinates do not vary. 


64 





X 



If ij, i 2 and i 3 are the usual unit vectors associated with the rec¬ 
tangular cartesian coordinate axes in Figure 40, then 


and so 


(30.7) 


ki = ^ cos 0+i 2 sin 0, 
k 2 = -ij sin 0 +I 2 cos 9 j 

k 3 =i 3 , 


—j~ = 0 = (-ii sin 0+i 2 cos 0) 0 = k 2 0, 

0 = (_i x cos 0 -i 2 sin 0) 0 = ~M, 


dk 3 

dt 


= 0 . 


From Figure 40 it follows that the position-vector x of the particle is 
given by the relation x = rk 1 +x 3 k 3 . Thus 


v 


da 

dt 


db 

rk 1 +^3 k 3 + r -^- 1 +^3 


A,. 

dt 


Because of (30.7) we then obtain 
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(30.8) v = rk 1 J t-rQk. 2 J rX 3 k 3 . 

Thus the desired components of v are 

(30.9) r, rd, x s . 

From (30.8) we find that 

dv .. .• 

a = — = rkj-f (r0+r0)k 2 +x 3 k 3 
+r dt +r0 dt +x3 dt 

Substitution in this equation from (30.7) then yields 

(30.10) a = (r - r0 2 )k 1 + (2r0+r0)k 2 +x 3 k3 . 
Hence the desired components of a are 

(30.11) r-r&, 2r0 + r0 = ^ (r 2 0), x 3 . 


Of course, if the particle is confined to the x x x 2 plane then x 3 = 0, 
r = x, and we find from (30.9) and (30.11) that the components ofv 
and a in the directions of the parametric lines of the plane polar coor¬ 
dinates x, 0 are 

(30.12) *0; x-xW, 1 -^ (A) . 


31. Newton’s laws. The concept of force is intuitive. We can de¬ 
fine a unit force as that force which produces a standard deflection of 
a standard spring. Hence we can assign a numerical value to the mag¬ 
nitude of any force. 

We know that forces have magnitude and direction. It has been 
verified experimentally to within the limits of experimental error that 
forces obey the law of vector addition. Hence we shall assume that 
forces are vectors. The sum of two or more forces is sometimes called 
the resultant of the forces. 

The term “mass of a body” refers to the quantity of matter present 
in the body. We can define a unit mass as that mass which, when sus- 
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pended from a standard spring at a standard place in the earth’s 
gravitational field, produces a standard deflection of the spring. 
Hence we can assign a numerical value to the mass of any body. 

We now introduce the laws governing the motion of a particle. 
These laws, which were first stated by Isaac Newton and are called 
Newton’s laws, are as follows: 

(i) Every particle continues in a state of rest or uniform motion in a 
straight line unless compelled by some external force to change that 
state. 

(ii) The product of the mass and acceleration of a particle is pro¬ 
portional to the force applied to the particle, and the acceleration is in 
the same direction as the force. 

(iii) When two particles exert forces on each other, the forces have 
the same magnitudes and act in opposite directions along the line 
joining the two particles. 

In the second law, the acceleration of the particle enters. This accele¬ 
ration depends on the frame of reference employed. It thus appears 
that Newton’s second law cannot apply in all frames of reference. 
Those frames of reference in which this law does apply are called 
Newtonian frames of reference. A frame of reference fixed with respect to 
the stars is Newtonian, and in making an accurate study of any motion 
such a frame of reference should be used. However, for many problems 
we may consider the earth as a Newtonian frame of reference, when 
effects due to the motion of the earth are negligible. 

Let us now consider a particle of mass m acted upon by a force F. 
Let a denote the acceleration of the particle relative to a Newtonian 
frame of reference. Then according to Newton’s second law 

F —km a, 

where k is a constant of proportionality. It is customary to choose units 
of length, mass, time and force so that k is equal to unity. We then 
have 

(31.1) F=ma. 

There are three such systems of units in general use. These are indi¬ 
cated in Table 1, together with abbreviations commonly used for these 
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units. Thus, for example, when a force of one pdl. acts on a particle 
with a mass of one lb., the acceleration of the particle is one ft./sec. 2 . 
The systems of units in the second and third columns of Table 1 are 
called foot-pound-second sytems, or simply f.p.s. systems. The sy¬ 
stem of units in the fourth column is called the centimeter-gram- 
second system, or simply the c.g.s. system. 



f.p.s. 

c.g.s. 

Unit of length 

foot (ft.) 

foot (ft.) 

centimeter (cm.) 

Unit of mass 

pound (lb.) 

slug 

gram (gm.) 

Unit of time 

second (sec.) 

second (sec.) 

second (sec.) 

Unit of force 

poundal (pdl.) 

pound-weight (lb.wt.) 

dyne 


Table 1. Systems of units used in mechanics. 


The lb. wt. is the force exerted on a mass of one lb. by the earth’s 
gravitational field. If G denotes the acceleration due to gravity, ex¬ 
pressed in ft./sec. 2 , then 

1 lb. wt. - G pdl., 

1 slug = G lb. 

At points near the surface of the earth, G is approximately equal 
to 32. 

Equation (31.1), which governs the motion of a particle, may also 
be written in the equivalent forms 

( 31 -2) ' F = m d J r F = 

32. Motion of a particle acted upon by a force which is a given function of 
the time. When the force F acting on a particle is a given function of 
the time, Equations (31.2) can be solved by integration for the velocity 
v and position-vector x of the particle. 

As an example, let us suppose that 

F = 12 p+q cos t, 
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where p and q are given constant vectors. Because of the first equa¬ 
tion in (31.2) we then have 

dx 

m —j^ = 12p+q cos t , 

whence 

mv = J (12p+q cos t) dt. 

We now carry out this integration in the manner outlined in § 13, 
obtaining 


(32.1) mx = 12p£+q sin f+r, 

where r is an arbitrary constant vector. 

Now v = dx/dt. Thus from (32.1) we have 

(32.2) mx = y"(12p<+q sin t-\-r) dt 

= 6p t 2 - q cos <+r<-|-s. 


where s is an arbitrary constant vector. The arbitrary constant vec¬ 
tors r and s can be found if the initial values of x and v are known. 
If these initial values are x 0 and v 0 , it is readily found that 

r=mx 0 , s = mx 0 +q. 


33. Simple harmonic motion. Let 0 be a point fixed in a Newtonian 
frame of reference. Let us consider a particle moving under the action 
of a force directed toward 0, the force having a magnitude proportional 
to the distance from the particle to 0. If x denotes the position-vector 
of the particle relative to 0, then the force F acting on the particle 
satisfies the relation 

F = -kx, 

where k is a constant. From Equation (31.2) we then have 


or 


—kx = 


m 


d 2 x 

~d¥’ 


<Px k n 
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This is a differential equation of the type considered in § 14. According 
to the procedure demonstrated there, the general solution of this dif¬ 
ferential equation is 


x 


= c i 


>Y 


cos 1/ - t + c 2 sin 



where cq and c 2 are arbitrary constant vectors. These arbitrary con¬ 
stant vectors can be found if the initial values of x and v are known. 
If these initial values are x 0 and v 0 , it is readily found that 


c, = x. 


whence we have 


x = x 0 cos 


]/|*+v 0 j/fsin)/^. 


It will be noted that x is a linear function of x 0 and v 0 ; hence it follows 
that the motion of the particle is confined to the plane P containing 
the given vectors x 0 and v 0 . This result could have been anticipated, 
for the force F acting on the particle has no component perpendic¬ 
ular to the plane P. 


34. Central orbits. Let 0 be a point fixed in a Newtonian frame of 
reference. Let us consider a particle acted upon by a force F directed 
toward 0, the magnitude of F being a function of the distance from 0 
to the particle. The path of the particle is called a central orbit. It will 
be noted that the problem considered in § 33 dealt with a special type 
of central orbit. 

Denoting the vector OX by x, we then have 

F = F{x). 

The equation of motion is 

F = ma. 

Let Xo and v 0 be the initial values of x and the velocity v . The entire 
path of the particle will be in the plane P containing the vectors 
Xo and v 0 . Let x and 0 be polar coordinates in this plane. The com¬ 
ponents of F in the directions of the parametric lines of these coordi- 
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nates are -F, 0. Also, the components of a in these directions are given 
in Equation (30.12). Hence we have 

(34.1) -F=m{x~x 0 2 ), 

(34.2) °=f I, <*’»>• 

These are two equations from which * and 0 can be determined as 
functions of the time t. It is more convenient, however, to determine 
from (34.1) and (34.2) a single equation by elimination of the time 
variable. This single equation will now be deduced. 

We first introduce a variable y defined by the relation 

(34.3) y ~ \/x. 

Then from (34.2) we have 

fir 2 = const. = h 

whence 


(34.4) 

Then 


0 = Aj 2 . 




(34.5) i_- 4 g 6= _ 4 y!g. 

By substitution in (34.1) for x, 0 and x from (34.3), (34.4) and (34.5), 
we finally obtain 

< 34 - 6 > 

Now F is a function oiy alone. Once the form of this function has been 
assigned, we can find the path of the particle by solving Equation 

(34.6) . 

Let us now consider the special case when F varies inversely as the 
square of x. Then we can write 

F = y my 2 , 

where y is a constant, and Equation (34.6) becomes 
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dW ' y h?' 

The general solution of this equation, expressed in terms of x, is 

1 Y 

(34.7) - = p + c x cos 0+4 sin 0 , 

where c l and c 2 are arbitrary constants. These constants can be found 
if the initial values Xq and v 0 are known. It can be shown that (34.7) 
represents either an ellipse, parabola or hyperbola, depending on the 
values of x 0 and r 0 . 


Motion of a system ofparticles 

35. The center of mass of a system ofparticles. Let us consider a system 
of N particles. We denote their masses by the symbols m 1 , m 2 , m 3y 
■ ■ ■, m^r. The total mass m of the system is then given by the relation 

w 

(35.1) m = E mj. 

i- i 

We denote the coordinates of the particle of mass m f by the symbols 
(Xjd x j2 , x j3 ). The position-vector x 7 of this particle then satisfies the 
relation 

*/ = xji i 1 -\-Xj 2 i 2 +x js i 3 (j = 1,2, • • •, TV). 

We have then a set of 3jVscalars x jk (j = 1, 1, 3, • • •, W; k = 1, 2, 3) 
which denotes the coordinates of the particles. 

The center of mass of the system of particles is defined to be the 
point C with position-vector x c determined by the equation 

w 

(35.2) m x c = 2 mj %j. 

i -1 

The center of mass is sometimes called the mass center or centroid 
or center of gravity. 

If the distances between the individual particles in a system remain 
unaltered, as mentioned previously the system is called a rigid body. 
A rigid body often consists of a continuous distribution of matter, and 
in this case the summations in Equations (35.1) and'(35.2) above must 
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then be replaced by integrations. Thus, if p is the density of matter 
in the body, V is the region occupied by the body, dV is the volume 
of an element of the body and x is the position-vector of a point in 
dV, then 

(35.3) m = y p dV, 

v 

(35.4) mx. c — f P x dV. 

v 

36. The moments and products of inertia of a system of particles. Let us 
first consider a single particle of mass m. Let l denote the length of the 
perpendicular from the particle to a line L, as shown in Figure 41. 



The moment of inertia of the particle about the line L is defined to be 
the scalar I given by the relation 

I = ml 2 . 

Let p and q denote the lengths of the perpendiculars from the particle 
to two perpendicular planes P and Q_, as shown in Figure 42. The 
product of inertia of the particle with respect to these two planes is 
defined to be the scalar K given by the relation 

K = mp q. 

The moment of inertia about a line of a system of particles is defined 
to be the sum of the moments of inertia about the line of the individual 
particles. Also, the product of inertia with respect to two planes of 
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Figure 42 


a system of particles is defined to be the sum of the products of inertia 
with respect to the two planes of the individual particles. 

We now consider a set of N particles and introduce a set of rec¬ 
tangular cartesian coordinate axes with origin at a point 0. As in 
§ 35 we denote the masses of these particles by rrij (j = 1,2,3, • • -, W) 
and their coordinates by the 3JV symbols x jk (j = 1, 2, •••, J\f; 
k = 1, 2, 3). The moments of inertia of this system of particles about 
the three coordinate axes are denoted by A> I 2 , and I 3 . It is easily 
seen that 

N 

A = S nij (x 2 j2 +x z j 3 ) , 

j = 1 

(36.1) A = £ rrij {x\+x\) , 

j = 1 

A = £ mj (x\+x\) . 
i= i 

The products of inertia of this system of particles with respect to the 
three coordinate planes, taken in pairs, are denoted by K x , K 2 and K z . 
It is easily seen that 

x 

K\ = £ mj xj 2 x ]3 , 
j - 1 
x 

(36.2) K 2 = £ mj x, 3 % , 

j -1 

N 

kC ;j = £ Ulj Xj-y Xj 2 . 
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Of course, if the system of particles forms a rigid body, the summations 
in Equations (36.1) and (36.2) must be replaced by integrations. 

If 7 1; / 2 , I 3 , K x , IC 2 and K z are known, the moment of inertia I of 
the system about any line L through the origin can be found easily. 



In order to prove this we let pj denote the length of the perpendicular 
from the particle of mass m.j to the line L, as shown in Figure 43. Then 

jv 

I = 2 mjpj 2 . 
i = 1 

But from the figure we see that 

pj = Xj sin 0 = |b X Xj\ , 

where x ; - is the position-vector of the particle of mass m.j, Xj is the 
magnitude of x ; -, b is a unit vector on L, and 0 is the angle between 
Xj and b. The components of the vector b x x ; - are 

K X j3 ~ ^3 X j2> ^3 X jl ~ j K X }2 ~ ^2 X jl , 

and hence />/ is equal to the sum of the squares of these three com¬ 
ponents. Thus 

JV 

I = h mj [(^2*J3 — ^3 X j2 ) 2 (^3 X j1 ~ ^1 X ii) 2J S~ (^1 X j2 ~ ^2 *;l) 2 ] 

J “ 1 

= V 2 mj (V+V)+V 2 m i (V+% 2 ) + 

; = i i = i 
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N N 

+ b z 2 S mj {xji~\-Xj 2 2 ) - 2 b 2 b z S rrij x j2 x jZ 

j= i j = i 

JV JV 

-2 b 3 b 1 H mj x j3 x jx -2 b 1 b 2 1 l mj x jl x j2 . 

i - 1 j = l 

Because of (36.1) and (36.2) we then have 

(36.3) / = /, b*+I 2 b 2 *+I 3 bf - 2 Kj, b 2 b 3 - 2 K 2 b 3 b x - 2 r, i 1 b 2 . 

It will be noted that ^ b 2 and b 3 , which are the components of the 
unit vector b on the line L, are also the direction cosines of L. Equa¬ 
tion (36.3) is the desired equation which permits a simple determi¬ 
nation of the moment of inertia / of a system about any line L through 
the origin, once / 1; I 2 , I 3 j K 1} K 2 and K 3 have been found. 

If it happens that K 1 = K 2 = K 3 = 0, the coordinate axes are said 
to be principal axes of inertia at the point 0. It can be proved that at 
every point there is at least one set of principal axes of inertia. 1 In 
many cases, principal axes of inertia can be deduced readily by con¬ 
siderations of symmetry of the system of particles. For example, at 
the center of a rectangular parallelepiped the principal axes of inertia 
are parallel to the edges of the body. 

We shall now state without proof two theorems the proofs of which 
are very simple and may be found in almost any text book on calculus. 

The theorem of perpendicular axes. If a system of particles lies entirely 
in a plane P, the moment of inertia of the system with respect to a line 
L perpendicular to the plane P is equal to the sum of the moments of 
inertia of the system with respect to any two perpendicular lines 
intersecting L and lying in P. 

The theorem of parallel axes. The moment of inertia I of a system of 
particles about a line L satisfies the relation 

I = I'+rn l 2 , 

where /' is the moment of inertia of the system about a line L' parallel 
to L and through the center of mass of the system, m is the total mass 

1 See J. L. Synge and B. A. Griffith, Principles of Mechanics, McGraw-Hill Book 
Co., 1942, pp. 311-321. 
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of the system, and l is the perpendicular distance between L and L'. 

In most books of mathematical tables there are listed the moments 
of inertia of many bodies with respect to certain axes associated with 
these bodies. By the use of such tables together with Equation (36.3) 
and the above two theorems, it is frequently possible to determine 
rapidly the moment of inertia of a body with respect to any given line. 

37. Kinematics of a rigid body. Let us consider a rigid body which 
is rotating about a line L at the rate of to radians per unit time. The 
body is said to have an angular velocity. We can represent this angular 
velocity completely by an arrow defined as follows: its length repre¬ 
sents the scalar to some convenient scale; its origin is an arbitrary 
point on L ; its line of action coincides with L ; it points in the direction 
indicated by the thumb of the right hand when the fingers are placed 
to indicate the sense of the rotation about L. To prove that angular 
velocity is a vector, it is then only necessary to prove that it obeys the 
law of vector addition. This will be done presently. We shall antic¬ 
ipate this result, and denote angular velocities by symbols in bold¬ 
faced type. 

We shall first determine the velocity v of a general point X in'a body 
which has an angular velocity to. Let 0 denote a point on the line of 
action of to, and let x denote the vector OX, as shown in Figure 44. 



Let 0 denote the angle between to and x, and let p denote the length 
of the perpendicular from X to the line of action of to. The displace¬ 
ment dx of the point X in time dt has the following properties: 

(i) its direction is perpendicular to both to and x; 

(ii) its direction is that indicated by the thumb of the right hand 
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when the fingers are placed to indicate the sense of the rotation 
0 from to to x; 

(iii) its magnitude is pcoi it, which is equal to xu>dt sin 9, x being the 
magnitude of the vector x. 

In view of the definition in § 8 of the vector product of two vectors, 
it then appears that dx = t oXx dt. Thus 


If the point 0 is fixed in a frame of reference S, the velocity v of the 
point X relative to S is then 

(37.1) v = to x x. 

On the other hand, if the point 0 has a velocity u relative to a frame 
of reference S, the velocity of the point X relative to S is 

(37.2) v = u+toXx. 

We shall now prove that angular velocity obeys the law of vector 
addition. Let us consider a body which is rotating simultaneously 
about two lines L and V which intersect at a point 0 fixed in a frame 
of reference S. These angular velocities can be represented by the 
arrows to and to' in Figure 45. Let X be a general point in the body, 



Figure 45 

with position-vector x relative to 0. The two angular velocities impart 
to X the two velocities to x x and to' x x which, being vectors, can be 
added to yield the resultant velocity 

(37.3) v = toxx+w'xx. 

To complete the proof we must show that (37.3) can be written in the 
form v = to"Xx, where to" is an arrow obtained by the application 
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of the law of vector addition to the arrows to and to'. Even though 
angular velocity has not been assumed to satisfy the law of vector 
addition, Equation (8.5) may be applied to the two products in (37.3) 
to yield 

V X = % 3 ~ t0 3 2 *3 ~ W 3 

= (w 2 d-0)' 2 ) X 3 - (c0 3 -j-(d' 3 ) x 2 , 

and two similar expressions for v 2 and v 3 . Hence we can write 

v = to" X x, 

where to " is an arrow having components tiq+ciq', <a 2 +w' 2 , oj 3 +o)' 3 . 
But these are the components of the vector obtained by application 
of the law of vector addition to the arrows to and to'. Hence to" is 
equal to the vector sum of to and to', and so angular velocity is a 
vector. 

It will be noted that two angular velocities can be added only when 
their lines of action have a point of intersection, and that the line of 
action of the sum passes through this point of intersection. 


38. The time derivative of a vector. Let us consider a set of rectangu¬ 
lar cartesian coordinate axes with origin 0 fixed in a frame of reference 
S, and with axes rotating relative to S with angular velocity to. Then 
the line of action of to passes though 0. 

If i-L, i 2 and i 3 are the usual unit vectors associated with these coor¬ 
dinate axes, then the velocities relative to S of the terminuses of these 
vectors are 

d i x di 2 di 3 

dt 5 dt ’ dt 


But by the previous section these velocities are 

to X ii, to X i 2 , to X i 3 . 

Hence 
(38.1) 


di± 

dt 


d i 9 


-WXl1 ’ dt ~~ w Xl2 ’ dt 


dU 


= to Xl 3 . 


Let a be a vector with components a x , a 2 , a 3 relative to the rotating 
coordinate axes. Then 
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a — di i 1 +a 2 * 2 - l~ a 3 *3 j 

and the time derivative of a, relative to 6', is then 

da da,. , da 2 . , da 3 . ddi 2 d i 3 


1 dt 


Because of Equations (38.1) we can write the last three terms in the 
form 

aqto X ij+a 2 w X i 2 + a 3 t» X i 3 , 
which reduces to w X a. Hence we can write 


(38.2) 
where 

(38.3) 


da 

dt 


8a . 

¥ +mXa : 


Sa 

8 t 


dci-t « da$* , da • 
**■+ 


dt 


! *2 + ' 




Equation (38.2) expresses dajdt as the sum of two parts. The part 
Sa/St is the time derivative of a relative to the moving coordinate 
system. The part wXais the time derivative of a relative to S when 
a is fixed relative to the moving coordinate system. 

When the origin of the coordinate system is not at rest relative to S 
but has a velocity u, Equations (38.1) still hold, and hence also does 
Equation (38.2). 


39. Linear and angular momentum. Let us consider a particle of mass 
m, with a position-vector x relative to a point 0 fixed in a frame of 


m 



reference .S'. Let v denote the velocity of the particle, as shown in 
Figure 46. The linear momentum of the particle is a vector M defined 
by the relation 
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M = m v. 


The angular momentum of the particle about the point 0 is by defi¬ 
nition the moment of M about 0. We shall denote it by the symbol h. 
Hence, by § 10 where the moment of a vector about a point is consid¬ 
ered, we have 

h - xxM - xXfflv — mxXv. 

Let us now consider a system of JV particles. As' before, we denote 
the mass and position-vector of th ej-th particle by and x,-, respect¬ 
ively. Also, we denote the velocity of this particle relative to S by v ; -. 
Then for this system the linear momentum M and the angular mo¬ 
mentum h about 0 are defined by the relations 
jv jv 

(39.1) M = S m;Vj, h = 2 mjXj Xv;. 

j = l ; = l 


Theorem. The linear momentum of a system of particles is equal 
to the product of the total mass of the system and the velocity of the 
center of mass of the system. 

Proof. The position-vector x c of the center of mass of the system 
is given by Equation (35.2). We differentiate this equation with 
respect to the time t, obtaining 


But 


dx c £ dxj 
m -j— —l*mj — . 
dt j= i dt 


dx c 

dt 


v Cj 


dxj 

dt 


Vj, 


where v c is the velocity of the center of mass C relative to S. Hence 

JV 

(39.2) mv c = S m-Vj = M. 

= i 

This completes the proof. 


Let us now suppose that the'system of particles constitutes a rigid 
body, and that the body is rotating about the point 0 which is fixed 
in the frame of reference S. The body then has an angular velocity to 
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with a line of action which passes through 0. The velocity relative to 
S of the j-tk particle in the body is then 

Vy = U) X Xj, 

and by Equation (39.1) the angular momentum of the system about 
0 is then 

x 

(39.3) h = S mjXjX(iaXxj). 

j = 1 

Because of the identity (9.3), we can then write 

x 

h = E rrij [co xf - x y (xy • co)] . 
j-i 

Now let us introduce coordinate axes with origin at the point 0 fixed 
in S. The directions of these coordinate axes need not be fixed in S. 
As before we denote the coordinates of th ej-th particle by (x yi , x j2 , 
x j3 ) . The component h x of h then has the value 
x 

hi = S mj [coi (xjS+xrf+Xjr?) ~ Xj-y (x J1 u 1 +Xy 2 u 2 +Xy 3 to 3 )) 
i — l 

X XX 

= <0i S ntj (x y2 2 +x y3 2 ) - co 2 S rrij x ;1 x j2 - co 3 S m.j x j3 x ;1 
j— i j = i j = i 

= — K s (x> 2 — K 2 C0 3 , 

where I x , K 2 and K 3 are moments and products of inertia defined in 
§36. There are similar expressions for h 2 and h 3 . We have finally 

— -fi^i — K 3 (x> 2 — K 2 o 3 , 

(39.4*) h 2 = -IC 3 a^-\-I 2 oi 2 — K ^ 3 , 

h 3 = —K 2 (a i — A' 1 co 2 + I 3 u> 3 • 

Let us now consider a rigid body which is moving in a general 
fashion relative to a frame of reference S. Let us introduce coordinate 
axes with origin at the center of mass C of the body. The directions 
of the coordinate axes need not be fixed in the body, however. We 
may consider the body as having a velocity of translation v c plus an 
angular velocity about a line through C. Then, as seen in § 37, the 
velocity of th ej-th particle can be expressed in the form 

v j = Vc+wXXy. 
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Hence the angular momentum h of the system about the center of 
mass C has the value 

JV 

h = ^i mj Xj X ( v c+*° x x j) 

/ * \ * 

= y . Xy J X V c +Sm, Xy X (w X Xy) . 

By Equation (35.2) the first sum is equal to mx c . Since the origin of 
the coordinate system and the center of mass C of the body coincide, 
x c = 0. Thus 

jv 

h = E ■kj Xj X (to X Xj) . 
i -1 

The right side of this equation is the same as the right side of Equation 
(39.3). Hence in the present case the components ofh are also given 
by Equations (39.4). 

We have then the important result: Equations (39.4) may be used for 
the determination of the components of the angular momentum h of a rigid 
body about either a fixed point 0 in the body or the center of massC of the body. 
In the two cases the origin of the coordinates is at 0 and C, respectively , the 
directions of the coordinate axes being quite general. Equations (39.4) 
cannot be used in the case of the angular momentum of a rigid body 
about a moving point which is not the center of mass of the body. 

40. The motion of a system of particles. Let us consider a general 
system of N particles. Let mj denote the mass of th ej-th particle, and 
let Vy denote its velocity relative to a Newtonian system. The forces 
acting on the j-th particle can be divided into two groups called 
internal forces and external forces. Internal forces are those due to other 
particles in the system. External forces include all other forces. Let 
F jt denote the internal force exerted on the j-th particle by the k-th 
particle, and let Fy denote the total external force exerted on the 
j-th particle. 

Theorem 1. The rate of change of the linear momentum of the 
system is equal to the sum the external forces acting on the system. 
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Proof. Applying to the j-th particle Newton’s Second Law as 
stated in § 31, we have 

(40.1) Wj-F,+fF*. 

We now sum the jV equations in (40.1), obtaining 

JV J X XX 

(40.2) £ = £ F,-+ £ £ F*. 

i = i at j = i j=i t = i 


Because of Newton’s Third Law, as stated in § 31, Fy* = -F) k . Thus 
the double sum in Equation (40.2) vanishes, and we can then write 

(40.2) in the form 


(40.3) 



where M is the linear momentum of the system and F is the sum of 
the external forces acting on the system. 

Theorem. 2. The center of mass of a system of particles moves like 
a particle with a mass equal to the total mass of the system acted upon 
by a force equal to the sum of the external forces acting on the system. 

Proof. In § 39 we saw that M = mv c , where m is the total mass of 
the system, and v c is the velocity of the center of mass of the system. 
Thus Equation (40.3) can be written in the form 

( 40 . 4 ) m ^ = F . 

This completes the proof. 

41. The motion of a rigid body with a fixed point. Let us now con¬ 
sider a system of particles which constitutes a rigid body with a point 
0 fixed relative to a Newtonian frame of reference. 

Theorem 1. The rate of change of the angular momentum of the 
body about 0 is equal to the total moment about 0 of the external 
forces. 

Proof. Let us introduce coordinates with origin at 0. Then 
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where v ; -, x,- and t have the usual meanings. By Equation (39.1), the 
angular momentum h of the body about the fixed point 0 is 


JV 

h = E m.j Xj x Vy, 

J = 1 


and so 


(41.2) 

dh AiR 

It - A + B > 

where 

X dx- 

A = E mj-jfxxj, B = : 

j' = i j- 

Because of Equation (41.1) we have 


JV 

A = E m ; - v ; - x Vy 

j - 1 

Equation (40.1) gives an expression for 
have 

B = E XjX (F y + E F jt ) 

j = i * = i 

where 

(41.3) 

JV 

G = 2 XyXFy 
; = l 


JV 


d,Vj 

dt 


JV JV 

H=E SxyXFy*. 

j = l Jt=l 


It will be recalled that F,- is the external force acting on the j-th par¬ 
ticle and F ) k is the internal force exerted on the j-th particle by k-th 
particle. We note that G is the sum of the moments about 0 of the 
external forces. Now 

JV JV JV JV 

H = S 2x,xF ft = 2 E x,xF tj , 

j = 1 t = 1 * = 1 3 - 1 

Thus 

JV JV 

(41.4) 2 H = E S (xyXF^+^xF*,-). 

j —\ k = 1 


85 



But F kj = -F jk . Thus (41.4) becomes 

2H = 2 2 (xj - x k ) X Fj* . 

7 = 1 A: = 1 


Since the lines of action of the vectors x ; - x k and F* coincide, their 
vector product vanishes. Hence H = 0, and B = G, so Equation 
(41.2) reduces to the form 


(41.5) 


dh 

1* =G ’ 


where h is the angular momentum of the system about the fixed 
point 0, and G is the total moment about 0 of the external forces. 
This completes the proof. 

We have placed the origin of the coordinate system at the fixed 
point 0. Let us now choose as coordinate axes a set of principal axes 
of inertia of the body at 0. (Principal axes of inertia are defined in 
§36.) Then the products of inertia K±, K 2 , K 3 all vanish, and from 
Equations (39.4) we obtain for the components of the angular momen¬ 
tum h of the body about 0 the expressions 

(41.6) h-L = COj , h 2 = /2 W 2 3 ^3 ~ ^3^3 3 


where I t , / 2 , 1 3 are the moments of inertia of the body about the coor¬ 
dinate axes, and oq, co 2 , co 3 are the components of the angular velo¬ 
city u) of the body about 0. 

In most cases the coordinate axes will be fixed in the body and will 
hence have an angular velocity to about 0. However, in a few special 
cases when the body has a certain symmetry it will be found possible 
and desirable to choose coordinate axes not fixed in the body. To 
include such special cases we denote the angular velocity of the axes 
about 0 by £2, which may or may not differ from u>. According to 
Equation (38.3) we then have 


or 

(41.7) 


dh Sh _ , 

It ^ lit +S2xh 


dh 


A 2 i 2 T ^3X3 •-■( QA3 q 

+ (G3A1 _ ^1^3) i 2 + (^A — ^2^1) *3 ■ 
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f rom this equation we can read off the components of the vector 
dh/dt. According to Equation (41.5) these components are equal to 
the components of G. Hence we have the equations 

H 3 A 2 = G^, 

— HjAg = G 2 , 
h 2 -\-Q-Ji 2 ^ 2^1 ~ ^3j 

where G 1 , G 2 , G 3 are the components of G. Because of Equations 
(41.6) these relations can be written in the form 

/lO)l —/ 2 w 2^3 - t - -^3 W 3^2 = ^ 1 > 

(41.8) -f 2 <b 2 — f3 (J 3^ld _ A W 1^3 = ^2> 

I 3&3 ~ ■fl tji 1^2 - t“-^2 W 2^1 = ^*3 • 

In the case when the coordinate axes are fixed in the rigid body, 
then £2 = w and so (12.8) reduce to the form 

/jWj — (/ 2 — / 3 ) co 2 u 3 = Gj, 

(41.9) / 2 w 2 — (/ 3 — Ii) w 3 ^i — G 2 , 

^ 3^3 — (A ~ ^ 2 ) w l w 2 = ^3" 

These equations are called Euler’s equations of motion. 

Theorem 3. The total moment about 0 of the gravity forces acting 
on a system of particles is equal to the moment about 0 of a single 
force equal to the resultant of the gravity forces and acting at the 
center of mass of the system. 

Proof. Let us introduce a coordinate system with origin at the 
point 0. Let k be a unit vector in the direction of the gravity forces. 
Then the gravity force acting on th ej-th particle is mygk, and the total 
moment about 0 of the gravity forces is 

N IN 

G' = 2 x,- X m,gk = I 2 
j — 1 \j = 1 

But by Equation (35.2) we have 

N 

2 Mj xj =mx c 



87 



where m is the total mass and x c is the position-vector of the center of 
mass. Thus 

G" = 7nx c Xgk = x c X (mgk). 

But x c x(mgk) is the moment about 0 of a single force mgk equal 
to the resultant of the gravity forces and acting at the center of mass 
C of the system. This completes the proof. 

Example 1. A sphere of radius a is placed on a rough plane which makes 
an angle a with the horizontal, and is then released. Find the distance the 
sphere moves down the plane in time t. 

Figure 47 shows the configuration of the system at a general time t. 



The center of mass of the sphere is at its geometrical center. The 
point Q_ is the initial point of contact of the sphere and plane, and in 
the time t the sphere has rolled through an angle <£>, as shown. The 
point on the sphere which is in contact with the plane is at rest. 
Hence the sphere has a fixed point, and we select this point as the origin 
0 of the coordinate system. We must select the coordinate axes to 
coincide with principal directions of inertia of the sphere at 0. This 
requirement is satisfied if the unit vectors ij and i 3 are chosen as 
shown in the figure. The unit vector i 2 is then perpendicular to i x and 
i 3 , and points up from the page. We note that the coordinate axes 
are not fixed in the body, and that consequently Equations (41.8) 
apply. 

The external forces acting on the sphere consist of gravity and the 
reaction of the plane. Because of Theorem 3 above, the forces exerted 
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by gravity on all the particles of the sphere may be replaced by a single 
force mgj acting at C, as shown, where j is a unit vector. The reaction 
of the plane is a force which may be resolved into a force M i 3 normal 
to the plane and a force -T i x along the plane, as shown. The moment 
G of the external forces about 0 is given by the relation 

G = OCXmgj, 

since the moments of N and T about 0 are equal to zero. But 
OC = ai 3 , j = ij sin a - i 3 cos a. 

Thus 

G = ai 3 X mg (i x sin a - i 3 cos a) 

= mgai 2 sin a, 


(41.10) G x = 0, G 2 =mga sin a., (? 3 =0. 

The coordinate axes have no angular velocity, so 

(41.11) £1 1 =Q 2 =Q 3 =0. 

The angular velocity co of the sphere about 0 is 

co = 6 i 2 , 

where the superimposed dot denotes differentiation with respect to t. 
Thus 


(41.12) Gq=0, w 2 =6, co 3 = 0 . 

The moments of inertia of the sphere about the coordinate axes are 

(41.13) I x = / 2 = £ma 2 , I 3 = \ma 2 . 

We now substitute in Euler’s equations (41.8) from Equations 
(41.10), (41.11), (41.12) and (41.13) to obtain the relation 

M 2 = mgai 2 sin oc. 


Thus 


* _ 5 S sin a 
la 


and two integrations then yield 
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■since $=6=0 when t = 0. 

If £ is the distance the sphere has rolled down the plane in time t, 
then z = and we have 

(41.14) £ = T \ gt 2 sin a. 

Example 2. Two shafts are attached to the comers A and C of a rectangular 
plate ABCD, in such a way that the axes of the shafts are continuations of 
the diagonal AC. The shafts turn in two bearings each at a distance c from the 
center of the plate. The system is made to rotate at a constant rate of (a radians 
per unit time. Find the forces exerted on the bearings. 

The system is shown in Figure 48. We choose the center 0 of the 



plate as the origin of the coordinate system. We also choose the unit 
vectors i x and i 2 parallel to edges of the plate, as shown. The unit 
vector i 3 is then perpendicular to the plane of the plate. The directions 
of these three vectors are principal directions of inertia of the plate 
at 0. The moments of inertia of the plate about the coordinate axes 
are 

(41.15) / 1 =ima 2 2 , / 2 = \ma*, f = \m(a^+af ), 

where m is the mass of the plate, and 2 a x and 2a 2 are the lengths of 
the edges. 



The plate has an angular velocity to the line of action of which is 
the diagonal AC and the magnitude of which is the given constant 


w. If a is the angle between to and i x , then 

(41.16) 

a 2 

tan a. — — 
ai 

and 

to = odq cos a+coia sin a. 

Thus 


(41.17) 

<«h = co cos a, cd 2 = to sin a, co. 


Since the coordinate axes are fixed in the body, Euler’s Equations 
(41.9) apply. We substitute in these equations from (41.15) and 
(41.17), obtaining the relations 

= 0, G 2 = 0, Gj=|»i (a j 2 - a 2 2 ) to 2 sin a cos a. 

Thus the moment G about 0 of the external forces is normal to the 
plate and rotates with it. Hence the forces exerted on the shaft by the 
bearings must be in the plane of the plate. Let us denote these forces 
by R and -R, as shown in Figure 48. We must then have 

2 cR = G 3 

whence we find that 


(41.18) 


R = -§ 7 -.m (a-f - a 2 2 ) to 2 sin a cos a 


- r mu>‘ 


gi (ai 2 ~ a 2 2 ) 

c{a 2 +af] 


By Newton’s third law (§31), the forces exerted on the right and 
left bearings are -R and R, with magnitudes R given in Equation 
(41.18) above. 

Example 3. A gyroscope is mounted so that one point on its axis is fixed. 
Investigate those motions under gravity in which the axis of the gyroscope makes 
a constant angle with the vertical. 

A gyroscope is a body with an axis of symmetry, the shape of the 
body being such that the moment of inertia of the body about its 
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axis of symmetry is large. For example, the disc and shaft in Figure 
49 constitute a gyroscope. 



In Figure 50 the line OA is the axis of the gyroscope, the fixed point 
being at 0 and the center of mass being at C. We introduce a fixed 
unit vector j pointing up from 0, and a set of moving orthogonal 



unit vectors i l9 i 2 and i 3 defined as follows: i 3 is along OA; i x is in the 
plane of j and i 3 , as shown; i 2 completes the triad. We note that i 2 is 
horizontal. 

The angle between OA and j is denoted by 0; it is constant. The 
plane of j and i 3 rotates about j at a rate of p radians per unit time; 
p is called the precession. The gyroscope spins about its axis at the 
rate of s radians per unit time; s is called the spin. We denote the 
mass of the gyroscope by m and the distance OC by l. 

The only external forces are the reaction of the pivot support at 
0 and the gravity forces. The former force has no moment about 0. 
The latter forces may be replaced by a single force -mgj at C. Hence 

G = (l is) X (-mgj). 
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But 

(41.19) j = i x sin 0+i 3 cos 0, 
whence we get 

G = -Img sin 0 i 2 . 

Thus 

(41.20) G x = 0 3 G 2 = -Img sin 0, G 3 = 0. 

The angular velocity «n of the gyroscope is 

w = si 3 +/)j. 

Because of Equation (41.19) we then have 

a) = p sin 0 i 1 +(s'+^ cos 0 ) i 3 , 

whence 

(41.21) cxi-y = p sin 0, co 2 =0, co 3 = s+p cos 0. 

The angular velocity £2 of the coordinate axes is given by the re¬ 
lation 

si = Pi 

— p (i x sin 0 +i 3 cos 0 ). 

Thus 

(41.22) Q 1 =/>sin0, 0 2 = 0, Q 3 =jf;cos0. 

The coordinate axes associated with i ls i 2 and i 3 are principal axes 
of inertia at 0, and we have 

(41.23) Ji=/ 2 , r 1 = r 2 = r 3 = o. 

We now substitute in Euler’s equations (41.8) from Equations 
(41.20), (41.21), (41.22) and (41.23) to obtain the relations 

(41.24) Iyp sin 0=0, 

(41.25) [/ 3 ^+(/ 3 - Iy) p cos 0] p sin 0 = Img sin 0, 

(41.26) / 3 {s-p cos 0) = 0. 

One solution of these equations is 0 = 0. In this case the axis of the 
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gyroscope is vertical, and the gyroscope is said to be “sleeping”. If 
0 is not equal to zero, then Equations (41.24) and (41.26) yield 


p = constant, s = constant, 
and Equation (41.25) takes the form 

(41.27) (. I z - / x ) cos 0 p^+fsp - Img = 0. 

This is a relation among the three constants p, s and 0. Hence it 
appears that we may assign arbitrarily values for two of these constants 
and there will exist a corresponding motion of the top with 0 a con¬ 
stant, provided of course the value of third constant, as obtained from 
Equation (41.27), is real. 

We note from (41.27) that 

( 41 . 28 ) !sl - < w .) tsal . 

hP ■‘3 


The quantities p and 0 may be observed readily. The corresponding 
spin s may be computed by means of this relation. If the precession 
is small, we note from Equation (41.28) that the spin is large and has 
the approximate value 



42. The general motion of a rigid body. We now consider a rigid body 
moving in a general manner. It may or may not have a fixed point. 
The motion of its mass center can be determined from Theorem 2 of 
§ 40, which applies to the motion of any system of particles. This 
theorem yields 

(42.!) *^ = F, 

where m is the total mass of the body, v c is the velocity of its center of 
mass, and F is the sum of the external forces acting on the body. 
Integration of (42.1) gives the position-vector x c of the center of mass 
C of the body as a function of the time t. 

To determine the complete motion of the body it is then only nec¬ 
essary to find the angular velocity of the body about its center of mass. 
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To do this, we choose the origin 0 of the coordinate system at the 
center of mass C of the body. We then consider the body as having 
a velocity of translation v c plus an angular velocity co with a line of 
action through C. The velocity v ; - of the j-th particle is then given by 
the relations 

(42.2) vy = v c+^ = vc+gj X xj. 

But by definition the angular momentum h of the body about the 
point 0 is 

JV 

h = 2 rrijXjX-Vj, 
j = i 

and we then have, just as in § 41, 

(«.3) f= A +». 


where 


JV J„. 

(42.4) A =S«/| X v i , 


3 -1 


Jf dv 

B = S myXy x ^ . 

i=l at 


From Equation (42.2) we then have 

N 

A = 2 OTy(Vy - V c ) X Vy 

j = 1 

X JV 

= E MyVy X Vy - V c X 2 OTy Vy 

J = 1 j = 1 

= 0 - v c xM, 

where M is the linear momentum of the body. But M = mv c by 
Equation (39.2). Thus A = 0. Just as in § 41 we find that B = G, 
where G denotes the total moment about 0 of all the external forces. 
Thus Equation (42.3) takes the form 


(42.5) 


dh 

dt 


= G . 


We have thus the result: the rate of change of the angular momentum of 
a body about its center of mass is equal to the total moment of the external 
forces about the center of mass. 
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We have placed the origin of the coordinate system at the center of 
mass C of the body. If we choose the coordinate axes to coincide with 
principal axes of inertia of the body at C, then just as in § 41 we obtain 
Equations (41.6) and finally Euler’s Equations (41.8) from which we 
can find the unknown quantities w 1; co 2 and co 3 which characterize the 
rotation of the body about its center of mass. 

In conclusion, it should be noted particularly that the equation 
dh/dt = G can be used only in the two following cases: (i) the body 
has a fixed point and the origin is at this fixed point; (ii) the origin 
is at the center of mass of the body. 

Example 1. A gyroscope with a constant spin is carried along a horizontal 
circular path at a constant speed , with its axis tangent to the path of its center 
of mass. Find the forces exerted on the axle of the gyroscope by the bearings in 
which the axle turns, neglecting gravity. 

We choose the center of mass of the gyroscope as the origin 0 of the 
coordinate system. The path of 0 is shown in Figure 51; it is a circle 



with center A and radius a. We introduce an orthogonal right triad 
of unit vectors at 0, as shown; ij points vertically up; i 2 points towards 
A ; i 3 is tangent to the circle and hence lies along the axis of the gyros¬ 
cope. 

Let us suppose that 0 travels at a speed v in the direction opposite 
to i 3 , and that s is the rate at which the gyroscope spins about its axis, 
s being positive when it is in the sense of the 90 rotation from i, 
to i 2 . 

The time required for 0 to go around the circle is 2 7c ajv. In this 
time the gyroscope has turned about through an angle of 2 tc 
radians. Hence the angular velocity to of the gyroscope is 
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Thus 


“> =si 3 + 


2tt . 
2n a\v 1;L 


v . , . 


(42.6) <o x = - > co 2 = 0 , « 3 = s . 

The angular velocity £2 of the coordinate axes, is 




o v • 

S2 = -1, 

a 

whence 



(42.7) 


n 2 = o, q 3 = o. 

Also 



(42.8) 

II 

o 

II 

k" 

II 

II 


We now substitute in Euler’s equations (41.8) from Equations 
(42.6), (42.7) and (42.8) to obtain the relations 

C I= 0. c, = o. 

Hence the forces exerted on the gyroscope by the bearings must be 
in the x 3 x x plane. If R and -R are these forces, and 2 b is the distance 
between the bearings, then 

2 bR =I 3 - 

3 a 

whence 


Problems 

1. A particle moves on the curve x 2 = k tan kx l3 x 3 = 0, where h 
and k are constants. The x 2 component of the velocity is constant. 
Find the acceleration. 
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2. A particle moves with constant speed. Prove that its acceleration 
is perpendicular to its velocity. 

3. A particle moves on an elliptical path with constant speed. At 
what points is the magnitude of its acceleration (i) a maximum, (ii) 
a minimum? 

4. A particle moves in space. Find the components of its velocity and 
acceleration along the parametric lines of spherical polar coordinates. 

5. A particle moves in space. Its position-vector x relative to the 
origin of a fixed set of rectangular cartesian coordinate axes is given in 
terms of the time t by the relation 

x = cos f+i 2 sin /+ i 3 1 ), 

where h is a constant and i 1; i 2 and i 3 are the usual unit vectors in the 
directions of the coordinate axes. Find the components of the veloc¬ 
ity and acceleration in the directions of (i) the coordinate axes men¬ 
tioned above, (ii) the principal triad of the path of the particle, (iii) 
the parametric lines of spherical polar coordinates. Find the speed 
and the magnitude of the acceleration. 

6. A particle describes a rhumb line on a sphere in such a way that 
its longitude increases uniformly. Prove that the resultant acceleration 
varies as the cosine of the latitude, and that its direction makes with 
the inner normal an angle equal to the latitude. 

7. Two forces A and B act at a point. If « is the angle between 
their lines of action, prove that the magnitude of the resultant R is 
given by the relation 

R 2 = A 2 +-B 2 +2A6 cos a. 

8. Four forces A, B, C and D act at a point 0 and are in equilib¬ 
rium, the forces C andD being perpendicular and having equal mag¬ 
nitudes. Find C in terms of A, B and the angle a between A and B. 

9. Forces with magnitudes 1, 4, 4 and 2V3 lb. wt. act at a point. 
The directions of the first three forces are respectively the directions 
of the positive axes of x 1 , 'x 2 and x 3 . The direction of the fourth force 
makes equal acute angles with these axes. Find the magnitude and 
direction of the resultant. 
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10. A force F acts on a particle of mass m. Find the magnitude of 
the acceleration, given that (i) F = 6 poundals, m = 3 lb., (ii) F = 

6 lb. wt., m = 2 slugs, (iii) F = 6 poundals, m = 2 slugs, (iv) F = 

6 lb. wt., m — 3 lb., (v) F = 5 dynes, m = 10 gm. 

11. A particle of mass m is acted upon by a force F given by the 
relation 

F = 16p sin 2f+qe'‘, 

where p and q are constant vectors and t is the time. Find the veloc¬ 
ity v and position-vector x of the particle in terms of t, given that v 
= 0 and x = 0 when t = 0. 

12. A particle of mass m is acted upon by two forces P and Q,. 
The force P acts in the direction of the x 1 axis. The force Q, makes 
angles of 45° with the axes of x 2 and x z . Also P — p sin kt and CL = 
q cos kt , where p, q and k are constants and t is the time. At time 
t — 0 the particle has coordinates ( b , 0, 0) and is moving towards the 
origin with a speed p/mk. Find the position-vector x of the particle. 
Prove that the particle moves on an ellipse, and find the center and 
lengths of the axes of the ellipse. 

13. A particle of mass m moves under the action of a force pe' qt 
and a resistance -Zv, where p is a constant vector, q and Z are positive 
constants, t is the time, and v is the velocity of the particle. Prove that 

Xoo- X„ — J~q (p+ m ? u ) 

where u is the velocity when t = 0, and x 0 and x M are respectively the 
position-vectors of the particle when t = 0 and when t becomes in¬ 
finite. Is the above result true when Z = mq ? 

14. A particle of mass m moves under the action of a force p cos qt 
- kx, where p is a constant vector, q and k are positive constants, 
t is the time, and x is the position-vector of the particle relative to a 
fixed point 0. When t — 0, the particle is at 0 and has a velocity u. 
Find x in terms of t when (i) k ^ mq 2 , (ii) k = mq 2 . 

15. A particle of mass m is acted upon by a single force y mjx 2 
directed towards a fixed point 0, where y is a constant and x is the 
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distance from 0 to the particle. At time t = 0 the particle is at a point 
B and has a velocity of magnitude u in a direction perpendicular to the 
line OB. Prove that the orbit is (i) an ellipse if bu 2 < 2y, (ii) a parab¬ 
ola if bu 2 = 2y, (iii) an hyperbola if bu 2 > 2y, where b = OB. 

16. Find the moment of inertia of a circular disk of mass m and 
radius a about (i) the axis of the disk, (ii) a diameter of the disk. 
[Answer: (i) -|ma 2 ; (ii) \ma 2 .'] 

17. Using the result of Problem 16, find the moment of inertia of 

a circular cylinder of mass m, length 21 and radius a about (i) the 
axis of the cylinder, (ii) a generator of the cylinder, (iii) a line 
through the center of the cylinder perpendicular to its axis, (iv) a 
diameter of one end of the cylinder. [Answer: (i) (ii) |ttkz 2 ; 

(iii) j 2 m (4/ 2 +3a 2 ); (iv) ^ m (16/ 2 +3a 2 ).] 

18. A circular cylinder has a mass m , length 21 and radius a. Rec¬ 
tangular cartesian coordinates are introduced, with origin 0 at the 
center of the cylinder, and the x 3 axis coinciding with the axis of the 
cylinder. Two particles each of mass m' are attached to the cylinder 
at the points (0, a, l ) and (0, -a, -1). Find the moments and products 
of inertia / 1; / 2 , / 3 , K 1} K 2 , and K 3 for the system consisting of the 
cylinder and the two particles. 

19. A circular disk of mass m and radius a spins with angular speed 
w about a line through its center 0, making an angle a with its axis. 
Find the angular momentum of the disk about 0. 

20. For the system of masses in Problem 18, find the angular 
momentum about the point 0 when the system has an angular speed 
co about (i) the x 1 axis, (ii) the x 2 axis, (iii) the x 3 axis. 

21. A circular cylinder of mass m, length 21 and radius a turns 
freely about its axis which is horizontal. A light inextensible cord is 
wrapped around the cylinder several times. A constant force F is 
applied to the end of the cord. If the cylinder starts from rest at time 
t = 0, show that at time t it has turned through the angle Ft 2 /ma. 

22. The circular cylinder of Problem 21 is again mounted with its 
axis horizontal, and has a light inextensible cord wrapped around it. 
A body with a mass m' is attached to the end of the cord. If the cylinder 
starts from rest at time t = 0 show that at time t the cylinder has 
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turned through the angle 
to gravity. 


mgt 2 

a{m J r2m ) 


where g is the acceleration due 


23. A circular cylinder of mass m, length 21 and radius a is placed 
on a rough plane which makes an angle a with the horizontal, and is 
then released. Find the distance the cylinder moves down the plane 
in time t. 

24. In Problem 18 there was introduced a system consisting of a 
circular cylinder of mass m with two attached particles each of mass 
m'. This system is mounted so it can turn about the axis of the cylinder 
in two smooth bearings each at a distance c from the center of the cy¬ 
linder. The system is made to rotate with constant angular speed <o. 
Find the reactions of the bearings. 

25. A circular disk of mass m and radius a turns with constant angu¬ 
lar speed <jj about an axis through the center 0 of the disk and making a 
constant angle a with the axis of the disk. The disk turns in two 
smooth bearings each at a distance c from the point 0. Find the reac¬ 
tions of the bearings. 

26. A uniform rod of length 21 is free to turn about an axis L 
perpendicular to it and through ics center. The center of the rod 
moves at constant speed v around a circular track of radius a, the 
axis L being always tangent to the track. Deduce the equations of 
motion of system. 
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Chapter IV 


PARTIAL DIFFERENTIATION 

43. Scalar and vector fields. Let V denote a region in space, and let 
A be a general point in V. Let x 1 , x 2 , x 3 denote the rectangular cart¬ 
esian coordinates of X, and let x denote the position-vector of X. Then 

* = x^+x^+x^. 

Let us now consider the case when there is associated with each 
point in the region V a value of a scalar f. Then we write/ = f (x l , 
x 2 , x 3 ) , or more compactly 

/=/(*)• 

The values of/associated with all the points in V constitute a scalar 
field. 

Let us now consider the case when there is associated with each 
point in the region V a value of a vector a. Then we write a = a(x 1 , 
x 2 s *s) } or more compactly 

a = a(x). 

The values of a associated with the points in V constitute a vector field. 

It is frequently necessary to consider scalar and vector fields which 
vary with a parameter, such as the time t. In such cases we write 

/=/(*, a = a(x, t). 

44. Directional derivatives. The operator del. Let C be a curve in a 
region V, as shown in Figure 52. Let A be a general point on C, with 
position-vector x as shown, and let s be the arc length of C measured 

from some fixed point Q on C. It was seen in § 28 that the vector ~ 

ds 

is a unit vector tangent co C in the direction of s' increasing. This 
vector is shown in Figure 52. 
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Figure 52 


We now introduce a function f [x t , x 2 , x 3 ) defined everywhere in 
the region V. Then f is defined at all points on the curve C, and the 
rate of change of f with respect to the arc length s at the general point 
X on C is given by the well-known formula of differential calculus, 


(4 


• 1 ) 


df __ 3f_ dxy ( dj^dx 2 , df dx 2 


ds 


_ 1 “ _i_ _3. 

dx 1 ds dx 2 ds dx 3 ds 


Let us now consider all curves through X having the same tangent at 
X as C has there. For each of these curves the value at X of the right 
side of Equation (44.1) is the same. Thus at each point X there is 

a unique value of -j- associated with each direction. This value is 

called the directional derivative of f. 

We now introduce an operator called del, which we denote by the 
symbol V and define by the relation 


„ • 3 . • 3 

V = 11 ^ + * 2 Yx 2 


d_ 

dx„ 


(44.2)- V =i i^T +1 ^ 2 + >3 „ A3 

If we wish to operate on a function with this operator, we write the 
symbol denoting the function immediately to the right of the symbol 
V. Thus, if/is a scalar field, then 

d , . d 

12 3x, 


or 


V/ = 

V/= ix 


+ *2 “ + *3 —V 


1 dx 1 

3L + i ^ +i 

dx x 2 dx 2 3 


13 ix 3 y 
3J_ 

dx. 
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The expression V/ is often called the gradient of f and is denoted 
by grad /. 

We note that the right side of Equation (44.1) is equal to the scalar 
product of V/and dxjds, by Equation (7.2). Thus 


(44.3) 

*1 

& 

II 

^143 

or x 


(44.4) 


where t is the unit vector in the direction in which the directional 
derivative is taken. Some theorems will now be proved. 

Theorem 1. The component of V/ in the direction of a unit vector 
t is equal to the directional derivative of f in that direction. 

Proof. The component of V/ in the direction of t is | V/ 1 cos 0, 
where 0 is the angle between the vectors V/ and t. But, since t is a 
unit vector, 

|V/| cos 0 = V/-t 
or 

(44.5) |V/|cos0=J^> 

by Equation (44.4). 

Theorem 2. The vector V/ points in the direction in which dfjds 
has a maximum value; also, this maximum value is equal to | V/|. 

Proof. Both parts of this theorem follow from Equation (44.5) 
which shows that when 0=0, dfjds attains its maximum value which 
is |V/|. 

Theorem 3. The vector field V/ is normal to the surfaces f = 
constant. 

Proof. Let X be a point on a surface f = constant, and let t be 
a unit vector tangent at X to this surface. Then the value of dfjds at 
X corresponding to the direction of t vanishes. By Equation (44.4) 
we then have 
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V/t = 0. 



This implies that V/ is perpendicular to all vectors at X tangent 
there to the surface f = constant, which completes the proof. 


45. Properties of the operator del. This operator, which is denoted 
by the symbol V, is defined in the previous section. For convenience 
we express it in the form 

< 45 -» 

Then 

< 45 - 2 > 

We shall now prove a number of theorems relating to this operator. 
Theorem 1. If/ and g are scalar fields, then 


(45.3) V (f+g) = V/+V£. 

Proof. Now, by Equation (45.2) we have 


r= 1 dx, r =l 

= V/+ Vg. 


dx, 


Theorem 2. If f is a function of a single scalar field u , then 
(45.4) 




Proof. By Equation (45.2), we have 

V/ = S ir^-' 


r=l dx r 

But/is a function of u, which is a function of x,. Thus 

Qf_ _ df du_ 
dx, du dx. 


(r= 1,2, 3), 


so 
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j 4[2L 

Zx'dudx, 

= # | . £« 
du r=\ ’ dx, 

= fvu- 

du 

Theorem 3. If/is a function of/z scalar fields z/ l5 h 2 > * • *, #», then 

(45.5) V/_^V« 1+ f i V«, + ...+f_V,, 

=i / vn,. 

Proof. Since/is a function of the n variables u l9 u 2 , • • •, «„ which 
are themselves functions of x 1 , x 2 , x 3 then 


dj_ = j. 

dx r j = i du s dx r 

Hence 

3 n df du 

V/=2 i, 2 f-p 

r=l 4=1 du, dx, 

= v y J 
j=i r=i r du, dx V 

v;4 

-^du.rZl^dKr 

= 2 fv*. 

4=1 du. 

As mentioned in the previous section, V operates only on functions 
written on its immediate right. Thus, if/and g are scalar fields, then 
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Theorem 4. If / and g are scalar fields, then 
(45.6) V (fg) =fVg+gVf. 

Proof. We have 

=U(/i 


f= 1 r y / dx r 5 <9* r 
=/V^+ iV/. 


46. Some additional operators. Let a be any vector field. We first 
consider the operator a-V- This operator has the obvious meaning 

(46.1) a-V = (a 1 i 1 +a 2 i 2 +a 3 i 3 )-(i 1 |- + i 2 + i 3 J-\ 


fll dx x ^ d ‘ l dx 2 
v 5 

= -T-- 


a* x ' 2 Px 2 

d_ 

dx 3 


r= 1 f 

This operator is a scalar, and can be applied to a scalar field or to a 
vector field. Thus, if/ and b are two fields, then 

(46.2) (»-V)/ = |,«,f r , (-V)b=lig. 

We note that 

(a-V)/=a-V/. 

The operator ax V can be considered similarly. We have 

(46.3) a X V = Kii+flaLj+agi;,) X f i x -j- + i 2 + i 3 ’ 


a ^- a3 ^) + 4 a ^ r al ^ 3 . 
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(46.5) 


ax V = 


• • • 

*1 *2 *3 

d-y &2 a 'i 

d_ l_ 3_ 

dx-y dx 2 dx 3 

The operator a X V is a vector operator. It can be applied to scalar 
fields. Thus, if/ is a scalar field, then 


(46.6) (a x V)/ = 

We note that 


a i 

£f 

dXy 


*2 

dx 2 


*3 

a 3 

dx 3 


(a X V) / = a X V/- 


In writing the expressions (a-V)/, (a-V)b and (a XV)/, one must 
exercise care in the matter of the order in which the symbols appear, 
since the operator V and all operators constructed from it operate 
only on the functions on their immediate right. Thus, for example. 


(a-V)b ^b(a-V). 


The left side of this expression is a vector field, while the right side 
is a vector operator. 

We now introduce the operators V • and V X. These operators can 
be applied to vector fields, the vector portion of the operator V 
operating on the vector field with scalar or vector multiplication. If 
b is a vector field, then 


(46.7) ™ = (i.|; + *4, + i -Jr»)' b - 

Thus we may write 

(46 - 8 > v b - (A* k} h =A b 's - kk < b ' b >- 

But 

Vb =b r (r = 1, 2, 3). 

Thus 

(46.9) V b = S ^ +^- 2 + ^- 3 . 

r= i ox r dx x dx 2 dx 3 
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The expression V ■ b is often called the divergence of b , and is 
written div b. 

We also have 

(46.10) Vxb-(i,i. + 1 ,^ + l,^xb. 

Thus we may write 

(46.11) Vxb = (UA) xb - i i,x g = (txb). 
But 

iiXb =iiX(i 1 i 1 + i 2 i 2 +A 3 i 3 ) 

= ^ 2*3 ~ ^ 3*2 

by Equations (8.9). Similarly we have 

i 2 xb = bfo-b x i 3 , i 3 xb=J 1 i 2 -i 2 ii. 

Hence (46.11) becomes 

/ 3( ~ 19') V7vli _i ^ 2 \ _i_i l^i 3b 3 \ . I3b 2 3 b A 

(46.12) Vxb-.,^- + i, ^ ^ - T \ ■ 


4 v -- dx3 j T *1 ^ ^ 

This can be written conveniently in the form 


Vxb = 


(46.13) Vxb= ~ E it- . 

C/ xV C/ A 2 C/ A 3 

^1 ^2 ^3 

The expression V X b is often called the curl of b or the rotation of 
b , and is written curl b or rot b . 

Theorem. 1. If a and b are two vector fields, then 

(46.14) V • (a+b) = V ■ a+V • b, 

(46.15) Vx(a+b) = Vxa+Vxb. 

Proof of Equation (46.14). From Equation (46.9), we have 

V-(a+b) =1 ^(Or+br) 

4 (t+t) 

= V- a-j-V-b. 
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Proof of Equation (46.15). From (46.13) we have 

*2 >3 

, d d d 

V X (a+b) - 3 d 


»i »2 

*3 

3 A 




a l + ^l a 2"t~^2 

«3 + ^3 

ii i 2 

»3 

i. -L 

d_ 

dx x dx 3 

dx 3 

a x a 2 

a 3 

ii 

i 2 

d_ 


dx x 

dx, 

K 

b 2 


= Vxa+Vxb. 

Theorem 2. If a and b are two vector fields, then 
(46.16) (axV)'b = a - (V xb). 

Proof. Because of (46.4-) we have 


(46.17) 


(axV)-b = U 2 


d d X,. , , 
dL * 3 ftj (vb) 


+ | fl3 h~ dl i) (ia ‘ b) 




ij-b = b x , i 2 -b = b 2 , i 3 -h = b 3 . 

Thus (46.17) can be written in the form 

l3b^\ (db x db 3 \ I db^ ~ db-\ _ 

(a X V) -b = - *-) + - aj + "3 at.) 

Because of (46.12), we see that the right side of this equation is 
equal to a - (V xb). 


(i 3 ‘b) . 
i 3 -b = b a . 


!db 3 db£\ I 3b x db 3 \ / dbi _ 
(axV)-b = a 1 [ 3 ~- 3 --)+a 2 l^--- 2> ) + a 3 l ? 
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47. Ivariance of the operator del. Let us consider two rectangular 
cartesian coordinate systems with a common origin 0. We denote 
them by the symbols S and S'. For the coordinate system S the coor¬ 
dinates and associated unit vectors will be denoted by x 1 , x 2 , x 3 and 
i x , * 2 ) * 3 ) respectively, while for the system S' the corresponding 



Figure 53 


quantities will be denoted by xf, x 2 ', x 3 ' and if, i 2 ', i 3 '. Figure 53 
shows the coordinate axes and the associated unit vectors. 

We now introduce the two operators 


Y7 


V 


~ 11 dx 1 + 12 dx 2 + 13 dx 3 


V'= 


*1' 





Let /be a scalar field, and let b be a vector field. We shall now con¬ 
sider proofs of the relations 


(47.1) 

< 

H 

II 

(47.2) 

V' . b = V • b, 

(47.3) 

V'Xb — VXb. 


The implications of these equations are the following. The operator 
del includes differentiations with respect to the coordinates of the 
system S and multiplications by the three unit vectors associated 
with the system S. It thus appears that when del operates on any 
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field, the resultant field is likely to depend on the particular choice 
of coordinates. However, Equations (47.1), (47.2) and (47.3) state 
that this is not the case. It is this property of del which is termed 
“invariance of the operator del”. 

Before proceeding to a proof of the above three formulas, we shall 
develop some preliminary formulas. Let us consider the three angles 
which the x x ' axis makes with the axes of the system S. These are the 
angles 0 U , 0 12 and 0 13 in Figure 53. We denote their cosines by a 11; 
a 12 and a 13 . These three quantities are the components of i x ' relative 
to the system S, and so 

3 

*i = a n *i~f" a i2 13 *3 = 21 a. ls v 

S = 1 

There are similar expressions for i 2 ' and i 3 '. The three expressions 
can be written in the compact form 

(47.4) i/=S«„i, (r = 1, 2, 3). 

s = 1 

We note that the constant a rs is the cosine of the angle between the 
two vectors i/ and i,. 

The components of the vector i 3 relative to the system S' are seen 
to be a n , a 21 and a 31 . Thus 

3 

ii = 2 i/. 

s= 1 

Proceeding similarly for i 2 and i 3 , we obtain the three equations 

(47.5) ir = E a„ i/ (r = 1,2,3). 

s = 1 

Throughout the remainder of this section we will adopt the conven¬ 
tion that latin subscripts range over the values 1, 2 and 3, as in the 
above equations. 

We now introduce a set of nine quantities called the Kronecker 
delta. This set is denoted by the symbol S rs , and is defined by the 
equation 

(47.6) 
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S„ = 1 if r = s 
— 0 if r =2= s. 



Thus Sji — §22 — ^33 — 1 and §23 — 8 31 — 8 12 — 8 32 — 8 13 — ^21 — 0. 
Since i x , i 2 and i 3 are unit orthogonal vectors, they satisfy the nine 
relations (7.3). These nine relations can now be written compactly 
in the form 

(47.7) i, i , = S„. 

We also have the nine relations 


(47.8) i/i/ = 8„. 

Let us now substitute in Equation (47.8) from (47.4). This gives 
the relation 

8„ = L) • 

3 3 

= S £ a rt a m i,-i a . 

1 = 1 U = 1 


But by (47.7) we have i,-i u = §,„ so 
3 3 


(47.9) 


§„= S £ a rt a su K 

' = 1 U = 1 


3 3 

= £ a rt £ a m 8 lu . 

1=1 U =1 


The last summation here is 
3 

E d su 8 tu = d s i ^2 8,2 +<2,3 §,3. 

U = 1 

The right side reduces to a sl when t = 1, to a s2 when t= 2, and to 
a s3 when t = 3. Thus we may write 

3 

(47.10) £fl IU S,„=a J( , 


and so (47.9) reduces to the relations 

(47.11) 8 „ = £ a rt a st . 

< = 1 

The nine equations in (47.1.1) are called the orthogonality conditions. 
We could also run through the above derivation of (47.11) but with 
the roles of the primed and unprimed quantities interchanged. This 
would entail substituting in Equations (47.7) from (47.5). In this 
way we would obtain the orthogonality conditions in the form 
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(47.12) 


8 „ — 2 o-t, o-is • 
i = 1 


Let b be any vector field with components b T relative to the system 
S, and with components b r ' relative to the system S'. Then 


(47.13) 


3 3 

b = £ Mr = 2 b/i/. 

r — 1 s — 1 


We now substitute here for i r from (47.5), obtaining the relation 

3 3 3 

£ b; i,' = £ b, £ a„ i/ 

r = 1 r= 1 r = 1 

3 3 

= £ hb r a„ i/. 

.1 = 1 r = l 


Hence we must have 


(47.14) 


bj = £ b r . 

r = 1 


If in an analogous fashion we substitute in Equation (47.13) for i/ 
from (47.4), we obtain the relations 


(47.15) 


b s = £ b/. 

r = l 


Equations (47.14) and (47.15) are the equations of transformation 
of the components of a vector field b. 

If we choose b = x, where as usual x is the position-vector of a 
general point X with coordinates x r relative to the system S and 
coordinates x/ relative to the system S,' then Equations (47.14) and 

(47.15) yield 

3 3 

(47.16) x/ =T,a„x r , x s = £ a„ x /. 

r =1 r=l 

These are the equations of transformation from a set of rectangular 
cartesian coordinates to a second set whose axes are obtained from 
those of the first set by a rotation about the origin. From Equations 

(47.16) we have the relations 
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(47.17) 


We are now in a position to turn to proofs of Equations (47.1), (47.2) 
and (47.3). 

Proof of Equation (47.1). This equation reads 

vy= v/. 


(47.18) 


v/= S i',i£ 


dx" i = i ibc, dx' r 

Because of this relation and (47.4) we can write (47.18) in the form 

3 3 3 df dr 

(47.19) V/ = S S Za„i,2-%L. 


,= i j=i 1 = 1 


dx, dx'. 


But, from Equations (47.17) we have 


dx'r - a,t • 


Thus (47.19) becomes 


3 3 3 

vy = s S UUt 

r= 1 j= 1 / = 1 

33 df 3 

= 2 S E 

j=l / = 1 oXi r = i 

Because of the orthogonality conditions (47.12), the last summation 
here is equal to 8„. Thus 

3 3 df 
vy=s Si, ^8,, 

,= i i=i dx, 

3 3 df 

= S i, S f - 8„ . 
j=i i=i dx. 

The last summation here reduces to dfjdx 1 when s = 1, to dfjdx i 
when s = 2, and to dfjdx 3 when s = 3. Thus 


115 



(47.20) 
and so 


3 

2 


<=i 


dx. 


8 * 


dx,' 


v/= 


3 

= £ i. 


& 


= V/. 


The truth of Equation (47.1) also follows from Theorem 2 of §44, 
as this theorem states that V/ points in the direction in which the 
directional derivative of/has a maximum value, and that the magni¬ 
tude of V/ is equal to this maximum value. This theorem then im¬ 
plies that both the direction and magnitude of V/ are independent 
of the coordinate system. 

Proof of Equation (47.2). This equation reads 


Now 


V'-b = V-b. 


V' 


•b=l|V 

s=i dx s 


Because of (47-14) we then have 


3 3 T 

V'-b = 2 2 -fy ( a„b T ) 

s=I r=l dx , 

V V 3b ' 

= 2 2 a„ -=7 

*=1 r= 1 dX , 

-V y y tor dj±_ 

~ -I r-t "l-l dx, dx', ' 

But by Equations (47.17), the last partial derivative here is equal to 
a,,. Thus 

3 3 3 db 

V'-b = 2 2 2 a„a«^ 

,= 1 r=l 1 = 1 dx, 

3 3 db 3 

= 2 2 ^2 a„a». 

r- 1 1 = 1 dX,,= 1 


Because of the orthogonality conditions (47.12) the last summation 
here reduces to 8„. Thus 
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3 3 db 

V'-b =S S ^-'S r , 

r-l 1 = 1 dX t 

l d A 

,-I a*’ 

by a repetition of the arguments leading up to Equations (47.10) and 
(47.20). This completes the proof. 

Proof of Equation (47.3). This is left as an exercise for the reader 
(Problem 13 at the end of the present chapter). 

48. Differentiation formulas. We shall consider here a group of 
well-known formulas involving the operator del. If f is a scalar field, 
a and b are vector fields, and x is the usual position-vector of a gen¬ 
eral point X, these formulas are the following: 


(48.1) 

V-(/a) =. 

/(V-a) + (V/)-a, 

(48.2) 

VX(/a) = 

=/(Vxa) + (V/)Xa, 

(48.3) 

V-(aXb) 

= b-(Vxa)-a-(Vxb), 

(48.4) 

VX (axb’ 

) =a(V-b) + (b-V)a-b(V-a)-(a-V)b, 

(48.5) . 

II 

S' 

A 

i> 

(a-V)b+(b-V)a+a < (V Xb) +b X (V Xa), 

(48.6) 

V x (V/) = 

= 0, 

(48.7) 

V-(VXa) 

= 0, 

(48.8) 

V X (V Xa 

) = v (v-a) - (v-v) a, 

(48.9) 

V-x =3, 


(48.10) 

VXx = 0 


(48.11) 

(a-V)x = 

a. 


Direct proofs of all eleven of these formulas follow similar lines. We 
shall present here only proofs of (48.1), (48.3), (48.5) and (48.8). 
The proofs of the remaining formulas are left as exercises for the 
reader (Problems 14, 15 and 16 at the end of the present chapter). 

Proof of Equation (48.1). Because of Equation (46.8), which gives 
some equivalent forms for V • b, we have 
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V-(/a)=! i r ~(/a) 


, 3f 


= Ei f ./^ + 


r= I f \ J 1 a*,y 

_ j-y . <?» , V • d .f 

r= 1 OX r r= 1 <7X r 

=/(V-a) + (V/)-a. 

Proof of Equation (48.3). Because of Equation (46.8) we have 
V-(axb) =2^,— (axb) 

= 2 v^xbUs vfaX^V 


~r\dxr~) \zr\ ''dx r j 

Because of the permutation theorem for scalar triple products (Theo¬ 
rem 1 of § 9), we can write this last equation in the form 

r= 1 OX, r=l OX, 

= b-(Vxa) - a-(Vxb). 

Proof of Equation (48.5). Because of Equation (46.11) we have 

3 

aX(Vxb) = ax2 i r X^— 

r-1 OX, 

We now apply to the right side of this equation the identity (9.3) 
for vector triple products, obtaining 

(48.12) ax(Vxb) =2 ifa—U 2 ^ (a-i,) 


dx,j r= I OX, 


= 2 i r ( / a-^\-(a-V)b. 
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Similarly, by an interchange of a and b we have 

(48.13) bx(VXa)=Si r |^-b|-(b-V) a. 

Addition of (48.12) and (48.13) yields 

3 d 

ax (Vxb)+bx(Vxa) =S i, -f- (a-b) - (a-V) b - (b-Va) 

r=l OX r 

= V (a-b) — (a-V)b - (b- V) a. 

If the last two terms here are taken to the left side of the equation, 
(48.5) results. 

In addition to this manner of proving formulas (48.1)-(48.11), 
there is another manner of proof which applies to some of these for¬ 
mulas, and which is quite expeditious. We mention it here because it 
affords one an opportunity to acquire additional facility in manipu¬ 
lations involving the operator del. This manner of proof consists in 
applying to expressions involving del the permutation theorem for 
scalar triple products (Theorem 1 of § 9), and the identity (9.3) for 
vector triple products which is 

(48.14) aX(bXc) = b(a-c) - c(a-b). 

Of course, when such an application is made, a symbol denoting a 
scalar or vector field must never be moved from one side of del to the 
other, since del operates only on those quantities on its right. Also, 
when such application is made, the resulting expression must have the 
same scalar or vector character as the original expression. 

For example, let us prove Equation (48.8) in this way. If we apply 
the identity (48.14) formally to V X (V Xa), we obtain the difference 
of two terms. The first term must be one of the following: 

(48.15) V(V-a), V(a-V), (V-a) V, (a-V)V. 

Now, in the expression VX(VXa), the symbol a appears to the 
right of both symbols V • Hence we must select from the terms in 
(48.15) the one with this same property. Hence we select V(V - a). 
As a check, we note that V X (V Xa) is a vector, and that the only 
term in (48.15) which is a vector is V (V-a). In a similar fashion 
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we find that the term corresponding to the second term on the right 
side of (48.14) is (V-V) a. Thus 


Vx(Va) = V (V-a) - (V-V)a. 

We note that the operator V • V satisfies the relation 

V r, * * * 

v v ~ (ao* + (^ 2 ) 2+ 

This operator is called the Laplacian operator, and is denoted in 
many books by the symbol V 2 - 

As a second example, let us prove Equation (48.3) by application 
of the permutation theorem for scalar triple products. To do this, 
it is convenient to introduce an operator V„ which is just the operator 
del with the added restriction that the partical differentiation is to be 
applied to no vector field other than a. The operator V* is defined 
similarly. Equation (12.4) states the differentiation rule for vector 
products. Because of this rule we have 

(48.16) V-(aXb) = V*-(axb)+V 4 -(aXb). 

We now apply the permutation theorem for scalar triple products 
to the first term on the right side of Equation (48.16), remembering 
that a must always remain on the right side of \7 a , and that b need 
not do so. We then get b- (V„ X a) . If the second term on the right 
side of Equation (48.16) is treated analogously, we obtain the relation 


But 


V-(axb) = b- (VjXa) — a - (V&Xb). 


Vo X a = V X a, VjXb=VXb. 


Thus Equation (48.3) is proved. 


49. Curvilinear coordinates. Let us write down the three equations 

(49.1) Zl =A (*!, *2, X 3 ), ' *2 =/ 2 (*!, *2, *3) > ^3 =/s(*l, X 2, X i) , 
where/x (* x , x 2 , x 3 ), f 2 (x 1 , x 2 , x 3 ) and/ 3 (x 1 , x 2 , x 3 ) are any functions 
which are single valued and differentiable throughout some region V. 
These equations prescribe for any point X with rectangular cartesian 
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coordinates , x 2 and x 3 a new set of coordinates Z\ , Z 2 and z 3 . These 
new coordinates are called curvilinear coordinates, and Equations 
(49.1) are equations of transformation of coordinates. 

We now propose to compute, in terms of quantities pertaining to 
curvilinear coordinates only, the expressions V/ V • b and V X b, 
where/and b denote respectively a scalar and a vector field. With 
this goal in mind, we shall devote the rest of this section to some pre¬ 
liminary considerations, and shall complete the final computations 
in the next section. 

The Jacobian of the transformation (49.1) is the determinant 


2zi 

2zi 

2z x 

dx x 

dx 2 

dx 3 

3z 2 

dz 2 

2z 2 

dx x 

dx 2 

dx 3 

2z 3 

2 z% 

2z 3 

dx x 

dx 2 

dx 3 


We shall consider only the case when /' does not vanish anywhere 
in V, so that Equations (49.1) may be solved 1 for x x , x 2 and x 3 to 
yield the relations 

(49.3) x 1 =g 1 {z 1 ,Z 2 ,Z 3 ), * 2 = gt(Zi,Z t ,Za), x 3 =g 3 {z 1 ,Z 2 ,Z 3 ). 

Let us write down the equation 

Zi =fi ( X 1 > x 2 , x 3 ) = constant. 

The locus of this equation is a surface. If we vary the constant in this 
equation, we get a family of surfaces called the parametric surfaces 
of Z x - Similarly, we have parametric surfaces of z 2 and z 3 - Let us 
now consider the curves of intersection of the parametric surfaces of 
z 2 and z 3 . Along each of these curves z 2 and z 3 are constant, and z x 
alone varies. These curves are called the parametric lines of z x ■ Simi- 

1 See almost any book on Advanced Calculus; for example, I. S. Sokolnikoff, 
Advanced Calculus, McGraw-Hill Book Co., New York, 1939, pp. 430-438. 
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larly, we have parametric lines of z 2 and z 3 ■ For example, in the case 
of cylindrical coordinates r, 0, x 3 , the parametric lines of r, 0 and x 3 
are respectively horizontal straight lines cutting the x 3 axis, horizontal 
circles with centers on the x 3 axis, and vertical straight lines. 

Let X be a general point in a region V. Through X there passes a 
parametric line of each of the curvilinear coordinates z 1; Z 2 , Z 3 - We 
now introduce three unit vectors k 1; k 2 and k 3 with origins at X, 
defined as follows: k x is tangent at X to the parametric line of z ,, 
and points in the direction of Z\ increasing; k 2 and k 3 are defined ana¬ 
logously with respect to the parametric lines of z 2 and z 3 - Figure 54 



shows these parametric lines and associated unit vectors. In the 
case of cylindrical coordinates, the associated unit vectors kj, k 2 and 
k 3 were introduced in § 30. 

In the remainder of the present section we shall consider only the 
case when the three vectors k x , k 2 and k 3 are mutually perpendicular. 
The curvilinear coordinates are then said to be orthogonal. We shall 
also require that k 1; k 2 and k 3 form a right-handed triad, the term 
„right-handed” being as defined in § 6. 

The distance ds between two adjacent points is given by the relation 

(49.4) (*) 2 = (dx^+(dx 2 y+(dx 3 ) 2 = 2 {dx r y, 

r= 1 

where dx 1 , dx 2 and dx 3 are the infinitesimal differences between the 
rectangular cartesian coordinates of the two points. From Equations 
(49.3) we have 
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Substitution from these relations in (49.4) then yields for (ds) 2 a 
homogeneous quadratic expression in the quantities dz x , dz 3 and 
dz 3 . Because the coordinate system z x , Z 2 , Z 3 is orthogonal, it turns 
out that no terms involving products of these differentials appear in 
this quadratic expression, and we can then write 

(49.5) {ds) 2 = {h x d Zl )*+ [h 2 dz 2 )*+{h 3 dz 3 )\ 

where h x , h 2 and h 3 are known positive functions of Z\, Z 2 and z 3 - 
The right side of (49.5) is called the fundamental quadratic form, or the 
metric form. 

Les j ls s 2 and j 3 denote the arc lengths of the three parametric 
lines in Figure 54. From (49.5) we then have 

(49.6) ds x — h x dzi, ds 2 = h 2 dz 2 , ds 3 =h 3 dz 3 - 

Let us consider z 3 - It is a scalar field. We shall now consider V^i- 
Because of Theorems 2 and 3 of § 44 we see that 

(i) V*i, being normal to the surface Z\ = constant, has the same 
direction as k x , 

(ii) |V^i| is equal to the maximum value of the directional deri- 
tive dzfds, 

(iii) this maximum value arises when the directional derivative is 
taken in the direction of k x , and is hence equal to dzfds 1 . 

Because of (ii), (iii) and Equation (49.6) it follows that |V2i| = 1/Ai, 
and because of (i) we then have SJzi = ki/Aj. Similar observations 
regarding z 2 and z 3 then permit us to write 

(49.7) k x = h x V^x, k 2 = h 2 \Jz 3 , k 3 = h 3 \/z 3 . 

From these equations we get the relation 

^•(kjaXks) = h x h 2 h 3 V Z%- (V4 2 'X V^ 3 ) • 

But the left side of this equation is equal to one since k 1( k 2 and k 3 
form a right-handed orthogonal triad of unit vectors. Thus 
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( 49 . 8 ) 


V^(V^XV^,)=^- 


50. The expressions XJf V-b and Vxb in curvilinear coordinates. 
We consider first the expression V/- Now/is a function of z 1} Z 2 and 
Z 2 - Hence by Equantion (45.5) we have 

Because of Equation (49.7) we may then write this relation in the 


form 



(50.1) 


1 

or 



(50.2) 

3 

V/= E 

r=. 

UL 

i h, dz. 


.1 3f_. 
h 2 dz 2 ' 


h 3 dzi 


We now turn to V-b. We first express b as a linear function of the 
unit vectors k x , k 2 and k 3 , in the form 


(50.3) 


b = b x ki+^2 k 2 +6 3 k 3 . 


It should be noted that b 1 , b 2 and b 3 are here the orthogonal projec¬ 
tions of b on the lines of action of the unit vectors k 1; k 2 and k 3 , 
respectively. In order to deduce conveniently the desired expression 
for V-b, we now make a rather unusual step. Since 

k 3 — k 2 X k 3 , k 2 — k 3 X k 3 , k 3 — k 3 x k 2 , 
then substitution from Equations (49.7) yields 
k 3 = A 2 A 3 V^ 2 x V£ 3j 

k 2 = h 3 h 3 V-^3 X V^i, 
k 3 = hji 2 V Z\ X V Z 2 • 

We now substitute from these expressions in Equation (50.3), and 
operate on the resultant expression with the operator V • to obtain 
the relation 
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(50.4) V-b = V- {bxh 2 h z X7z 2 X Vz 3 ) + V • (6 2 Mi V-t 3 X V^x) 

4“ V • (A 3 AxA 2 V* x X V* 2 ) • 

For the first of the three terms on the right side of (50.4) we can write, 
because of relation (48.1), 

(50.5) V • [b x h 2 h 3 V-j 2 X V£ 3 ) = V(AxA 2 A 3 ) • (V£ 2 X V£ 3 ) 

+ b 1 h 2 h 3 V-(V-C 2 XV z 3 ) • 

But because of Equation (45.5) we have 


V (W 3 ) — ^ (b\b 2 b z ) V*x + ( b 3 h 2 h 3 ) X7z 3 + ( b x h 2 h 3 ) V.j 3 . 

Thus, since V*i, V£ 2 and V£ 3 are mutually perpendicular, we obtain 

V (b 1 h 2 k s )-(Vz 2 x S7z s ) = ^r(MA) Vq'(Vt 2 xV^: 3 ) 
or, by Equation (49.8), 

(50.6) ^A b A) ■ {Vz 2 xVz 3 ) = (MA) ■ 

Also, because of Equation (48.3) we have 

(50.7) V • (V£ 2 x Vz 3 ) = S7z 3 - (V x Vz 2 ) - Vz 2 - (V x X/z 3 ) = 0 

because of Equation (48.6). The Equations (50.6) and (50.7) now 
permit us to write (50.5) in the form 

1 d 

V-AkA Vt 2 XVt 3 ) = i b Ah) ■ 


This relation and two similar relations involving the second and third 
terms on the right side of Equation (50.4) then permit us to write 
(50.4) in the form 


(50.8) 


V-b = 


Ah 2 h 3 ) + ~ (b 2 hA) + (MA) 


We note that if the curvilinear coordinates happen to be rectangular 
cartesian coordinates, then = h 2 = h 3 = \ and (50.8) reduces to 
(46.9), as expected. 

Finally, we turn to the expression Vxb. Because of Equations 
(49.7) and (50.3) we have 
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or 


Vxb = Vx(6 1 k 1 +i 2 k2+£ 3 k 3 ), 


(50.9) V Xb = V X (Mi V*i) + V X (bJi 2 Vz 2 ) 

+ VX {b 2 h 3 Vz 3 ) • 

For the first term on the right side we can then write 

(50.10) v x (Mi = v (M0 x vci+Mi (v x v^) , 


because of Equation (48.2). But VxVti = 0 by Equation (48.6), 
and by Equation (45.5) we have 


v (Mi)x v*]. = 


^(Mi)v^ + ^- (Mi) v* 2 


+ ^ (Mi) Vz 3 


x V^. 


Now Equations (49.7) yield 


Vz x x Vz x = 0, 

V^ 2 x Vzi = ~k 2 xk 1 =~^k 3 , 
Vz 3 XVz t =^-k 3 xk 1 =^k 2 . 


Thus Equation (50.10) reduces to 


V x (^i^i V^i) 


k 2 <7 

b 3 bi dz 3 


/ L L \ ^3 3 

1 ^ Mi dz 2 


(Mi) • 


This relation and two similar relations involving the second and third 
terms on the right side of Equation (50.9) permit us to write Equation 
(50.9) in the form 


(50.11) 


V x b = 


bi 

h 2 h 3 


_L -Mi 

h 3 h-i 


+ 


h\h 3 


JL A 

dz 2 - dz 3 

c) c) 

d ^ 3 (Mi) - dz x ( Ms) 

' d . d 

dz x (M2) ~dz 3 (Mi) • 
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We note that if the curvilinear coordinates happen to be rectangular 
cartesian coordinates, then (50.11) reduces to (46.12), as expected. 


Problems 

1. If f = {x 1 ) 2 x 2 -\-(x 2 ) 2 x 3 - x 1 x 2 x 3 , find the directional derivative 
of f at the point 4(1,-4,8) in the direction of the position-vector 
a of A. 

2. If / = x 1 sin (r.x 2 ) +x 3 tan (r.x 2 ), find the directional derivative 
of f at the point A (1, 0, -2) in the direction of the vector drawn from 
A to the point B (3, -3,4). 

3. Find a unit vector normal to the surface x i x 3 - x 3 x 1 +x 1 x 2 -1=0 
at the point 4 ( 1, 2, -1). 

4. Two surfaces x x x 2 - (x 3 ) 2 +15 = 0 and (x 2 ) 2 - 3x 3 +5 = 0 inter¬ 
sect in a curve C. At the point A (3,-2,3) on C find (i) the angle 
between the normals to the two surfaces, (ii) a unit vector tangent 
toC. 

5. If /= (x 1 ) 2 +2{x 2 ) 2 +2(x 3 ) 2 , find the maximum value of the 
directional derivative of f at the point (1,-2, -4). 

6. If x = prove that Vx = x/x, and that 

V-x" = rcx"' 2 x, 

where re is a constant. 

7. Find SJr and V0, where r and 0 are the usual plane polar coor¬ 
dinates. Also, find the magnitudes and directions of \/r and V0- 

8. If f = r 3 - cos 2 0, where r and 0 are plane polar coordinates, 
find V / in terms of r, 0 and the unit vectors i x and i 2 associated with 
the corresponding rectangular cartesian coordinates. 

9. If f and g are scalar fields, prove that 

V(//*) =£' 2 (.g Vf-fVg). 

10. If f = (x 1 ) 2 +x 3 V(x 1 ) 2 -|-(x 2 ) 2 andg = x 1 x 2 x 3 , find at the point 
.4(3, 4, 5) the expressions V {fg) and V (f/g ). Note Theorem 4 of 
§ 45, and Problem 9 above. 

11. If f = x 1 x 2 x 3 , a = Xjii - x 2 i 2 and b = x 3 X]i 2 - , compute 

the following: (i) (a-V)/, (ii) (a-V)b, (iii) (aXV)/, (iv) (aXV)-b, 
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(v) (aXV) xb, (vi)V-b, (vii) VXb, (viii) a-(V xb). 

12. Let S and S' be two rectangular cartesian coordinate systems. 
The axes of S can be moved into coincidence with the axes of S' by 
a positive rotation of Jrc radians about the x 3 axis, followed by a 
second positive rotation of At r about the bisector of the angle between 
the positive axes of x 1 and x 2 . (i) Express the coordinates of S' in terms 
of the coordinates of S, and conversely, (ii) If / = x 2 x 3 -\-x 3 x 1 and 
b = (A^+Xg) i i +(x 1 - x 2 ) i 2 +x 3 i 3) express / and b in terms of quan¬ 
tities pertaining to the system S'. 

13. Prove Equation (47.3). 

14. Starting from Equations (46.8) and (46.11), verify Equatiors 
(48.2), (48.4), (48.6), (48.7), (48.8), (48.9), (48.10) and (48.11). 

15. Using the identity (48.14), verify Equations (48.4) and (48.5). 

16. Using the permutation theorem for scalar triple products, 
verify Equation (48.7). 

17. If a is a constant vector, prove that V (a-x) = a. 

18. Prove that ax (V X x) =0. 

19. If a is a constant vector, prove that V X (a X x) = 2a. 

20. Prove that 

V (a• b) =a(V-b)+b(V-a) + (axV)xb+(bx V)xa. 

21. Prove that 

ax(VXb) - (ax V) Xb = a(V -b) - (a- V)b. 

22. Prove that 

(a-V)a = A Va 2 - ax (V Xa). 

23. Prove that 

a(V-a) = l Va 2 -(ax V) Xa. 

24. If f(x l3 x 2 ,x 3 ) is a homogeneous polynomial of degree n, prove 
that (x-V)/= nf. 

25. Prove that if a is a constant vector, then 

(aX V)'(bxc) =(a-b) (V-c) + (c- V) (a-b) 

- (a-c) (V-b) - (b- V) (a-c). 

26. If r, 6, z are cylindrical coordinates, describe their parametric 
surfaces and show that 

(ds)* ■= [dr) 2 +r 2 (dQ)*+(dz) 2 . 
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27. If r, 0, 9 are spherical polar coordinates, describe the paramet¬ 
ric surfaces and lines, and show that 

(■ ds ) 2 = (/r) 2 -|-r 2 (/ 0 ) 2 4 -r 2 sin 2 0 ((ftp) 2 . 

28. Prove that for the transformation from rectangular cartesian 
coordinates to orthogonal curvilinear coordinates z x , Z 3 , Z 3 for which 
the metric form is as given in Equation (49.5), the Jacobian / satisfies 
the relation h^h t h z I — 1 . 

29. Express Equation (50.11) in terms of a determinant. 

30. Write out the expressions V/, V-band VXb in the case of 
cylindrical coordinates r, 0 , z • 

31. Write out the expressions V/, V-b and Vxb in the case of 
spherical polar coordinates. 

32. By setting b = V/ in Equation (50.8), deduce an expression 
for V 2 / in terms of general orthogonal curvilinear coordinates. 

33. Show that, in cylindrical coordinates r, 0, z, we have 

V 2 f=^ + I^+^ + I^. 

J dr 2 ' r 2 , 90 2 ^ dz 2 ^ r dr 


34. Solve the differential equation V 2 /' = 0 when / is a function 
only of the cylindrical coordinate r. 

35. Show that, in spherical polar coordinates r, 0, 9 , we have 

<? 2 / 1 <? 2 / 1 2df cot0# 

^ dr 2 r 2 <90 2 ' r 2 sin 2 0 dtp 2 r dr r 2 <90 

36. Solve the differential equation S7 2 f = 0 when / is a function 
only of the spherical polar coordinate r. 

37. If k 1; k 2 and k 3 are the unit vectors associated with the spheri¬ 
cal polar coordinates r, 0 and 9 , prove that 


_ 0 - k 

dr -U) dQ ~ k2> 


^ = 0 
dr ’ 


<9k 2 

dd 


= k i 


dk t 

d(f> 

dk 


= k 3 sin l 


dks 

dr 


= 0 , 


^? = 0 
<90 


^ = k 3 COS0, 

= -kj sin '0 -k 2 cos 0 . 


Note: the corresponding problem in the case of plane polar coordi¬ 
nates is worked out in Chapter III, § 30. 
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Chapter V 

INTEGRATION 


51. Line integrals. A curve is called a regular arc it can be represented 
in some rectangular cartesian coordinate system by the equation 

x =* 1 (a)i 1 +* 2 (a)i 2 -h* 3 (a)i 3 , 

where x is the position-vector of a general point X on the curve, a 
is a parameter with the range a< a <(3, and for this range of a the 
functions * 1 (a), x 2 (u) and * 3 (a) are continuous with continuous first 
derivatives. A curve which consists of a finite number of regular arcs 
joined end to end, and which does not intersect itself, is called a 
regular curve. Throughout this chapter we shall consider only regular 
curves, and shall refer to them simply as curves. 

Let us consider a curve C with terminal points A and B, as shown in 
Figure 55. Let f {x 1 ,x 2 ,x 3 ) be a function which is single valued and 



Figure 55 


continuous on the curve C. We divide C into jV parts by the JV+1 
points Q_ 0 , Q_ x , Q_ 2 , • • •, Q,^, as shown. The length of the line seg¬ 
ment Q.^iQ.>) {p= 1,2, • • -, N) is denoted'by A s p . Let X p beapoint 
on the arc Q W Q^, and let its coordinates be (x pl , x p2 , x p3 ). The line 
integral of f over C is then defined to be 

N 

(51.1) lim E f(x pl ,x p2 ,x p3 ) &s p = [fds. 

jV->- oo P= 1 
As p + 0 

This limit is independent of the manner in which the curve C is divided 
into parts, since / is continuous and single valued on C. If /= 1 
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everywhere on C, then Equation (51.1) defines the arc length of C. 
Let X be a general point on the curve C, as shown in Figure 56. 



Let j denote the arc length of C measured from the end A of C. The 
vector dx/ds was seen in § 28 to be a unit vector tangent to C in the 
direction of s increasing. Denoting this vector by t, we have 


(51.2) 



Let b(* ls x 2 , x 3 ) be a vector field defined over C. The orthogonal 
projection of b on the unit tangent vector t is called the tangential 
component of b. If we denote it by b t , we have 

(51.3) b t = b t. 

The line integral of b over G is defined to be 

(51.4) Jb, ds = fh-t ds. 

c c 

Because of (51.2), we have 

(51.5) t ds = dx = 


Thus 

(51.6) Jb t ds = Jh-dx = J{b 1 dx 1 -\-b 2 dx 2 -\-b 3 dx 3 ). 

c c c 

We may also consider line integrals over C with integrands which 
are vectors. The integration of vectors was defined in § 13. Following 
this definition, we have 


(51.7) 


Jh ds 
c 


i x J b x ds+i 2 J b 2 ds+i 3 J b 3 ds, 

c c c 
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(51.8) fhxtds = fhxdx 

c c 

= i if {b 2 dx 3 -b 3 dx 2 ) + i 2 J {b 3 dx 1 -b 1 dx 2 )+i 3 J (b^-b^Xj) . 
c c c 

To apply the above considerations to two-dimensional problem)! 
involving line integrals along curves in the x 1 x 2 plane, it is only neces¬ 
sary to set b 3 = x 3 = 0 in the above formulas. A few examples will 
now be worked out. 

Example 1. Let f = (*i) 2 -Hx 2 ) 3 , an d let us evalute the line integral 
of f along the straight line x 2 = 2x x in the x 1 x 2 plane from the origin 
to the point B (2,4). This problem is two-dimensional. There are 
several ways of solving this problem, since a curve may be represented 
parametrically in many different ways. We present here two solutions. 

(i) The curve C can be represented by the equations 

x 1 = u, x 2 = 2 u, 0<a<2. 

Thus 

dx = i 1 dx 1 -\-i 2 dx 2 = (ij^ —1- 2i 2 ) du, 
ds = |rfx| = du. 


Also, / = u 2 +8u 3 , whence 
2 

Jf ds = J (u 2 - f-82/ 3 ) V5 du 
c o 


104VS 

3 


(ii) On C we have 

*2=2*!, /= (*l) 2 + 8(*!) 3 , 

Thus 

2 

Jfds = /[(*!) 2 +8(*!) 3 ] V 5 dx ! = 

C 0 

Example 2. Let b = A: 2 i x A: 1 ) 2 i 2 -f-ATii 3 , and let us evaluate 
the line integral of b over the curve C in Example 1 above. We 
present two solutions. 
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(i) We have 

(51.9) b-ifx = b 1 dx 1 -\-b 2 dx 2 -\-b 3 dx 3 

= * 2 <fei+(*s+*i) 2 dx 2 +x 1 dx 3 _ 
_ L I I .. u dx 2 , .. dx 3 




But on C we have x 2 = 2a: 1j * 3 = 0, so 


fh-dx = j [ 2 z 1 + 2 (z 1 ) 2 ] dx j = : 


(ii) From (51.9) we have 

b = *2 ^ + (*3+*l) 2 +*l ^ ^2 

But on C we have ^ x 3 — 0, so 

4 

fh-dx =f [|* 2 +i(x 2 ) 2 ] dx 2 = ^-. 

C 0 


Example 3. Let b = a 2 x 1 i 1 -\-ax 2 x 3 i 2 -\-x 1 (x 3 ) 2 i 3 where a is a constant, 
and let us evalute the line integral of b along the curve C in Figure 57 



Figure 57 

from the point A to the point B, the curve C being a portion of the 
intersection of the cylinder ( a : 1 ) 2 + ( a : 2 ) 2 = a 2 and the plane x 1 +x 3 = a. 
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Now 


h-dx — 


bl+b *^ + b3 S' dx ' 

a?x x +ax 2 x 3 ~ + *iW 2 


dx. 


But on C we have 


(* 2 ) 2 = a 2 - Oj) 2 , x 3 = a - x x (a>x x >0 ), 
dxz _ _*i <& 3 _ t 

rfx x x 2 <A x 


Thus 

o 

J b ■ dx = J [a 2 x x - a (a - x 1 )x 1 - x x (a- Xj) 2 ] dx x = - i a 4 . 
c a 


52. Surface integrals. A regular surface element is defined to be a 
portion of a surface which, for some orientation of the coordinate 
axes can be projected onto a region S' in the x x x 2 plane enclosed by a 
regular closed curve, and which can be represented by the equation 
x 3 = g(x x , x 2 ), where g(x 1; x 2 ) is continuous and has continuous first 
derivatives in S'. In this chapter we shall consider only surfaces com¬ 
posed of a finite number of regular surface elements. 

Let us consider a surface S bounded by a closed curve C, as shown 
in Figure 58. Let/(x 1 , x 2 , x 3 ) be a function which is single valued and 



continuous on S. We divide the region 5 into JV parts with areas 
A S p (p = 1,2, • • •, JV). Let X p be a point on the element of area A S p , 
as shown, and let its coordinates be (x #1 , x fi , x p3 ). The surface 
integral of f over S is defined to be 


134 



(52.1) 


lim S (xp 1} Xp^, Xp S ) A Sp — /' J dS . 

JV-»- OO />= 1 

AS,-0 

This limit is independent of the manner in which S is divided into 
parts, since / is continuous and single valued throughout the region 
S. If / = 1, then Equation (52.1) yields the surface area of S. 

Let us suppose that S is a regular surface element, as shown in Figure 
59. To evaluate the surface integral in Equation (52.1), we let X be 



a general point on S and let dS- be an element of area at X. We now 
project S onto the region S' in the x x x 2 plane, X and dS projecting into 
X' and dS', respectively. Let n be the unit vector normal at X to the 
surface S, making an acute angle with the x 3 axis. Then 

dS' = n 3 dS. 


Let the equation of the surface S be 
(52.2) * 3 -«(*!, * 2 ) = 0- 


Denoting the left side of this equation by G, we have by Theorems 
and 3 of § 44, 


„|VC[-VG-!|i 1 +gi,+ 


dX' 


3xo 


_ 3 1 ; 1 • 

3 ~ 1 dx, 2+ 


dx x 


dG . 

13 


2 
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Thus 


2 


(52.3) 


We then have 
(52.4) JfdS 

s 


n 3 




-i 



S' 




i 

dS‘. 


To evaluate the integral on the right side, we may write dS' = dx x 
dx 2 , and then perform a double integration with respect to x x and 
x 2 over the region S'. Or we may use polar coordinates r and 0 in 
the x x x 2 plane, writing dS' = r dr dQ. 

If S is not a regular surface element, we divide it into regular sur¬ 
face elements. The surface integral of/ over S is then found as the 
sum of the surface integrals of/ over these regular surface elements. 

To distinguish between the two sides of a surface S , let us designate 
one side as the positive side and the other as the negative side. Let 
n be the unit vector normal to S at a general point X, and lying on the 
positive side of S. Let b(* 1; x 2 , x 3 ) be a vector field defined over S. 
The orthogonal projection of b on n is called the normal component 
ofb. If we denote it by b„, we have 

(52.5) b n = bn. 

The surface integral ofb over S is defined to be 

(52.6) f b n dS = fh-ndS. 

s s 

It is often convenient to introduce an infinitesimal vector dS defined 
by the relation 

(52.7) dS =n dS, 


so Equation (52.6) may take the form 

(52.8) fb n dS=fb-dS. 

s s 

It is sometimes necessary to consider surface integrals with inte¬ 
grands which are vectors. For example, we have, following the de¬ 
finition of integration of vectors in § 13, 
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(52.9) fb dS = ii f b 1 dS+i 2 f b 2 dS+ i 3 Jb 3 dS, 

s s s s 

(52.10) JhXndS = *if (b 2 n 3 -b 3 n 2 )dS-\-i 2 J (b^—b^n^dS 

s s s 

+i 3 J {b x n 2 -b 2 n^dS. 
s 

We shall now work out some examples. 

Example 1. Let / = (x 1 ) 2 + 2 x 2 +x 3 -l, and let us evaluate the 
surface integral of/ over a region S consisting of that part of the plane 
2x 1 -)-2a: 2 +x 3 = 2 lying in the first octant. 

The region S is shown in Figure 60. It is a regular plane element. 



Figure 60 


We have < 2*3 = _ 2 = - 2 

dx l ’ dx 2 

whence (52.3) yields n 3 = 1/3. Thus dS = 3 dx 1 dx 2 , and so 
J fdS = J 3/ dx x dx 2 . 

s s 

But on S, we have 

* 3 = 2 - 2x x - 2x 2 , 

f = (*i) 2 +2* 2 +(2 — 2x x - 2 * 2 ) ~ f 
= (^i-l) 2 - 
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Thus, if S' denotes the projection of S on the x x x 2 plane, we have 

1 1-X t 

JfdS =3 if (%- \) 2 dx x dx 2 =3 <J J (x x - 1 ) 2 dx 2 dx x = 

S S' 0 0 

Example 2. If b = x 2 i 1 +x 3 i 2 , evaluate the surface integral of b 
over the region S in Example 1, the origin being on the negative side 
of S. 

We have 

71 — — 71 — — Tin — ! 

— 3 J ' L 2 3 5 "3 3 3 

b n = b x n x +b 2 n 2 +b 3 n 3 = f (x 2 +x 3 ), 
dS = 3 dx x dx 2 . 

Thus 

Jb n dS = 2 J{x 2 +x z )dx 1 dx 2 . 
s s 

But on S we have x 3 = 2 - 2x x - 2x 2 , whence 

1 1-Xj 

Jb n dS = 2 J(2-2x 1 -x 2 )dx 1 dx 2 = 2 J f (2-2x 1 -x 2 )dx 1 dx 2 = l. 


53. Tn/>/« integrals. Let V be a region in space enclosed by a sur¬ 
face S as shown in Figure 61. Let f (x 1} x 2 , x 3 ) be a function which is 



Figure 61 

single valued and continuous throughout V. We divide V into JV 
parts with volumes AF # {p = 1, 2, • • •, JV). Let be a point in the 
element of volume AV P , as shown, and let its coordinates be (x pl , 
x t2 , x p3 ). The triple integral of f over V is defined to be 

(53.1) lim 2 f{x n ,x p2 ,x p3 ) AV p = f fdV. 

JV-». CXD /’“l y 
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This limit is independent of the manner in which V is divided into 
parts, since/is single valued and continuous in V. 

We may also consider triple integrals with integrands which are 
vectors. Thus, if b is a vector field which is single valued and contin¬ 
uous in V, we have, following the definition of integration of vectors 
in §13, 


(53.2) 


/WF 


— *i J bidV-j~i 2 I b 2 dV -f -13 J~ b 3 dV. 

v v v 


To evalute triple integrals, one may divide the region V into ele¬ 
ments by means of three systems of planes parallel to the coordinate 
planes. The value of the integral is then found by the performance of 
three integrations with respect to the rectangular cartesian coordinates 
x-l , x 2 and x 3 . Or we may divide V using parametric surfaces of a 
curvilinear coordinate system, in which case we evaluate the triple 
integral by performing three integrations with respect to the three 
curvilinear coordinates. Since most readers will have had considerable 
experience in evaluating triple integrals in connection with elementary 
calculus, no example need be given here. 


Problems 

1. If x is the prosition-vector of a general point on a circle C of 

radius a, and t is the unit tangent vector to C, evaluate J t-dx. 

c 

2. If f = (xj) 2 - (x 2 ) 2 , evaluate the line integral of/ along the line 
x 1 J r 2x 2 = 2 from the point A (0,1) to the point 5(2,0). 

3. Evaluate the line integral in Problem 2 when the curve C consists 
of (i) the two line segments AD and DB. where D has coordinates 
(1,1), (ii) the two line segments AO and OB, where 0 is the origin. 

4. If / = 81*!- 9, evaluate the line integral of / along the curve 
( x 2 ) 2 = (x-l) 3 from the origin to the point (1,1). 

5. If / = x 2 x 3 /XjXj+XjXj, evaluate the line'integral of/from the 
origin 0 to the point 5(1,2,3) along the path consisting of (i) the line 
segment 05; (ii) the three line segments 0D, DE and EB, where D 
and E have coordinates (1,0,0) and (1,2,0), respectively. 
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6 . If/ = x 1 -\-x 2 -{-x 3 , evaluate theline integral of/ along the curve 

x = a cos u i^a sin u i 2 +a u cot a i 3) (0<u<i7i), 
where a and « are constants. 

7. If b = (*i - ^ 2 )i 1 +Ar 2 i 2 , evaluate the line integral of b along the 
following curves: (i) the path in Problem 2 above, (ii) the two 
paths in Problem 3 above, (iii) the path in Problem 4 above. 

8 . If b = 2x 1 x 2 i 1 +[(x 1 ) 2 - (x 2 ) 2 ]i 2 , evaluate the line integral of b 
from the point .4(0,0) to the point 5(1,1) along the following curves: 
(i) * 2 = x l5 (ii) (x^ 2 = x 2 , (iii) x x = (x 2 ) 2 . 

9. If b is as given in Problem 8, evaluate the line integral of b from 
the point A (a, 0) to the point 5(0, a) along the circle (x 1 ) 2 +(x 2 ) 2 = a 2 , 
using the polar angle 0 as the parameter varying along the curve. 

10. Evaluate the line integral in Example 3 of § 51, using for the 
curve C the parametric representation 

x = a cos u i x +a sin u i 2 +a(l-cos u) i 3 , (0<«<in). 

11. If b = - x x i 3 , evaluate the line integral of b along 

the curve in Problem 6. 

12. If C is the curve in Problem 6, evaluate the following: 

(i) Jxds, (ii) JxXtds. 
c c 

13. If / = x 1 -)-x 3 , evaluate the surface integral of f over the region 
S consisting of the triangle cut from the plane 6XJ+3X2-I-2X3 = 6 by 
the three planes x x =0, x 2 = 0, x 1 -\-x 2 = 1. 

14. If b = x 1 i 1 +(x 2 ) 2 i 3 , evaluate the surface integral of b over the 
region S in Problem 13, the origin being on the negative side of S. 

15. If/= (xj) 2 +(x 2 ) 2 , evaluate its surface integral over the region 
S consisting of the part of the surface x 3 = 2 - (x x ) 2 - (x 2 ) 2 in the first 
octant. Use polar coordinates in the x ± x 2 plane. 

16. Find the area of the region S in Problem 15. 

17. If b = (x 1 ) 2 i 1 +x 3 i 3 , evaluate the surface integral ofb over the 
region S of Problem 15, the origin being on the positive side of S. 

54. Green’s theorem in the plane. This theorem is as follows. Let S 
be a closed region in the x x x 2 plane bounded by a curve C, as shown in 
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Figure 62 


Figure 62, and let b be a vector field continuous and with continuous 
first derivatives in the region S. Then 

(54.1) yi 3 -(Vxb) dS = Jh-t ds, 

s c 

where the integration over C is carried out in the positive direction, 
that is, in that direction of travel around C in which the interior of S 
lies on the left. This direction is indicated by the arrow on C in Figure 
62. In Equation (54.1), i 3 is as usual a unit vector which forms with 
ij and i 2 a right handed triad, and t is a unit vector tangent to C in 
the positive direction. 

Equation (54.1) is often written in the form 




+bJx , 


We first prove that 


= -f V*i> 


in the case when C can be cut by a line parallel to the x 2 axis in two 
points at most. Figure 63 illustrates the situation. On C there are 
two points D and E where the tangent to C is parallel to the x 2 axis. 
Let d and e be the abscissas of D and E, respectively. These points 
divide C into two parts C and C". At a general point X(x 1 , x 2 ) in 5 
we introduce an element of area lying in a strip parallel to the x 2 axis, 
the left edge of the strip cutting C and C" at the points X' (x 1; x 2 ) 
and X"(x 1 , x 2 '), as shown. Then 
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/ 


dbi 

dx 2 



dx 2 dx x 




e < 

= / M*i> x 2 ') dx x -f b x (x x , x 2 ") dx x 

d d 

d t 

= -Jh(x 1 , X 2 ') dx x -f b x (x x , x 2 ") dx x 

e d 

— - f b x dx x . 
c 


Let us now consider the case when C can be cut by a line parallel to 
the x 2 axis in more than two points, such as the case of the curve C 
in Figure 62. Here we have only to join the points F and G where 
there are tangents parallel to the x 2 axis by a curve K which is contained 



in S and which cannot be cut by a line parallel to the x 2 axis in more 
than one point. The curve K divides S into two portions for both of 
which the above proof of (54.3) holds. Hence, is we apply (54.3) to 
both portions, and add, the two line integrals over K cancel, and we 
hence establish (54.3) for the entire region S. In a similar fashion we 
can prove Equation (54.3) for the case when several curves such as 
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K are required. And we can do likewise when S is multiply connected, 
that is, S has holes in it and C then consists of several isolated parts. 

In an analogous fashion we can prove that 

(54.4) fgtS-f 

s c 

Subtraction of (54.4) from (54.3) then yields (54.2). This completes 
the proof. 


55. Green's theorem in space. This theorem is as follows. Let V be 
a closed region bounded by a surface S, and let b be a vector field 
continuous and with continuous first derivatives in V. Then 


(55.1) 


fV-b dV = Jb-n dS , 


where n is the unit outer normal vector to S. 

This theorem is also called the divergence theorem, and can be 
written in the form 


(55.2) 


/(!I 


I ^2 , 

i r\.. i 


3x 0 


dx- 


dV = J (Mi+Vs+Ms) d S - 


We shall first prove that 

(55.3) f dV = J b 3 n 3 dS, 

v * 3 s 

in the case when S can be cut by a line parallel to the x 3 axis in two 
points at most. Figure 64 illustrates the situation, T being the projec¬ 
tion of S on the x x x 2 plane. On S there is a curve C consisting of points 
where the tangent plane to S is parallel to the x 3 axis. The curve C 
cuts S into two portions S' and S". At a point X(x ly x 2 , x 3 ) in V we 
introduce an element of volume lying in a prism parallel to the x 3 
axis, the vertical line through X meeting S' and S" at the points X' 
(*!, x 2 ,x 3 ) and X"(x 1} x 2 , x 3 ), as shown. Thus 
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Let n' be the unit outer normal vector at X', and let dS ’ be the area of 
of the element cut from S' by the vertical prism. Let us define n" and 
iS" at X" similarly. Then we have 

dx 2 dx | — ti 3 dS — j dS . 

and so we may write (55.4) in the form 

fllx dV= I ' ds '+f b s( x i> x 2 > x 3") n 3"dS"=Jb 3 n 3 dS. 

V s S' S" s 

Let us now consider the case when S can be cut by a vertical line in 
more than two points. In such cases we can always divide V into a 
number of regions V 1} F 2) ■ • • by cutting V by a number of surfaces 
so chosen that the boundary of each of the regions V l3 
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V 2 , ■ ■ ■ can be cut by a vertical line in two points at most. The above 
proof of (55.3) then applies to the regions V 1} V g , ■ • •. If we apply 
(55.3) to V 1; V 2 , ■ ■ ■ and add, the surface integrals over K u K 2 , 
cancel, and we hence establish (55.3) for the entire region V. 

In a manner similar to the above we can prove that 

(55 ' 5 >' f&r=!>*<*• 

V S V d S 

Addition of Equations- (55.3) and (55.5) then .yields (55.2). This 
completes the proof. 

If/ is a scalar field with continuous second order derivatives, then 
we can set b = V / and substitute in Equation (55.1) to obtain 

/ V- (V/) dV = f SJf-ndS, 

V s 

or 

(55.8) f(V-V)fdV=f d £dS, 

v s 

where V • V is the familiar Laplacian operator, often denoted by V 2 , 
and dfjdn is the directional derivative of/in the direction of the outer 
normal to the surface S. 

56. The symmetric form of green's theorem. Let / and g be scalar fields 
with continuous second derivatives in a closed region V bounded by 
a surface S. We may then apply Green’s theorem as stated in Equa¬ 
tion (55.1), but with the vector b replaced by / Vg. This yields 

(56.1) f V • {fVg) dV = f f'tfg-n dS. 

V s 

But 

V-(/V*) =/(7-V)*+V/-Vf 
=/V 2 £+V/7£, 

because of Equation (48.1). Also, Vg-n is equal to the directional 
derivative dgjdn of g in the direction of the outer normal n to S. Thus 
Equation (56.1) becomes 
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(56.2) 


/ (/V«*+V/-V*)rfF= f/f n dS. 

V $ 

Similarly, by an interchange of/and g in the above, we obtain 

(56.3) f(g V»/+V*-V/) dV=fgf n dS. 

V s 
Subtraction of (56.3) from (56.2) then yields 

(56.4) /(/- S V</) ir~j{f‘£-S 

V s 

This equation is called the symmetric form of Green's theorem. 

57. Stokes's theorem. This theorem is as follows. Let S be a closed 
region on a surface, the boundary of S being a curve C. We choose 
a positive side for S, and let n be the unit vector normal to S on the 
positive side. The positive direction on C is defined to be that in which 
an observer on the positive side of S would travel to have the interior 
of S on his left. Let t be the unit vector tangent to C in the positive 
direction, and let b be a vector field with continuous first derivatives 
in the closed region S. Then Stokes’s theorem states that 

(57.1) J n-(VXb) dS — Jh-tds, 

s c 

where the integration around C is carried out in the positive direction. 
This theorem can also be written in the form 



= y\v*i+M* 2 +M* 3 ) • 

c 


We shall first prove that 

(57.3) / n- (Vx6 1 i 1 ) dS = Jb 1 dx 1 , 

s c 

in the case when S is a regular surface element and the positive side of 
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S is that side on which the unit normal vector n points in the direction 
of increasing x 3 . Figure 65 illustrates the situation, and shows the 



unit tangent vector t of C and also the region S' in the x x x 2 plane into 
which S projects. 

Now 

(57A) /»'(Vxi 1 i 1 ) dS =/n-(i, ^-i 3 ^jdS. 

Let the equation of the surface S be x 3 = i g(x 1 , x 2 ). Then on S, we have 


bl[ X l> X 2> X s( X l> X 2)] — C l( X l> X 2) J 

dc x _ db x db x dx 3 

dx 2 dx 2 dx 3 dx 2 

We now substitute from this equation for db x ldx 2 in Equation (57.4) 
to obtain the relation 


(57.5) 


(57.6) /»• (V X MO 4S=-fn i, | JS+fn- (i, + i, ^ ^ dS 
s s s 

= - A+4> 

where I x and I 2 denote the two integrals on the right side of this 
equation. 

Let us consider I x . We have n-i 3 dS = n 3 dS = dS\ where dS' is 
the projection of dS on the x x x 2 plane. Since c x is a function of x x and 
x 2 only, we can then write 


147 




By Green’s theorem in the plane, as stated in Equation (54.2), we 
then have 


(57.7) I± — f Cl ( Xl , X2 )d Xl — fb l[XlM M]d Xl —J b^dx^. 
c c c 

We now consider / 2 . The position-vector of the general point X 
on S is 

Hence 


and so 


— *>+ ~ 


dx. 


dx 0 


< 57 - 8 > '•“/•'It'®- 

But the vector <3x/ 3x 2 is tangent at X to the curve of intersection of S 
and a plane parallel to the x 2 x 3 plane. Hence this vector is tangent to 
5 and is then perpendicular to the unit normal vector n, so that 


n 


dx 

ar 2 = 0 ‘ 


Thus Equation (57.8) yields / 2 = 0, and from Equations (57.6) and 
(57.7) we can then conclude that Equation (57.1) is true. 

When the positive side of S is chosen so that the unit normal vector 
n points in the direction of decreasing x 3 , the proof of Equation (57.3) 
is similar to the above, the only differences in the proofs being that 
in the present case n 3 is negative and the direction of integration around 
the curve G is opposite that in the above proof. 

When the surface S is not a regular surface element, we divide it 
into a number of regular surface elements S 1 , S 2 , ■ ■ ■ by a number of 
curves , T 2 , The above proof of (57.3) then applies to the 
regions S 1; S 2 , • ■ If we apply (57.3) to these regions, and add, the 
line integrals over L 1 , Z, 2 , • • • cancel, and we hence establish Equation 
(57.3) for the entire region S. 


148 



In a manner similar to the above we can prove that 

(57.9) / n- (V X b 2 i 2 )dS = j b 2 dx 2 , J n-(V X b 3 i 3 )dS — J b 3 dx a . 

s c s c 

Addition of Equations (57.3) and (57.9) then yields (57.1). This 
completes the proof. 

58. Integration formulas. Green’s theorem in space and Stokes’s 
theorem were considered in §§ 55 and 57, respectively. These theo¬ 
rems are integration formulas which we may write in the form 


(58.1) 

f V-b dV= f n-bdS, 

V s 

(58.2) 

J (n x V) • b dS = Jt-h ds 


s c 


Both of these theorems involve a vector field b; Equation (58.1) 
presents a transformation from a triple integral to a surface integral, 
while (58.2) presents a transformation from a surface integral to a 
line integral. We shall now introduce four other integration formulas, 
which we shall state as two theorems. 

Theorem 1. Let V be a closed region in space bounded by a sur¬ 
face S with the unit outer normal vector n, as in the case of Green’s 
theorem in space. Let/and b be two fields with continuous first deri¬ 
vatives in V. Then 

(58.3) / V /dV = / n fdS, 

V s 

(58.4) f V X b dV = j nxb dS. 

v s 

Theorem 2. Let S be a closed region lying on a surface and bounded 
by a curve C, n being the unit positive normal vector to S and t being 
the unit positive tangent vector to C, as in the case of Stokes’s theorem. 
Let f and b be two fields with continuous first derivatives in S. Then 

(58.5) J(nxV)fdS=ftfds, 

s c 
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(58.6) 


J (nxV)xb(f5’= J txb ds. 

s c 

Proof of Equation (58.3). Let c be a constant vector field. If in 
Equation (58.1) we set b =/c, we obtain 

(58.7) f V • (/c) dV = f n-fc dS. 

v s 

But we have 

V-(/c) = V/-c, 

by Equation (48.1), since c is a constant vector. Thus Equation 

(58.7) may be written in the form 

c fvfdV-fn/dS =0. 
v s 

Since c is an arbitrary constant vector, the expression in the square 
brackets must vanish, whence Equation (58.3) is proved. 

Proof of Equation (58.4). We again introduce the constant vector 
field c, but in Equation (58.1) we replace b by b X c to obtain the re¬ 
lation 


(58.8) f V-(bxc) dV=fn-{hxc) dS. 

v s 

Since c is a constant vector, we have by the permutation theorem for 
scalar triple products the relations 

V-(bxc) =c-(Vxb), n-(bXc) =c-(nxb.) 

Thus Equation (58.8) may be written in the form 

c. J \7 Xh dV - J nxh dS =0. 

v s 

Since c is an arbitrary constant vector, we conclude that (58.4) is 
true. 

Proofs of Equations (58.5) and (58.6). To prove these two equa¬ 
tions we replace b in Equation (58.2) first by/c and then by bXc, 
where c is a constant vector field. The procedure then follows that in 
the previous two proofs. The details are left as exercises for the reader 
(Problem 11 at the end of this chapter). 
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The six integration formulas (58.1)-(58.6) may be written com¬ 
pactly in the form 

(58.9) fv*TdV=f n* T dS, 

v s 

(58.10) J (nxV)*TdS =J t*Tds, 

s c 

where T can denote a scalar field or a vector field, and the asterisk has 
the following meanings: if T is a scalar field, it denotes the multipli¬ 
cation of a vector and a scalar; and if T denotes a vector field, it 
denotes either scalar or vector multiplication. Thus, for example, 
if T denotes a vector field b and the asterisk denotes scalar multipli¬ 
cation, then Equation (58.10) becomes (58.2). 

59. Irrotational vectors. A vector field h(x 1; x 2 , x 3 ) is said to be 
irrotational in a region V in space if everywhere in V we have 

(59.1) V xb = 0. 

Let cp be any scalar field with continuous second derivatives; and 
let us write b = V<p • Then 

Vxb = .V X V<p = 0, 

so a vector b defined as the gradient of a scalar field is irrotational. 

We shall now show that an irrotational vector field b has the fol¬ 
lowing properties: 

(i) Its integral around every reducible circuit in V vanishes. 

(ii) When V is simply connected, b is the gradient of a scalar 
field. 

To verify the first property, we consider a general circuit in V 
which is reducible, that is, it can be contracted to a point without 
leaving V. Let S be a surface entirely in V and bounded by C. If we 
assume that b has continuous first derivatives, then Stokes’s theorem 

(57.1) yields 

J bt ds = J n- (V xb) dS = 0 
c s 

by Equation (59.1). 
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To verify the second property, we let X be a general point in V, 
and let X 0 be a given point. We also let C' and C" be any two paths 
in V from X 0 to X, as shown in Figure 66. Because of property (i) 


x 



above, the line integral of b from X a to X is the same for paths C' and 
C" and hence has the same value for all paths in V from X 0 to X. Thus, 
if we write 

x 

(59.2) cp = J b-dx, 

x. 


then cp depends only on the coordinates {x 3 , x 2 , x 3 ) of X. 
tiation of Equation (59.2) yields dp =b -dx, or 


(59.3) 


dtp _ , dx 

U ~ dT 


Differen- 


But dp/ds is the directional derivative of cp, and by Equation (44.3) is 
equal to V<p- (dx/ds). Thus Equation (59.3) may be written in the 
form 


(V„-b).£-,0. 


Since dx/ds is an arbitrary vector, then 
(59.4) ■* b = V<p. 

This completes the proof. 

The function cp is called a scalar potential function. 


60. Solenoidal vectors. A vector field b(x l} x 2 , x 3 ) is said to be sole- 
noidal in a region V in space if everywhere in V we have 
(60.1) V-b = 0. 
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Let tp be any vector field with continuous second derivatives, and 
let us write b = V X tp . Then 

V-b = V-(VXtp) =0. 

We shall now show that, if b is any solenoidal vector field, there exists 
a vector field tp such that b = V X tp. 

To prove the proposition, we must solve the scalar equations 

(60.2) 

c/x 2 OX 3 

for 9 2 and cp 3 , where b 1} b 2 and b 3 are given functions subject to 
the condition 

db x db 2 db 3 

( 60 - 5 ) ^ + ^2 + ^3 _ 

Let us choose <p x = 0. Then we have from Equations (60.3) and 
(60.4) by partial integrations with respect to x u 


*i 

= f b 3 dx 1 +<\i 2 (x 2 , x 3 ), 

x i 

= -J + s, x s), 


where a x is a constant and <\i 2 and are arbitrary functions of x 2 
and x 3 . To satisfy Equation (60.2) we must have 




_^3 _ 

?x« dx„ 


Because of Equation (60.5) we can then write 




^3 _ ^2 

dx 2 dx 3 


— x 2 , * 3 ) bx{a x , x 2 , x 3 ) 2 X 
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This equation is satisfied if we choose = 0 > 

x a 

= J i J J X 2 , #3) dx 2, 

where a 2 is a constant. The final result is then 

*1 

cpi = 0 , 92 = y* ^s(*i> ^2? *3) ^*15 



93 = -/&*(*!,'*2> *3) ^1+ / *2, *3) *2, 

Oj 

where all integrations are partial integrations, and and <z 2 are con¬ 
stants. The function tp is called a vector potential function. 

In the above proof, several arbitrary selections have been made. 
This indicates that a given solenoidal vector field b does not possess a 
unique vector potential function. In order to see this more clearly, 
we let «p be one vector potential function corresponding to the sole¬ 
noidal vector field b, and let /be any scalar field. Then 

Vx(cp + V/) = VX<p+VXV/ = VXtp = b. 

Thus, cp + V/ is also a vector potential function corresponding to the 
field b. 

If b is any vector field having continuous second derivatives in a 
region F, then b can be expressed as the sum of an irrotational vector 
and a solenoidal vector. The proof of this will not be given here. 

Problems 

1. Let C be a closed curve in the x x x 2 plane. Prove that the area 
A of the region S enclosed by C is given by the relation 

A = 1 J (jq dx 2 - x 2 dx-,) 
c 

where the integration over C is carried out in the direction of travel 
around C in which the interior of S is on the left. 

2. If x is the position-vector of a general point X on a closed surface 
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S, and n is the unit outer normal vector to S, prove that the volume 
V of the region enclosed by S is given by the relation 

V = ±fn-xdS. 
s 

3. If S is a closed surface with a unit outer normal vector n, prove 
that 


JnXxdS = 0. 


4. If V is a region bounded by a surface S, and n is the unit outer 
normal vector to S, prove that 

4 = if [(*i) 2 +(* 2 ) 2 ] (xA+x£t)-n dS, 

s 

where I 3 is the moment of inertia of V about the x 3 axis. 

5. If b has continuous first derivatives in a closed region V bounded 
by a surface S, prove that 

/ n-(VXb) dS = 0. 

6 . If b = a x (x 1 ) 2 i 1 -\-a 2 (x 2 ) 2 i 2 + a 3 (x 3 ) 2 i 3 , where a x , a 2 and a 3 are 
constants, evaluate the surface integral of b over the sphere through 
the origin with center at the point A (a x , a 2 , a 3 ). 

7. If b = (Ar 1 ) 2 i 1 +x 1 x 2 i 2 +x 3 i 3 , evaluate the surface integral of b 
over the cube bounded by the planes x x = 2 , x 2 = 2, x 3 = 2 and the 
coordinate planes. 

8 . If b = [(*0 2 - * 2 ] i 1 +[2(x 1 ) 2 + 3* 2 ] i 2 - 2 x ^3 i 3; evaluate the 
surface integral of b over the sphere 5 with center at the point E 
(1, 0, 2) and passing through the point F( 3, -2, 1). 

9. If C is any closed curve, prove that J dx = 0 . 

10. Let C be the circle with the equations (x 1 ) 2 +(x 2 ) 2 = 4, x 3 = 0, 
and let 

b = [M 2 +x 2 ] ii+[(*i) 2 +* 3 ] i 2 +^2»3- 
Evaluate the line integral of b around C in the direction indicated by 


155 



the fingers of the right hand when the thumb points in the direction 
of the positive x 3 axis. 

11. Prove Equations (58.5) and (58.6). 

12. A vector field b has continuous first derivatives in a closed 
region V. On the bounding surface S of V, b is normal to S. Prove 
that 

f VXb dV = 0. 
v 

13. If / and b are fields with continuous first derivatives in a region 
V bounded by a surface S, prove that 

nxb^= f /V Xb dV+ f V/Xb dV. 
v v 

14. If b and c are irrotational vector fields, prove that b X c is 
solenoidal. 

15. If b = x 2 x 3 i 1 -\-x 3 x 1 i 2 -\-x 1 x 2 i 3 , show that b is solenoidal, and 
find its vector potential. 

16. Show that the vector field 


If 


* 3*1 ~H *2*2 


is solenoidal in any region which does not contain the origin. 

17. If r, 0 and z are cylindrical coordinates, show that V9 and V In r 
are solenoidal vectors, and find their vector potentials. 

18. If V -b = 3> and V Xb = f, prove that 



dx 2 ‘ + ’ dx 3 


and derive similar equations involving b 2 and b 3 . 
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Chapter VI 


TENSOR ANALYSIS 

61. Introduction. Tensors are mathematical or physical concepts 
which have certain specific laws of transformation when there is a 
change in the coordinate system. As we shall see, vectors are just one 
of the many types of tensors. 

Tensor analysis is a study of tensors. It has many applications. One 
of these is in the field of classical differential geometry of the curve and 
surface in our ordinary space, as well as the more significant generaliza¬ 
tions of this geometry to spaces of higher dimensionality; this generali¬ 
zation is often referred to as Riemannian geometry. Another applica¬ 
tion is to mathematical physics. Here tensor analysis permits us to 
express easily in terms of curvilinear coordinates the fundamental 
equations of the various subjects such as hydrodynamics, elasticity, 
electricity and magnetism. It also aids in the formulation of the 
natural laws of mathematical physics, since such laws when expressed 
in terms of tensors are independent of any one particular coordinate 
system. 

62. Transformation of coordinates. Let us consider a set of curvilinear 
coordinates which we denote by the symbols z 1 , Z 2 and -j 3 . We use 
the supercripts to agree with a convention to be introduced later. 
We shall also refer to these coordinates by writing z' f = 1, 2, 3). 

Let us write 

(62.1) =f(z\z\z?), 

where the three functions / r (r = 1, 2, 3) are single valued and 
differentiable for some range of values of z 2 and z 3 - These equa¬ 
tions represent a coordinate transformation to new curvilinear coordi¬ 
nates z’ 1 , z ' 2 and ^' 3 . The Jacobian /' of a transformation such as 


157 



this one was defined in § 49. We express it conveniently here by writing 


(62.2) 



the ranges for r and s being 1,2,3. We assume that /' does not vanish, 
whence, as mentioned in §49, Equations (62.1) may be solved to 
yield 


(62.3) 


€ = g r {z'\ z'\ z’ 3 ). 


From Equation (62.1) we have 

(62.4) dz' r = £ dz? (r = 1, 2, 3). 

We now introduce two conventions, as follows: 

Range convention. A small latin suffix which occurs just once in a 
term is understood to assume all the values 1, 2, 3. 

Summation convention. A small latin suffix which occurs just twice 
in a term implies summation with respect to that suffix over the range 
1,2,3. 

Because of ihese conventions, we may now write Equations (62.4) 
in the form 

(62.5) dz' T = ^ dzf, 

the repeated suffix s implying the summation. In this equation, the 
suffix r is called a free suffix, since we obtain from (62.5) a different 
equation for each value of r; on the other hand, the suffix s in (62.5) 
is called a dummy suffix, since a change of both of the suffixes s into 
any other latin suffix does not alter Equation (62.5). 

The Kronecker delta 8„ was introduced in § 47. We shall require 
it in the present chapter, but shall find it convenient to denote it by 
the symbol §(, so we have 

K = 1 if r = r 

(62.6) 

= 0 ifr^r. 
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We note the identity 
(62.7) 


Si = 


dz r 3z' s 
Witt’ 


which is true because its right side is equal to 3^/3z‘, which, because 
of the independence of the coordinates z r , satisfies the relations 



= 0 


if r = t, 
if r # t. 


Of course the primed and unprimed coordinates in Equation (62.7) 
can be interchanged. 

Throughout this section we have considered three curvilinear coor¬ 
dinates z 1 , Z 2 , Z? in space. The considerations here can be applied 
equally well to the case of two curvilinear coordinates z 1 , Z 2 on a 
surface, the only modifications being in the range and summation 
conventions above, where it becomes necessary to assume that latin 
suffixes have the range 1, 2. This same modification, when applied to 
the rest of this chapter, yields results applicable to two-dimensional 
problems. In an analogous manner, by the assumption that latin 
suffixes have the range 1, 2, ■ • •, N, where X is any positive integer, 
the results in this chapter can be generalized to yield the Riemannian 
geometry of the W-dimensional Riemannian space. There are many 
excellent books on this subject. 1 


63. Contravariant vectors and tensors. Let X be a general point in 
space with curvilinear coordinates z r ■ Let A r denote a set of three 
quantities associated with the point X, such as for example compo¬ 
nents of a velocity. We now introduce a second curvilinear coordinate 
system relative to which the coordinates of X are z' T , and denote the 
transforms of A' by A". We then make the definition: a set of quan¬ 
tities A’ associated with a point X is said to be the components of a contra¬ 
variant vector if they transform , on change of coordinates, according to the equation 

1 See, for example: J. L. Synge and A. Schild, Tensor calculus, University of 
Toronto Press, Toronto, 1949. 
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where the partial derivatives are evaluated at X. 

Because of Equation (62.5) we see that the quantities dzp are the 
components of a contravariant vector. Also, we shall see later on that 
when the transformation is from rectangular cartesian coordinates to 
rectangular cartesian coordinates, the components of a vector as defi¬ 
ned in § 6 have the law of transformation (63.1) and are hence also 
the components of a contravariant vector. 

We shall now define contravariant tensors. Let A” be a set of nine 
quantities associated with the point X with curvilinear coordinates 
<?, and let A" transform into A" 1 when the coordinates are transfor¬ 
med to z". We then have the definition: a set of quantities A' 5 is called 
the components of a contravariant tensor of the second order if they transform 
according to the equation 

dz’ r dz' s 

(63.2) ^=gr§r^“- 

Contravariant tensors of higher order are defined analogously. Contra¬ 
variant vectors are often called contravariant tensors of the first order. 


64. Covariant vectors and tensors. The definition of covariant vectors 
is as follows: a set of three quantities A r associated with a point X is said 
to be the components of a covariant vector if they transform on change of coordi¬ 
nates, according to the equation 

(64.1) A\ = ^r r A s . 

As an example, let us consider a function O^ 1 , f, Z 3 )- We have 

dz' r ~Jz'W r ' 

whence we note that the three expressions are the components 

of a covariant vector. 

The definition of a covariant tensor of order two is as follows: 
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a set of nine quantities A rs is said to be the components of a covariant tensor 
of order two if they transform according to the equation 


(64.3) 


A' 


rs — 


dz‘ dz u . 
dz" dz' 1 


Covariant tensors of higher order are defined analogously. Covariant 
vectors are often called covariant tensors of order one. 

In § 63 we have used superscripts in writing symbols denoting contra- 
variant tensors, and in the present section we have used subscripts in 
writing symbols denoting covariant character. This convention will 
be followed throughout the present chapter. 


65. Mixed tensors. Invariants. We can also define tensors whose law 
of transformation involves a combination of both contravariant and 
covariant properties. Such tensors are called mixed tensors. Thus, 
for example, a set of twenty-seven quantities A', is said to be the components 
of a mixed tensor of the third order with one contravariant suffix and two co¬ 
variant suffixes if they transform according to the equation 


(65.1) 


'r — ifL idL 
■* “ dz u dz' 1 dz' 1 


It should be noted that (65.1) represents twenty-seven equations, 
and that the right side of each of these equations consists of the sum 
of twenty-seven terms. 

A single quantity A is said to be an invarinat if it transforms according 
to the equation 

A' =A. 


Invariants may be called contravariant tensors of order zero, or co¬ 
variant tensors of order zero. An example of an invariant is the tem¬ 
perature in a fluid. 

We shall now prove that the Kronecker delta 3) has the tensor char¬ 
acter indicated by its suffixes, that is, that it is a mixed tensor of the 
second order. We must establish the equation 


(65.2) 


s'; 


dzf djf_ 
dz‘ dz' 1 
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Now 


3jf_ # ff s , » 3z‘ , 
3z! s “ - 3z' s 1 + 3z s 


3z 3 

4 . — S' 

+ a?" 31 


When t is given the values 1, 2 and 3, the right side of this relation 
reduces respectively to the expressions 

3zf 3zf 
3 Z ' S ’ 3 Z ' S ’ 3z' s ‘ 

Thus we may write 

3z u s , 3z‘ 

3z! s “ 3 Z ' S ’ 

and so the right side of Equation (65.2) reduces to 

3zTdz 1 
3z‘ 3z! s ' 

But this is equal to the left side of Equation (65.2) because of Equation 
(62.7). 


66. Addition and multiplication of tensors. If we have two tensors of the 
same order and type, their sum is defined to be the set of quantities 
obtained by the addition of corresponding components of the two ten¬ 
sors. Thus the sum of the tensors A', and B’ si is a set of twenty-seven 
quantities C’ st given by the relations 

(66.1) C r si = A r sl +B' st . 

It is easily proved that the sum of two tensors is a tensor of the same 
order and type as the two tensors added. 

There are two products of tensors, called the outer product and the 
inner product. The outer product of two tensors is defined to be the 
set of quantities obtained by multiplication of each component of the 
first tensor by each component of the second tensor. Thus, for exam¬ 
ple, the outer product of the tensors A r s and B r Jt is the set of 243 quanti¬ 
ties C," given by the relations 

(66.2) C’f = A\F m . 

In writing this equation, we have been careful to keep each particular 
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suffix at the same level on both sides of the equation. When this 
convention if followed, it is easily proved that the outer product of 
two tensors has the tensor character indicated by the number and 
positions of its suffixes. Thus, for example, in Equation (66.2), C'f is 
a mixed tensor of the fifth order, with two contravariant suffixes and 
three covariant suffixes. 

We now introduce an operation called contraction. It consists in 
identifying a superscript and a subscript of a tensor. Thus, for example, 
if from the tensor A r sl we form the set of quantities B, defined by the 
relation 

(66.3) B t =A’ rt , 

we are performing a contraction. We note that each component 
of B, is equal to the sum of three components of A',. It is easily 
proved that contraction of a tensor yields a tensor Thus, for example, 
B, in Equation (66.3) is a covariant vector. It should be noted that 
we do not perform contractions by identifying suffixes at the same level, 
since such operations do not yield tensors, in general. For example, 
if A', is a tensor, A r ls is not a tensor, in general. 

We now define the inner product of two tensors. To obtain it, we 
form an outer product of the two tensors and then perform a con¬ 
traction involving a superscript of one tensor and a subscript of the 
other tensor. Thus, for example, an inner product of two tensors A r s 
and B rs is the set of quantities C st given by the relation 

C st =A'A, 

Of course the inner product of two tensors is not unique. Since outer 
multiplication and contraction of tensors both yield tensors, the inner 
multiplication of tensors also yields tensors. 

67. Some properties of tensors. One of the most important properties 
of tensors is the following: if a tensor equation is true for one coordinate 
system, it is true for all coordinate systems. To prove this, let us consider, 
for example, a tensor having components A r sl for a specific coordinate 
system z r , and let us suppose that 
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(67.1) A r sl = 0. 

Let z" be any other coordinate system, and let A denote the com¬ 
ponents of this tensor for the coordinate system z' r ■ We must show 
that 


(67.2) 


A" 

A si 


= 0 . 


By the laws of tensor transformation we have 


A';, = 


K r d JL d ^L a* 
dz u 3z' s 3z n vw 


= 0 , 


by Equation (67.1). The proof is similar in the case of other tensor 
equations. 

Now let A r st and B T st denote the components of two tensors for a 
specific coordinate system £ r , and let us suppose that 

A' t = B' t . 


Then, by the previous section, A r st - B r st is a tensor which vanishes 
for the coordinate system z', and hence for any other coordinate sys¬ 
tem z! r we have 

A”* = B” sl . 

Tensors also have the property of being transitive. To explain this 
property, we introduce three coordinate systems z\ z' r and z!”■ The 
property of being transitive is then the following: if a certain set of 
quantities is a tensor for the transformation from coordinates z to 
Z ,T , and is a tensor of the same order and type for the transformation 
from coordinates z' r to z" r \ then for the overall transformation from 
the coordinates z r to z" r the set of quantities is a tensor of this same 
order and type. It is easily proved that a tensor of any order or type 
is transitive. 


68. Tests for tensor character. We shall now demonstrate by some 
examples a useful test for establishing the tensor character of a set 
of quantities. 

Example 1. Let A rs be a set of nine quantities such that A"X S 
is a contravariant vector, where X T is an arbitrary covariant vector. 
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We shall now prove that A rs is a contravariant tensor of the second 
order. Since A"X S is a contravariant vector, we have 

(68.1) A' rs X' s = ~ A lu X u . 

But X„ is a covariant vector, so 


( 68 . 2 ) x u =^X\. 

We now rewrite Equation (68.1), replacing the dummy suffix s on 
the left side by v, and substituting for X u on the right side from (68.2), 
to obtain the relation 


or 

(68.3) 


A ,rv X\ 


dz ,r A , u 
3z‘ d z u 


dz” dz v 

(A'"-££a*)X'. = 0. 


Since X r is arbitrary, so is X' v . Hence the expressions in the brackets 
in Equations (68.3) vanish, whence we conclude that A” is a contra¬ 
variant tensor of the second order. 

Example 2. Let A rs be a set of quantities such that A rs X r X s is 
an invariant, where X' is an arbitrary contravariant vector. We shall 
now prove that A rs is a covariant tensor of the second order provided 
it is symmetric in all coordinate systems, that is, provided we have 
for every coordinate system relations of the form A„ = A sr . We 
proceed much as in Example 1. Since A rs X r X‘ is an invarint, we have 

(68.3) ' A' rs X"X ,s =A„XX. 

Since X r is a contravariant vector we can express the X T and X s in 
this equation in terms of X", obtaining 


A' rs X' r X' s 


dz T dzf 
3z‘ dz" 1 


X' l X' u 


We now replace the dummy suffixes r and r on the left side by t and u, 
whence we obtain the equation 

(68.4) b, u X' t X' tt = 0, 
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where 


(68.5) ^ = A' lu -~^^- u A rs . 

Equation (68.4) can be written in the form 

(68.6) b n (xy+b M ( xy+b 33 (xy+(b M +b n ) x'*x' 3 

+ (*«+*») X'°X n +(b l3 +b 31 ) X n X' ! = 0. 

Since X’ is arbitrary, so are the three quantities X". Hence the coef¬ 
ficients in Equation (68.6) must vanish, which leads to the equation 

(68.7) b tu +b ul = 0. 


From (68.5) we then have 


A' tu + A' u , = 


dz! dz s 

WlkX 


A rs + 


dz! dz s 

W u H7‘ 


In the last term on the right side, we interchange the dummy suffixes 
r and 5, obtaining the relation 


( 68 . 8 ) 


_dz r dz! . . . . 

tu i A. ut dz >u ' w * 


Since we are given that A rs = A jr , A' lu = A' an then Equation (68.8) 
reduces to the form 


a' — a 

A — dz’ 1 3z ' u • 

which establishes the tensor character of the set of quantities A rs . 


69. The metric tensor. Let x, be the three rectangular cartesian 
coordinates. Then the distance ds between two adjacent points is 
given by the relation 

(69.i) ( d s y = (d Xl y+(dx 2 y+(dx 3 y = dx,dx r . 

If we introduce curvilinear coordinates X defined by relations z? = 
* 2 , *,) or x r = g r (z\ z\ Z*), then 
. dx r 
dXr = dz? d ^’ 

and substitution in Equation (69.1) yields an expression of the form 
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( 69 . 2 ) (dsy = b^dzr+b^dey+bjdzr+bjz'dz* 

+b al dz 3 dz 1 +b u dzdz\ 

where b rs are six functions of z 1 , z 2 , We now define nine quantities 
g rs by the relations 

gll J b 2 2 - g 22 , ^33 — £33 , 

(69.3) 

2^23 ^23 &2 J ^31 = ^31 = £l3 3 2^12 gli = gil ' 

We may now express (69.2) in the form 

{dsy = g l i(dz i y+g M {dz‘y+gjdz*y+{g i3 +g 3 ,) dydy 
+ (g n +gi 3 ) dz'dz 1 + {g 12 +g n ) ^‘^ 2 , 
or 

(69.4) [dsy = g jydy. 

Let us observe the right side of this equation. It is an invariant, since 
ds is an invariant; also dz‘ is an arbitrary contra variant vector; also, 
from Equations (69.3) defining g rs , we see that g rs is symmetric. 
Hence, by the test for tensor character presented in Example 2 of 
§ 68 we can conclude that g rs is a covariant tensor of the second order. It 
is called the metric tensor. 

70. The conjugate tensor. Let g denote the determinant whose ele¬ 
ments are the components of the metric tensor. Then 

(70.1) g = \g„\- 

In the expanded form of this determinant, the coefficient of any one 
element g rs is called the cofactor of g rs . We denote it by the symbol 
A". We note in passing that the minor of g, s is equal to 

{-iy +s a ". 

We shall now prove that 

(70.2) J?„A"=^, 

(70.3) g„ A" = g%. 

Proof of Equation (70.2). Let us suppose that, when we write the 
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determinant |g„|, the subscript r varies over the rows, while the sub¬ 
script s varies over the columns. We now consider Equation (70.2) 
when s = t = 1. The right side is equal to g. The left side is 

£riA ri = g u A“+foA“+^A“; 

which is just the expansion of the determinant g by elements of the 
first column. The proof of Equation (70.2) when r=s= 2 and 
r=j = 3 is similar. Now let us consider Equation (70.2) when s= 2 
and t = 1. The right side is equal to zero. The left side is 

g» Arl = gu A n +g n A 21 +^ 32 A 31 , 

which is just the expansion of a determinant obtained from g by a 
replacement of the elements of the first column by those of the second 
column; in this determinant the first two columns are identical, and 
so the value of the determinant is zero. The proof of Equation (70.2) 
in all other cases when r and s differ is similar. 

Proof of Equation (70.3). This proof is quite similar to the proof of 
(70.2), differing from it only in that the various expansions of deter¬ 
minants are by elements of rows rather than columns. 

We now define the conjugate tensor g rs by the relation 



Then Equations (70.2) and (70.3) yield the expressions 

(70.5) g„/‘ = 8J, 

(70.6) g„£ = Sj. 

We shall now prove that g” is a contravariant tensor of the second 
order. Even though Equation (70.5) indicates that g rs g n is a mixed 
tensor of the second order, and g rs had known tensor character, we 
cannot make any direct deductions from this as to the tensor character 
of g rt , since we cannot assign g rs arbitrarily without changing g rs . Let 
X r be an arbitrary contravariant vector. If we write 
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g*X = Y„ 



then Y, is an arbitrary covariant vector. Hence we have 

g p, Y r = fgJP = BfX 1 = X*, 

and so g pr Y r is a contravariant vector, Y r being an arbitrary covariant 
vector. Thus, by the test for tensor character presented in Example 
1 of § 68, we can conclude that g pr is a contravariant tensor of the 
second order. 


71. Lowering and raising oj suffixes. In order to avoid ambiguity 
in certain steps to follow, we shall adopt the convention of placing 
but one suffix in any one vertical line. We shall use dots to denote 
vacant spaces in the rows of suffixes. Thus, for example, we might 
write A[ s to denote a mixed tensor of the second order. Other exam¬ 
ples are B r '* n C ",. 

If we form inner products of any given tensor with the metric 
tensor or the conjugate tensor, we obtain a number of new tensors. 
It is customary to use the same principal letter in writing these tensors. 
For example, from a contravariant tensor A”, we can generate another 
second order tensor A[ s as follows: 

a:, =g sl r. 

In appearance, this new tensor differs from the original one only in 
that the suffix in the second position has dropped from the upper level 
to the lower one. We say that we have lowered a suffix. Other tensors 
obtained from A n by this lowering process are the following: 

A; s = g ri A\ A rs =g rlgsu A‘\ 

In an analogous fashion, the inner multiplication of a tensor by the 
conjugate tensor results in a raising of a suffix. Thus for example, 
from a given tensor B rs we can generate other tensors by this raising 
process, as follows: 

B[ s = fB ts , B\ ! =g l B rl , Er=ffB m . 

If a certain subscript is raised, and is then lowered, it is easily seen 
that the final tensor is the same as the original one. The convention 
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of writing no two suffixes in any one vertical line was introduced solely 
to ensure this property. 

72. Magnitude of a vector. Angle between two vectors. Let A r be a 
contravariant vector. Its magnitude A is defined by the relation 

(72.1) A =V]J 7 A s . 

From this relation it readily follows that 

(72.2) A = VAj 1 = Vg" A,A s . 

A unit vector is one whose magnitude is equal to one. We note that, 
if the coordinates are rectangular cartesian, then g rl = 8„, and 

(72.1) then reduces to the familiar form 

(72.3) A =V{A 1 y+(Ay+(Ay. 

Let X’ and Y r be two unit vectors at a point. The angle 0 between 
them is defined by the relation 

(72.4) cos 0 = g„X r Y ! . 

If the coordinates are rectangular cartesian, this equation reduces to 
the familiar form 

(72.5) cos 0 = X 1 Y 1 +X i Y , +X a Y\ 


73. Geodesics. The parametric equations of a curve C may be 
written in the form 


(73.1) 




the arc length s of the curve being used as the parameter. Because of 
Equation (62.5) we see that df/ds are the components of a contra¬ 
variant vector. We denote it by p r , so 


(73.2) 



Further, by a definition of the previous section, the magnitude p of 
this vector is given by the relation 
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r-i-rr-i-tt 


dz m dz n (dsY __ 

W~ 

Thus p r is a unit vector. It is called the unit tangent vector of the 
curve C. 

A geodesic may be defined as the curve of shortest length joining 
two points. In our three-dimensional space, the geodesics are straight 
lines. If we consider surfaces, which are of course two-dimensional 
spaces, the geodesics are not necessarily straight lines. For example, 
in the case of a spherical surface, the geodesics are the great circles, 
that is, those circles on the sphere whose centers coincide with the 
center of the sphere. 

Let X and Y be two points. The distance L between them, measured 
along some curve, is given by the line integral 

y y _ y _ 

(73.3) L =J ds = jVg mn d^dz 11 = JVw ds , 

xx x 

where 


(73,4 w =g mn p n p\ 

According to the Calculus of Variations, L is an extremum if the path 
joining X and Y is such that 


(73.5) 


d j 3w\ dw 
ds \dp j dg 


The actual derivation of these equations is beyond the scope of this 
book. 

Now in is a function of p r and g mn which are known functions of zf. 
Thus w =w(z\ Z\ e,p\p\p % ). The two partial derivatives in 
Equation (73.5) are computed as though the quantities p r and zf 
are independent. Thus 


^r=&- W + K 

= gmp" + gmrp m = gmp n + grmp m 

- gmp" + gmP" = 2g rn p", 
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dz r dz ,rr 

Hence equations (73.5) become 

/75 c\ „ dp | dg n i 3g mn _ _ 

(73.b) im-jp+^Pp-s-^-pP — U. 

But, by an interchange of the dummy suffixes m and n we have 

(73.7) j^p m p n = d -j~p n P m - 

Substituting for one half of the middle term in Equation (73.6) from 
Equation (73.7), we then obtain the relation 

(73.8) g m ^ + [mn, r]p m p n = 0, 
where 

(73.9) 

The quantities \mn, r] are the Christoffel symbols of the first kind. We 
note that 

(73.10) \mn,r\ = \nm,r\ . 

If we multiply Equation (73.8) by/', we obtain 

(73.11) g rl gm d £ + F‘ mn pT = 0, 
where 

(73.12) F' mn =g ,r [mn,r]. 

The quantities F‘ mn are the Christoffel symbols of the second kind. 
We note that 

(73.13) F‘ mn =F‘ m . 

The first term in Equation (73.11) reduces to 

dp” dp‘ d 2 z‘ 



Thus we may write the differential equations of a geodesic in the form 


(73.14) 


d*f Ft d^_d? 

ds 2 + mn ds ds 


= 0 . 


74. Transformation of the Christojfel symbols. Let z r and z’ r be two 
sets of curvilinear coordinates. In terms of these coordinates the equa¬ 
tions of the geodesics are 


(74.1) 

d l Z r pr dz m dz n n 
ds 2 + ds ds ~ ’ 

and 


(74.2) 

<Pz” , p,s dz’ p d’z? n 
ds 2 + * pq ds ds - U - 

But 

dz r dg dz p 


Ts ~Jff p 1T' 


d 2 z r dz’ d 2 z' p d 2 zd dz' p dz q 

ds 2 ~ d z ' p ds 2 + dz' p dz!" ds ds 

Substitution in 

Equation (74.1) then yields 


dz r d 2 z' p ( d 2 sf 
dz' p ds 2 + \ dz p dz' q 


UF' \dzfdjf _ 

+ Pmn dz' p dz’ q J ds ds ~ 


We now multiply this equation by dz' s /dz r , obtaining the relation 
dfzf [djf_ d 2 Z r djf_ djff djf \ dzf dzf _ 

ds" \dz r dz' p dz! q mn dz T dz p dz q ) ds ds 


Comparison of this equation with (74.2) then yields 


™ _ h" dz m 2 Z n ^ , d z ' s d‘z r 
tq d z r d z ' p 3 z ' q mn + dz’ 3 Z ' P d z " 


This is the equation of transformation of the Christoffel symbols of the 
second kind. We note that these symbols are not tensors, but would 
be if the last term on the right side were missing. 
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75. Absolute differentiation. Let A p be a covariant vector defined 
over a curve C with equations £ =/ r («), u being a parameter on C. 
Then A p = A p (u). The absolute derivative of A p along C is 


(75.1) 


ftfjg _ dAp ps a dz q 
8u du Pq ' du 


We shall now prove that 8A p /8u is a covariant vector. We have 

4' d jl. 

1 du 


(75.2) 

W p 

8u 

3A p .-./j 

= l^~ Fpq - 

Now 

dA' p 

d / 3z r . ' 

(75.3) 

du 

“ du \3z p \ 



3z r dA, 



= d z ' p du + 


3*z r dz' q , 

J.. 


dz' p dz' q du 


Also, because of Equation (74.3), the second term on the right side of 
Equation (75.2) satisfies the relation 

rv, r _ W i? K. F' 4' 

H ‘ du - di 3z' p 3z' q " ' du 

dzff 3 y , 

+ d* dz! p 3z' q J du 


(75.4) 



F’ A 

x mn lx r 


dz n 

du 


+ 


d‘z r dz’ q 
dz' p 3z' q r du ' 


Substitution from Equations (75.3) and (75.4) in (75.2) then yields 

8A' P _ 3zff_ / dA„ _ \ 

8u ~ dz ,p \du " 7 du) 

_ dz m 8A m 
3z ,p S u ' 

This completes the proof. 

The covariant vector A r is said to be propagated parallelly along 
Cif 

(75.5) g - 0. 
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In this case the components A r satisfy first order differential equations, 
and can hence be assigned arbitrarily at any one point on C. When 
the space is three dimensional and the coordinates are rectangular 
cartesian, Equations (75.5) reduce to dA r jdu=Q, so that A r are 
constant along C. 

If A r is a contravariant vector defined on a curve C, its absolute 
derivative along C is 


(75.6) 


8J' _ tf 
S u ~ du 


,d? 

du 


We could prove that 8A r /Su is a contravariant vector in a manner 
similar to that used above for SAJSu. However, it is more convenient 
to proceed as follows. Let B r be any vector propagated parallelly 
along C. Then 


dB r 

du 


— p m 

— rn ■*-' m 


djd__ 

du 


But A m B m is an invariant, and hence so is 

!(**)«£* + *£ 

= d -f B r + A' 

au du 

t£ +F » A .'X\ B , 

du du E ■ 


Since B m can be assigned arbitrarily at any one point, the coefficient 
of B m here, which is SA m /Bu, is then a contravariant vector. 

The contravariant vector A T is said to be propagated parallelly 
along C if 



Let A mn , A m n and A m „ be any second order tensors. Their absolute 
derivatives along a curve C are 


(75.7) 


_ FT — 
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(75.8) 

(75.9) 


8u 

Su 


J pm _ pp Jpn dz?_ 

du V tpqAn du tnqA ' p du 


dA m 

du 


- FL A t 


^ - pp A 

r nq A mp du 


du 


The patterns exhibited by the suffixes in these equations should be 
noted. Just as in the case of the absolute derivatives of the vectors, 
each of the above three absolute derivatives has the same tensor 
characters as the original tensor. We can prove this in a manner 
similar to that used above. For example in the case of 8A m J8u, we 
form the invariant A"!„X m Y", where X m and Y n are any vectors pro¬ 
pagated parallelly along C. 

The absolute derivative of any tensor of higher order is defined 
similarly, and has the tensor character of the original tensor. For 
example, we have 


§ A r sl 

8 u 


ds 


+ F' pq A p sl a f - F% A'.„ j^-F%A : sp 


dt 


If A is an invariant, we make the definition 


(75.10) 


SA _ dA 
S u du 


It can be proved in a direct manner that 


(75.11) 


^§mn _ q 

S U 



8 u 


= 0 . 


The proofs of these are left as exercises for the reader (Problem 28 at 
the end of this chapter). 

The rule for the absolute derivative of a product of two tensors is the 
same as for the ordinary derivative of a product. Thus, for example 

(75.12) + 


This can be proved in a direct manner which is somewhat lengthy. It 
can also be proved concisely by the use of special coordinates called 
Riemannian coordinates. 
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76. Covariant derivatives. Let A r be a covariant vector defined over 
some region V in space, and let C be a curve in V. If a is a parameter 
on C, then 


8 Ap _ dA p dz? 

8 u du pq ’ du 



fla, 


dz q 

du 


The expression on the right side is a covariant vector, and dz q /du is 
a contravariant vector which may be assigned arbitrarily. Hence the 
coefficient of dz q /du is a covariant tensor of the second order. We 
denote it by A p , q and call it the covariant derivative of A p , so we have 

A A 

(76.1) A PI , = ^-F' P ,A S . 


In a similar manner we can arrive at the covariant derivatives of ten¬ 
sors of all orders and types. We have, for example, 


(76.2) A>„=^+Ft t A', 

(76.3) A p . ms = ~ A\, + Ft A\ v - F‘ v A». tr - F‘ s A p . qt . 


These two covariant derivatives have the tensor character indicated 
by their suffixes. 

Of course, we have from Equations (75.10) and (75.11,), 

dA 

(76.4) A l} = 

(76.5) =°» S ^=°- 


Because of these relations, the operation of raising and lowering suf¬ 
fixes can be permuted with the operation of taking the covariant 
derivative. Thus, for example, we have 

grsA f |( = {grsA 1 ) |, = A, 

It can be shown that covariant differentiation has the same product 
rule as ordinary differentiation. 


77. The curvature tensor. Now 



Since A r |„, is a covariant tensor of the second order, then 


(77.2) 


(d r |, n ) |„ 


dt 


n (d r |m) Frn Aq\ m A r \q . 


By use of Equations (77.1) and (77.2), we can arrive after a straight¬ 
forward but lengthy calculation at the relation 

(77.3) A^ mn -A r ^„ m = R‘ rm „A s , 


where 

/77 4.\ b< -fL Fi _ F s \ pp F* - fp Fi. 

\ / /• */ *'--rmn — rn * 1771 ' 2 1 rm pn • 

From Equation (77.3) it follows by the tests for tensor character that 
R’ rmn has the tensor character indicated by its suffixes. It is called the 
mixed curvature tensor. 

The covariant curvature tensor is 


(77.5) R rmn =g rl R[„ 

This tensor plays an important role in mathematical physics. For 
any coordinate system in three dimensional space, this tensor vanishes, 
since it vanishes for rectangular cartesian coordinates. On the other 
hand, if for any one coordinate system on a surface this tensor vanishes, 
then there exists for this surface a curvilinear coordinate system such 
that 

g rl = 1 if r = s 
= 0 if r = s. 


78. Cartesian tensors. A tensor is said to be cartesian when the trans¬ 
formations involved are from one set of rectangular cartesian coordi¬ 
nates to another. In § 47 the special case of such transformations was 
considered when the transformation is a rotation about the origin. 
In the general case, the transformation is a rotation plus a translation. 
If we superimpose a translation on the rotation considered in § 47, 
the equations of transformation are 

(78.1) 
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% s — ^sr^r H“ ^s > — tlrsX r -f* d s , 



where a„ are the constants considered in § 47 satisfying the orthogon¬ 
ality conditions 


(78.2) 


^rt a sl $ rs , &t T Q-is 8 rs . 


Also, a s and a' s are constants such that a', = - a TS a.,. Just as in §47, 
we have the relations 


(78.3) 


dx s _ dx r 

d7~ asr = dl s 


The Jacobian I of the transformation is 

T dx r 

I= -W- 

Hence by the rule for the multiplication of determinants, we have 
(78.4) F = \a„\-\a mn \ = \a, r a tn \ = |S m | = 1, 

so / = ± 1 • 


Theorem 1. For cartesian tensors there is no distinction between 
contravariant and covariant character. 

Proof. Because of Equation (78.3) it follows that the laws of trans¬ 
formation of contravariant components and covariant components 
are the same. Further 

P = S a = 1 a mn = S 

Smn KJ mn j 5 } 6 u mm 

so that the raising or lowering of a suffix does not change the values 
of the components. 

Because of this theorem, when dealing with cartesian tensors we 
do not need both superscripts and subscripts, so subscripts will be used 
exclusively. 

In § 6 we introduced the orthogonal projections of a vector on rec¬ 
tangular cartesian coordinate axes, calling these projections the com¬ 
ponents of a vector. Throughout the rest of this book we shall refer 
to these as the physical components of a vector for rectangular cartesian 
coordinates. 

Theorem 2. The physical components of a vector for rectangular 
cartesian coordinates constitute a cartesian tensor of the first order. 
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Proof. The nine constants a u in Equations (78.1) are the cosines 
of the angles between the axes of two rectangular cartesian coordinate 
systems. Hence, if b is a vector with physical components b T and V, 
for these two systems of coordinates, then just as in § 47 we have 

b s a sr b r 



This is the law of transformation of a cartesian tensor of the first order, 
so the proof is complete. 

When the coordinates are rectangular cartesian,, the Christoffel 
symbols vanish, and so the absolute derivative becomes the ordinary 
derivative and the covariant derivative becomes the partial derivative. 
For example, SbJSu becomes db r \du , and b^ s reduces to db T \dx s . 
Thus we conclude that the directional derivatives of the physical 
components of a vector for rectangular cartesian coordinates constitute 
a cartesian tensor of the first order, and the nine partial derivatives 
of the physical components of a vector constitute a cartesian tensor 
of the second order. 

79. Oriented cartesian tensors. Equation (78.4) shows that the Jaco¬ 
bian I of a transformation from one set of rectangular cartesian coor¬ 
dinates to another is equal tO'either plus one or minus one. The former 
case arises when the transformation is between coordinates whose 
axes have the same orientation (both right-handed or both left- 
handed), and the latter case arises when the orientations are opposite. 
A set of quantities is said to constitute an oriented cartesian tensor if 
it is a cartesian tensor when I = 1 and is not a cartesian tensor when 
I — - 1 . 

In two-dimensional problems, suffixes have the range 1, 2. For 
such problems we introduce a permutation symbol c rs defined as follows: 

(79.1) c u = c 22 = 0, c l2 = 1, c 21 = -1. 

In the three dimensional case, the permutation symbol c rst has the defi¬ 
nition 
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c rsl — 0 if two suffixes are equal 

(79.2) = 1 if (rst) is an even permutation of (123) 

= -1 if (rst) is an oad permutation of (123), 

a single permutation of (rst) being an interchange of any two of r, s 
and t, and an even or odd permutation meaning an even or odd 
number of single permutations. Thus, for example, 

C 123 = C 312 = 1 J ^321 = — 1 , £u3 = 0 ' 

In Theorem 4 of the next section we shall see that both of these per¬ 
mutation symbols are oriented cartesian tensors. 

We can now express in tensor notation many of the formulas and 
equations of the earlier chapters of this book. Thus, the scalar and 
vector products of two vectors a and b are respectively 

a r b, •, c rsl a s b t . 

We note that this vector product is an oriented cartesian tensor. The 
scalar triple product a ■ (b X c) is 

C rst a r b s C t . 

The expressions V/, V • b and V X b become 
df_ dK db, 

dx r ' dx r ’ C "‘ dx; 

respectively. The differentiation formulas (48.1)—(48.11) may also be 
expressed in this notation; when this is done (Problem 36 at the end 
of this chapter), some of the formulas become trivial, and the truth of 
others follows at once from the identity 

(79.3) c rmtl c rst = 

the proof of which is left to the reader (Problem 34 at the end of this 
chapter). 


80. Relative tensors. The Jacobian / of a transformation from one 
set of curvilinear coordinates z' to another set z" is given by the 
relation 


/ = 


dj_ 

dz! s 
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We now introduce the following definition: a set of quantities A' is a 
relative contravariant vector of weight W, if it transforms according 
to the equation 


A" 


d z ‘ A 


rw 


W being an integer. Relative contravariant tensors of higher order 
are defined analogously, as are relative covariant tensors, relative 
mixed tensors, and relative invariants. For example, A[ st is a relative 
mixed tensor of weight W if it has the law of transformation 


A\ 


Am TW 

d Z m d Z " dz’ 1 " f ‘ 


A relative tensor of weight one is called a tensor density. 

When it becomes necessary to distinguish between the kind of 
tensors considered heretofore and relative tensors we refer to the former 
type as absolute tensors. The following properties of relative tensors 
are easily established: (i) If a relative tensor vanishes in one coordinate 
system, it vanishes in all coordinate systems, (ii) The transformation 
of relative tensors is transitive, (iii) The sum of two relative tensors 
of the same order and weight is a relative tensor of this same weight, 
(iv) Both the inner and outer products of two relative tensors yield 
relative tensors whose weights are equal to the sum of the weights 
of the two tensors, (v) Contraction does not change the weight of a 
relative tensor. 

Theorem 1. The determinant g — |g„| is a relative tensor of weight 
two. 

Proof. Because of the rule for the multiplication of determinants, 
we have 


( 80 . 1 ) g = | g'r. 


which completes the proof. 


dz 1 de_ 
d Z ' r dz" g,u 


dz‘ 


dz u 

dz" 


dz" 


- Tg, 
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Let us now introduce a covariant permutation symbol c rst and a 
contravariant permutation symbol c rs ‘, both defined numerically 
just as in Equation (79.2). 

Theorem 2. The covariant permutation symbol is a relative covari¬ 
ant tensor of weight -1. 

Proof. We must show that 


(80.2) 


£ rst 


dz' r dz s 3z n Cmnp 


Let us denote the right side of this equation by I 1 . Then 


(80.3) 


_ d ? 3 ** 

- dz n 3 Z ' 2 dz' z Cmnp ‘ 


But the expression on the right side here is just the expansion of the 
determinant I, so 0 123 = I. In a similar manner, we see that each 
component of ®„, is equal to a determinant. If the suffixes on any 
one component are subjected to a single permutation, two rows or 
columns of the corresponding determinant are interchanged, and so 
the sign of the component is changed. We thus have the results: 
(i) = 0 if two suffixes are equal, (ii) ®„, = / if (rst) is an even 

permutation of (123), (iii) = -/if (rst) is an odd permutation of 
(123). Thus we may write 


I C rsi) 


and so the proof is complete. 

Theorem 3. The contravariant permutation symbol is a relative 
contravariant tensor of weight one. 

Proof. This proof is quite similar to that of Theorem 2 above, 
and is hence omitted. 


Theorem 4. The permutation symbols are oriented cartesian tensors. 
Proof. When the coordinates are rectangular cartesian, we see 
from Equation (80.2) that c rsl is an absolute tensor only when / = 1. 
Hence it is an oriented- cartesian tensor. The proof is similar for c rl> . 

Because of Theorem 1 above, we may change the weight of a tensor 
by multiplying the tensor by a power of g. Let us assume that g is 
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positive. We can then construct the absolute permutation symbols as 
follows: 

(80.4) -n« = ^~gcrs„ — 

Because of Theorems 1, 3 and 4 above, we see that i)„ t is an absolute 
covariant tensor, while vf s ‘ is an absolute contravariant tensor. 

Let A'.'.', denote a general relative tensor of weight W. Then g* W AV.\ 
is an absolute tensor. The covariant derivative of A'.'.', is defined to be 

(80.5) Av:. v = g iW {g~ iW A\:.) v . 

We note that A'.'.'.i, has the same weight as A'.'.'. . The absolute deri¬ 
vatives of relative tensors are defined analogously. It can be proved 
that both the covariant and absolute derivatives of relative tensors 
obey the same product rule as do ordinary derivatives. 


81. Physical components of tensors. In Theorem 2 of § 78 it was 
stated that the physical components b r of a vector b for rectangular 
cartesian coordinates constitute a cartesian tensor of the first order. 
We shall now define the contravariant and covariant component of 
the vector b for curvilinear coordinates f. Denoting these compo¬ 
nents by E and B ,, we make the definitions 


(81.1) 



bs, 


B,= 


dx s 


b s . 


It is easily proved that B, — g rs B‘. The quantities E and B, describe 
the vector b in a certain manner, and have the tensor character indi¬ 
cated by their suffixes. 

In § 72 of the present chapter we defined abstractly the magnitude 
of a contravariant vector. According to this definition, the magnitude 
of B r is 

B = Vg rs EE. 


Since the right side of this equation is an invariant, and since S rs is 
the metric tensor of rectangular cartesian coordinates, we have 


B = Vs mn b m b n = Vb m b m = b. 
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Thus the abstract definition of magnitude of a vector given in § 72 
agrees with the definition of magnitude given in § 1. In a similar 
fashion, the abstract definition given in § 72 for the angle between 
two vectors agrees with our physical notions of angle. 

Let us consider a unit vector having components X r for a curvili¬ 
near coordinate system z'- The physical component of B" in the direc¬ 
tion of X r is defined to be the invariant 

= g, s B r X s . 

If 0 is the angle between B' and X r , then by § 72 we have 

B r 

cos 8 = g,s y X s . 

Thus B x — b cos 0, as might be expected. 

The physical components of a vector B r in the directions of the 
parametric lines of the coordinates z r are called the physical compo¬ 
nents of B r for the curvilinear coordinates z r - Let us denote them by 
f? (r ). When the curvilinear coordinates are orthogonal, then as men¬ 
tioned in § 49, we have 

(81.2) ( dsy= (Kdzy+{Kdzy+(Kd Z y. 

If X { [) is a unit vector in the direction of the parametric line of z 1 , then 
dz 1 1 

(81.3) X ( 1 l) =^ = -, X < i, = X& = 0. 

By lowering suffixes, we also have 

(81.4) ^( 1)1 - ^i> ^( 1)2 ^( 1)3 

Thus 

-5(0 = irrB r \0 — B% )r = B'kt. 

In a similar dashion we get 5 (s) and f? (3) , so that the physical compo¬ 
nents of the vector b relative to the curvilinear coordinate system 
Z T are 

B\, B\, B%. 

By lowering suffixes, we may also express these physical components 
in the form 
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(81.5) 


Bj_ B% B 3 

h-y h% h 3 

There is a similar procedure for tensors of higher order. For exam¬ 
ple, let t rs be a set of nine quantities which are the components of a 
cartesian tensor. The stress components of elasticity form such a set. 
We define the tensorial components of this tensor for curvilinear coor¬ 
dinates £ in terms of t„ by the appropriate laws of tensor transfor¬ 
mation. Let us denote these components by the symbols T", T%, 
T r s and T rs . The physical component of this tensor along the direc¬ 
tions of two unit vectors X r and p r is defined to be 7'"x r p J . The 
physical components of this tensor for the curvilinear coordinates £ 
are defined to be the physical components obtained by taking X r and 
[T in the directions of the parametric lines of the coordinates. Let us 
denote these nine physical components by T (rs) . Then 

T( rs) = r-x^x^, 

\' s) being the three unit vectors in the directions of the parametric 
lines of the coordinates. When the coordinates z r are orthogonal, then 
X (1 r } is as given in Equation (81.3), and there are two similar relations 
for Xpj and X ( J). We then get for T, rs) the expressions 


Kr\ 

hAT'\ 

hXT 13 , 

(81.6) Kh x T\ 

, h\T\ 

KKT l \ 

hJifT* 1 

, KhX\ 

h\T'\ 

There are other expressions 

for T (rs) in 

terms of the covariant compo- 


nents T rs , as well as in terms of mixed components. 

Physical components of tensors of higher order can be defined 
analogously. 

82. Applications. We shall consider now the problem of expressing 
the fundamental equations of mathematical physics in terms of quan¬ 
tities pertaining to curvilinear coordinates. 

Dynamics of a particle. Let x r be the rectangular cartesian coordi¬ 
nates of a particle at time t. Its velocity and acceleration are then 
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(82.1) 


v, = 


dx r 

dt’ 


dv r 

a ' ~ li' 


If f is the force acting On the particle and m is its mass, then by 
Newton’s second law we have 


(82.2) m a r =f r . 

Let zd be curvilinear coordinates of the particle, V, A' and F r being 
contravariant component of velocity, acceleration and force for this 
coordinate system. We have by definition 


dg _ 3f_ dx s _ <hf 
dx s Vs 3x s dt dt 
We could define ^4 r similarly in terms of a r , but it is simpler to write 

^ = ?IZ +F r V ^A. 

St dt + ” ra dt 


(82.3) 


(82.4) 


V r 


A r = ~ = 


We could define F r in terms of f r , but it is easier to use the relation 
(82.5) dW =fdx r =F,dz r , 

where dW is the work done in an infinitesimal displacement with 
components dx r and dz r . By substitution in this equation for dx r in 
terms of dz\ we get F, by equating coefficients of dz r - 

We now write the following expressions tentatively for the equa¬ 
tions of motion in terms of curvilinear coordinates: 


(82.6) 


mA r =F r . 


To check these we note that they are tensor equations, and are true 
when the coordinates are rectangular cartesian since in this case they 
reduce to (82.2). Hence they are true for all coordinate systems. 

The mathematical theory of elasticity. Using rectangular cartesian 
coordinates x r , we have the displacement u r , the strain components 
e ts and the stress components T rl . For the determination of these, 
we have the equations 


(82.7) 


— 1 (dUr du\ 

rs ~ z + w 


(82.8) 


3 2 e rn = 3 2 e rm + 3%„ , 

3x s dx,„ 3x,3x„ 3x s 3x n 3xf)x m 
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<82.9) 

(82.10) 


T rs = k, 

3T r 


dx s 


+ %r — 


&mn j 


d i u. 



together with certain boundary conditions. In the above, k rsmn are 
elastic constants, X r is the external force per unit volume, p is the 
density and t is the time. If the body is at rest, the right side of (82.10) 
vanishes. If the body is isotropic, that is, it has no preferred directions 
elastically, Equation (82.9) reduces to 


(82.11) 


T„ — X 08 „+ 2 p«„ 


where X and jj. are elastic constants, and 


(82.12) = 

The quantities in the above equations are all cartesian tensors. 

Let us now introduce curvilinear coordinates z r ■ For this coordinate 
system the quantities corresponding to a,, e„, t rs , k rsmn and X r are 
defined'by use of the laws of tensor transformation, as before. There 
is no difficulty if we use the same principal letter to designate these 
quantities for both coordinate systems. We now wish to convert 
Equations (82.7)—(82.12) to curvilinear coordinate. We write ten¬ 
tatively : 

(82.13) e rs = \(u rV +u Ar ), 

(82.14*) ^mjsm X ^sm\rn ^rm\sn~\~ ^sn\rmj 

(82.15) T„ = k rsmn e m \ 

(82.16) T% + JP = P 

(82.17) T rs = X0£„d-2pe„, 

(82.18) 0=4= u [ v . 

These equations are tensor equations, and reduce to Equations 
(82.7)—(82.12) when the coordinates are rectangular cartesian coor¬ 
dinates. Hence Equations (82.13)—(82.18) are true for all coordinate 
systems, and are the desired equations. (Note that the term d^u'jdf 
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in Equation (82.16) is a tensor; see Problem 45 at the end of this 
chapter.) 

Hydrodynamics. In terms of rectangular cartesian coordinates, the 
fundamental equations for a perfect fluid are 

< 82 - 19 » % + k w = °- 


(82.20) 


3Vr , dVr y 1 3p 

dt dx s Vs r p dxj 


where p is the density, t is the time, v, is the velocity, X T is the external 
force per unit volume, and p is the pressure. We convert these to 
curvilinear coordinates in a manner analogous to that used for elas¬ 
ticity, obtaining the equations 


(82.21) 


dp 


(p» r )ir = 0, 


(82.22) 


dll 

In 


x r --g” 

? 


dp 

a?' 


In a similar manner, we can obtain in terms of curvilinear coordi¬ 
nates the fundamental equations of other branches of mathematical 
physics, such as electricity and magnetism, geometrical optics and 
heat conduction. 


Problems 


1. Prove that 


(■ A„+A„)£z'! = 2A rs z r z s . 

2. If A = A rs z r z\ where A rs are constants, prove that 

3A , . . . d*A . . 

^ rs J ^ sr ‘ 

3. Evaluates), S;S), S;S)S‘ r . 

4. Prove that 

KA h =a m , s;s :a‘, = a“. 

5. By differentiating Equation (62.7) with respect to z us prove that 
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2 2 z' T 3z" dz" dz w d l ? 

2z , h‘ ~ dz u 2z‘ dz'° 2z' w ‘ 


6. Let b r denote the covariant components of a vector for rectangu¬ 
lar cartesian coordinates. Find the covariant components of this 
vector for cylindrical coordinates r, 0, £ in terms of b r and r , 0, z. 

7. If g are curvilinear coordinates and ® is an invariant, do the 


expressions 




constitute a tensor? 


2z r 2zf 

8. Prove that the sum of two tensors A', and B', is a tensor. 

9. Prove that the outer product of two tensors A' and B r sl is a tensor. 

10. If A r st is a tensor, prove that A' rt is a tensor. 

11. In three-dimensional space, how many different expressions are 
represented by A™B n pq G q J ? When each such expression is written out 
explicitly, how many terms does it contain? 

12. Prove that the tensor A' t s is transitive. 

13. Let A' be a set of nine quantities such that A r s X s r is an invariant, 
where X s r is an arbitrary mixed tensor of the second order. Prove that 
A' is a mixed tensor of the second order. 

14. Let A Isi be a set of quantities such that A rsl X‘ is a covariant 
tensor of the second order, where X‘ is an arbitrary contravariant 
vector. Establish the tensor character of A rst . 

15. If g rt = 0 for t ^ r s, prove that 


g n 


—, g 22 = — > 

g 23 = r = g 12 



g* 3 


16. Find the components of g rs for cylindrical coordinates r, 0, z- 

17. Find the components of g rl for spherical polar coordinates 
r, 0, cp. 

18. Let z l and z 2 be plane oblique cartesian coordinates whose 
axes have an angle 0 between them. For these coordinates find g„ 
and g". 

19. Prove that g rs g rs = 3. 

20. Prove that 




2gmn 

2z r 
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21. Using the results of Problem 20, prove that 


F s = 

- 1 rs 


Vg 3? yg ' 


22. Prove that 

\rm, n ] + \rn, m ] = [mra, r] = F‘ mn . 

23. Compute the Christoffel symbols for cylindrical coordinates 
r, 0, £. 

24. Compute the Christoffel symbols for spherical polar coordinates. 

25. By using rectangular cartesian coordinates and Equations 
(73.14), show that the geodesics in three-dimensional space are straight 
lines. 

26. Deduce the differential equations of the geodesics in three di¬ 
mensional space, in terms of cylindrical coordinates. 

27. Deduce the law of transformation of the Christoffel symbols 
of the first kind in the form 


[ pq, = 


2z m 2z n dz! 3 2 zT dz n 

2z' p 2z' q 2z' s {mn ’ rJ + §mn dz p 3z' q 3z s ' 


28. Prove that 


^grrm _ q 

S U 




29. Prove that 


S_ 


(A: s B t ) 


s a:. 


Bt + A\ s 


8 B, 

S u 


30. Using the results of Problem 21, prove that 

31. If f is an invariant, use the results of Problem 30 to prove that 
the invariant g r f rs = V 2 /, where V s is the Laplacian operator discussed 
in the second last paragraph of § 48; hence show that 


v/ = — — (wg p" . 
J VgteV g g 2z r ) 
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32. If f is an invariant, use the results of Problem 31 to evaluate 
Vy in terms of cylindrical coordinates r,d,z, comparing the result 
with that in Problem 33 of Chapter IV. 

33. Prove that the permutation symbols c rs for two dimensional 
problems satisfies the identity 

Crfri = • 


34. Prove Equation (79.3). 

35. Express in cartesian tensor form the vector identity 

aX (b xc) = b(a-c)-c(a-b), 


and prove it by use of Equation (79.3). 

36. Express the eleven equations (48.1)—(48.11) in cartesian tensor 
form, and then verify them. 

37. If A rs is an absolute tensor, prove that the determinant \A rs \ is a 

relative invariant of weight -2. Establish the tensor characters of \Af. s \ 
and \A rs \. 

38. Prove that c rsl c r “ = 6, where c rst and c rst are the permutation 
symbols. 

39. Prove that ^ = 0, where g is the determinant of the metric 

oil 

tensor. 

40. If A r is a relative tensor of weight one, prove that 


a: 


i>- — 




41. If B r are the covariant components of a vector for orthogonal 
curvilinear coordinates z r whose metric form is given in Equation 
(81.2), show that the physical components of this vector for the cur¬ 
vilinear coordinates z r are as given in (81.5). 

42. If T” are the components of a tensor for the coordinates z' in. 
Problem 41, find the physical components of T ” for these coordinates 
in terms of mixed components. 

43. Let r, 0 and cp be the spherical polar coordinates of a particle. 
Prove that for these coordinates the contravariant, covariant and 
physical components of velocity are respectively 
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/ dr dQ dy \ /dr 2 d% , 

\Jtli' liy \ Ji r Jt’ r 


sin 2 


dt / 


/dr (f0 
\dt 1 dt 


>r sin 9 

dt 


\ 


44. Let r, 6 and £ be the cylindrical coordinates of a particle. Find 
the physical components of acceleration for these coordinates, in terms 
of r, 6 , z and their time derivatives. 

45. Let A . V. be a general tensor which is a function of a parameter 
t as well as of the coordinates. Prove that, for any transformation 
independent of t, the derivatives of A with respect to t are tensors 
with the same tensor character as A'.'.'.. 

46. In the theory of electricity and magnetism, Maxwell’s equations 
for free space are 


V x H 



V 


X 



V -E = 0, V-H = 0, 

where H is the magnetic intensity, k is the dielectric constant, E is the 
electric intensity, (j. is the permeability and t is the time. Express these 
equations in tensor form for general curvilinear coordinates z r • 
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